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Résumé. — Un plasma complètement ionisé est en équilibre thermodynamique. On tente une 
mise au point du calcul des corrélations de position dans ce plasma. On utilise les équations de 
récurrence et le principe de superposition. On expose d’abord l’approximation linéaire. Dans le 
cas d’un plasma neutre homogène et isotrope l’étude est poursuivie à l’approximation suivante. 


Abstract. — This paper treats of a fully ionized plasmajin thermodynamic equilibrium. An 
attempt is made at reviewing the calculation of spatial correlations in such a plasma. The equa- 


tions of recurrence and the principle of superposition are used. 


The linear approximation is first 


treated. The next higher approximation is studied in the case of a neutral homogeneous and 


isotropic plasma. 


I. Le problème. — On considère un milieu cons- 
titué par des ions positifs et des électrons. Des 
atomes et même des molécules, et aussi des ions 
complexes, peuvent se former dans ce milieu. Nous 
supposerons que la température est assez élevée 
pour qu'il n’en soit rien. Le milieu constitue donc 
un plasma complètement ionisé. La température 
n’est néanmoins pas assez élevée pour provoquer 
Papparition de phénomènes relativistes. 

Divers champs statiques peuvent être appliqués 
au plasma : un champ de gravitation, un champ 
électrique, un champ magnétique. Il est parfois 
commode — ici le point de vue est purement théo- 
rique — de supposer que le plasma est confiné 
dans un volume déterminé par des parois appro- 
priées. L'effet de ces parois peut se traduire par un 
champ répulsif converiable. On pourrait aussi envi- 
sager un système en rotation uniforme qui serait en 
équilibre dans un système d’axes tournants. 

L’ionisation complète du plasma suppose une 
température telle qu’il y règne un rayonnement 
intense. Pour que le plasma puisse rester en équi- 
libre il est nécessaire qu'il ne perde pas d’énergie 
par rayonnement. Le plus simple est de supposer 
que le rayonnement est confiné comme le plasma 
par les paroïs dont il vient d’être question : ellessont 
parfaitement réfléchissantes. Dans ces conditions 
on admet que les échanges d’énergie entre le plasma 
et le rayonnement sont faibles et on néglige l’inter- 


action entre les particules chargées et le champ 
électromagnétique. Le traitement de l’équilibre 
ainsi défini paraît à première vue se prêter à un 
traitement non quantique. C’est ce traitement que 
nous allons tenter. Le plasma en équilibre ne 
comporte aucun courant de matière. Les conditions 
où nous nous sommes placés, les approximations 
faites sont probablement surabondantes pour 
entraîner cette conclusion. Un champ magnétique 
appliqué n’a aucun effet sur l’état du système si 
bien qu’il n’en sera pas fait mention davantage (1). 

Les forces coulombiennes qui s’exercent entre les 
particules créent des corrélations entre leurs posi- 
tions. Ces corrélations ne peuvent être calculées, 
même dans le cas d’un milieu dilué, par la formule 
élémentaire de la théorie cinétique des gaz par ce 
que les forces coulombiennes sont à longue portée. 
L'importance de ce problème a été reconnue autre- 
fois par Debye et Hückel qui en ont donné la solu- 
tion : en fait ces auteurs s’intéressaient aux cas des 
électrolytes mais le cas des plasmas n’est pas essen- 
tiellément différent. 

Je me propose d’exposer ici un traitement de la 
question qui permette d’en faire complètement le 


(2) En fait les raisonnements classiques sur la question 
(Becker [38]) négligent le rôle du champ magnétique induit 
et ne tiennent compte que du champ magnétique appliqué. 
Ce point semble demander quelques recherches supplé- 
mentaires. 
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tour. Le principe de la méthode a déjà été indiqué 
ailleurs (Yvon [36]). On se limite à un plasma 
d'hydrogène. 


IT. Notations. — Celles-ci dérivent des nota- 
tions utilisées pour le traitement d’un fluide à un 
seul constituant. Nous désignons par : 


Ni = Ni(%1, Ya 31)» 


la densité moléculaire des particules de l’espèce 2 
au point de coordonnées %,, Y1, Z, — point que 
nous désignerons par x. seul pour abréger. Nous 
désignerons de même par 


Nioiÿ — Niÿ(Ti Us 715 Lo, Vos 22); 


la densité de présence simultanée de deux parti- 
cules, l’une d’espèce z en x, l’autre d’espèce 7 en æ2. 
Une relation de symétrie est évidente : 

L7 


Mjoij — Moyji- (2). 


L'absence de toute corrélation se traduirait par 
la relation : 


Mioiÿ — lui Moje (2) 


Ces relations sont toujours correctes à grande 
distance (en toute rigueur si le système est ouvert 
au sens thermodynamique, ce que nous supposons): 
c’est une évidence physique. On envisagera éga- 
lement des densités triples, telles que : 


Ti2sik- 
Nous désignerons par 
Ge 
la charge de la particule d'espèce 1, en posant : 


CET. 


Nous désignerons par 
Fr 


la composante suivant Ox de la force subie de la 
part du champ appliqué résultant par une parti- 
cule z située en x,. La force électrique, au même 
point, sera désignée par 


CG Xe 


De même la force coulombienne exercée par une 
particule 7 située en x, sur une particule & située 
en Z, S’écrira : 

G (@] 9 


avec : 


Ty — Lo 
Xe -œ" 2 1—+, 


(3) 


T12 


III. Équations de récurrence. — Les densités 
simples et multiples sont liées par les relations de 
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récurrence suivantes (nous écrirons seulement celles 
qui concernent l’axe Ox) : 


xr d 


; Ni = Fri hui + té, D 6 [rss X12 ATas (1) 
Li 0) 


d 
kT d Nyoÿ = (Fi À G CG X1o) Nroi 
Le >> GE Xis dre (2) 
k 


ae 
KT yo = (Fe; — CG GG Lie) roi 
2 


4 + 2 Ok | Nassisx Xas ts. (3) 
à \ 


T est la température absolue, À la constante de 

Boltzmann, dr un élément de volume. Les inté- 

grations portent sur tout le volume occupé par le 

plasma. 
On désignera par : 


Vis >> Niis (4) 
la densité résultante et par 
C1 = ZGnui ‘ANNE 


la densité de charge électrique. En fait il manque 
à cette grandeur, pour répondre à sa définition, un 
facteur e. Aussi vaut-il mieux l’appeler la « densité 
électrique abstraite ». 

Revenons sur les intégrales qui figurent dans les 
équations de récurrence. Considérons par exemple 
la suivante : 


Ci 2 ta [ rrsso X15 dTs 


Les points 1 et 2 étant fixés, le point 3 décrit tout 
le volume offert. Il est en général éloigné de l’un et 
l’autre des deux premiers points. C'est-à-dire que 
presque partout : ; 


Nuosiiÿ — Moi Mis 4 Jyogik = Tsij ske (6) 


Utilisons l'identité : 
u2 ta frsssie Xis dTs 
= > tu J (sosie — Nauaiÿ sk) À15 ÂTs 
SEE mois [53 13 ds (7) 
Le premier terme du second membre correspond 
maintenant à une intégration dans un petit volume, 


cependant que le second terme fait apparaître la 
force électrique moyenne sur la charge + 1 : 


A = [os Xisdrs. (8) 
Il y a intérêt à incorporer la force électrique 
moyenne dans la force appliquée, en introduisant la 


force résultante : 


Fa = Fi + GXu. (9) 


OPA ST SRE DEAR IR NS ENER IIE 


ses UE to ton. À Rue di 


lits. 


sono ie D) dé si né 


N° 10 
Les équations de récurrence s’écrivent main- 
tenant : 
(10) 


d ’ 
Nos Ni — Fi Ni st G : € f (ras — Ni 5) x dr», 


») , 
RP ni = (Fu + G G Lio) Nyois 


dÆr 


+ 2 tu J (rsasise — Naoi sk) Xis dTs, (11) 


») 2 
KT Nyoiÿ — (Foi — G C; X 12) Nioij 
T1 û 


Ft > Cu [ (rsasine — Nyoij A3) Âos ATs (12) 


Il est commode d’exprimer les corrélations à 
laide de grandeurs sans dimensions. Nous pose- 
_rons : 

Mani — Mai Noÿ aie (13) 
C’est simplement une définition de la grandeur 
sans dimensions 4:25. Nous poserons de même : 


Niogik — ii Noj Ngk Quoi (14) 


En combinant les équations (10), (11), (12) on 
fait disparaître des deux dernières la force résul- 
tante. Il vient : 


M) 
KT — ayaii = Gi C5 Ayo Arai 


dt 
+ G 2 ta froutarssise — Givi dusik) À1s ATs, (15) 
n) 
Ray = — Ci GC Ai Gros 


+ ÿ; 2 ta frsxlarssur — Gyoij os) Xos Ts. (16) 


On notera que dans ces équations la masse des 
particules n'intervient nulle part. 


IV. Approximation. Principe de superposition. 
— Nous supposerons désormais que la densité du 
plasma est modérée, quitte à préciser plus tard ce 
que nous entendons par là. Si les forces entre les 
particules étaient à courte portée il suffirait alors de 
négliger dans les équations (10), (11), (12) du para- 
graphe précédent les termes intégraux pour obtenir 
des relations satisfaisantes. Une telle simplification 
est excessive dans le cas des forces coulombiennes. 
Nous essaierons une approximation meilleure, qui 
consiste à appliquer le principe dit de superposition 
(Kirkwood). à 

Le principe de superposition — il s’agit en fait 
d’une approximation et non d’un principe — con- 
siste à admettre que les corrélations des particules 
trois à trois se déduisent des corrélations entre les 
particules deux à deux par la formule : 


Œuosiik — Auoii Gysik osik- (3) 
Cette formule paraît intuitivement assez satis- 
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faisante. Elle devient notamment exacte lorsque 
l’un des points s’éloigne des deux autres. Si c’est le 
cas du point %, Ys, Z3 par exemple ; les deux 
membres deviennent en effet égaux à 41%; simul- 
tanément. 

Compte tenu du principe de superposition, les 
équations I11(15) et III (16) s’écrivent maintenant: 


n) 
ee Gyois = Ci G Ayo 
PAIE f Givik(Gose —1) Xis dre, (1) 


d 
Ne 108 dois = — Gi GC Âu0 


+ GC 5 Cr rx Jansr(asie — 1) Xos dtg (2) 


V. Situation homogène. — Examinons ce-qui se 
passe dans une région où le fluide est : 1° homogène, 
20 isotrope, 39 neutre. 

dix; ne dépend que la distance r,,. La symétrie 
sphérique fait disparaître un certain nombre de 
termes. De plus, étant donné que la masse n’inter- 
vient pas, les électrons se comportent comme des 
ions changés de signe, ce qui conduit à admettre 
sans risquer de contradiction que : 


Gia=—= — A12++) (1) 
ordis + 
Nous désignerons pour abréger la première de ces 


quantités par a, et la seconde par b.,2. Il y a fina- 
lement à retenir seulement deux équations : 


d 
ET UT 
= Xyw +n fa Gps — Dis Dos) X13 Ts (2) 


ù 
KkT du log bo 


= — Xe + n [tes Dos — bis Gos) Aya des. (3) 


VI. Approximation linéaire. — Les équa- 
tions IV (1) et (2), qui ne sont pas linéaires, pré- 
sentent d’évidentes difficultés de résolution même 
après les simplifications du paragraphe précédent. 
Mais, comme 4, tend vers l’unité quand la dis- 
tance augmente, on imagine de poser : 


Gyoiÿ = 1 + Eyoÿ (1) 


et d’essayer l’approximation où l’on considère les # 
comme des quantités petites du premier ordre. Il 
vient, en négligeant les termes du second ordre, 


n] LA 
KT nu = Fini + GD Gifs Eysiÿ Âu2 ÂTo, (2) 
1 


dT: 
(3) 
ù + = 
kT 5x, Eyaÿ — G(C Xe + D Ck | Psk Cosik X13 Ars). 


IT Ego = G(— G Lio + > Ge fra Eysik X22 Ta). 
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De la première des équations (3) on déduit aisé- 
ment : 


ù 
da (E1245 + Exo—5) = 0, (&) 


ce qui veut dire que la somme 
En2+5 À E1o—5 


ne dépend que des coordonnées de 2, c’est une 
certaine fonction f;; : 


Eyo+5 À Eyo—ÿ = Joj, 


mais comme les e,, sont nuls, lorsque la distance 1-2 


est grande, cette fonction ne peut être que nulle. 
Donc : n 


Ey2+5 + Eyo—s = 0. 


De même, à l’aide de la dernière équations (3) on 
montre que : 


Eyoÿ+ + Eux = 0. 


Les quatre grandeurs €,,;; peuvent donc se 
déduire de l’une d’entre elles : 
POS NE PRET PUSNTRS He (5) 
GE E a En 
Nous désignerons désormais e°+4 par la nota- 
tion abrégée &,,. Le système (2)-(3) se ramène main- 
tenant aux équations suivantes : 


KT Mi = F'i Ai + üf ve €12 À39 AT», (6) 
CA 
n 
se E1o — À10 + [vs 3 19 dTs, (7) 
1 
ù ; 
LES Ego = — ÀÂy0 + fvses X23 ATa. 
2 


Il est remarquable que les corrélations, repré- 
sentées par e,,, ne dépendent que de la densité : 
elles ne dépendent pas de la charge électrique. On 
étudiera en particulier le cas où le fluide est homo- 
gène et isotrope. e;, ne dépend que de la dis- 
tance r3,+ Le système (7), intégré membre à 
membre, donne : 


Cette équation a une solution qui s'écrit, en 
supprimant les indices, 


£ = AN O Eee 


(9) 


e? kTy e? 
PER ne AN TE TT 


A titre d'indication calculons ces paramètres r, 
et « pour un plasma tel que (?) : 


avec : 


(10) 


PONTOISE 
20. Ti= 105 °K. | 
(2) Rappel: e—4,80.10—1cgs, k—1,38.1071,6 erg deg, 
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Nous donnons au lieu de « la longueur, dite 
longueur de Debye, À — «1 et aussi la distance 
qui définit la densité : 


d = v—uWs8, 


Les résultats sont les suivants : 


cm 
ro 1,66.10—11 
MORT AD TES 
h 


__2,18.10—% 


Les corrélations s'étendent assez loin dans le 
fluide, à des distances de l’ordre de h, cependant 
qu'elles ne deviennent importantes — et alors 
l’approximation n’est plus valable — qu’à des dis- 
tances très courtes, de l’ordre de rs. 


VII. Corrélations à faible distance. — D’après V 
nous pOsOns ICI : 


Ayo+ + — Big = Ayo = À + Es 


Ua + = Ajo+— = Dis = 1 + N12: 
La formule obtenue en milieu homogène 
ai2 = 1 — (rofr) er, (1) 


est justifiée aux grandes distances. Aux petites dis- 
tances elle conduit à une densité simultanée néga- 
tive, ce qui est un non sens. Imaginons un instant 
qu'aux faibles distances les effets coulombiens dus 
aux actions ternaires s’cffacent ; deux ions très 
voisins réagissent l’un eur l’autre sans trop s’oc- 
cuper de ce qui se passe autour d’eux. Il est alors 
naturel d'appliquer à la densité simultanée la for- 
mule qui convient pour les forces à courte portée, 
dans le cas des faibles densités : 


da = EXP — rofr. (2) 


Il n’est pas difficile de trouver, empiriquement, 
une formule qui couvre (1) et (2). C’est : 


die = eXp (— (rofr) e%r), = (3) 


Le même raisonnement, appliqué à la corrélation 
ion-électron, conduit à la formule suivante : 


b,9 — exp [(rafr) ea]. (4) 


Mais cette formule est inacceptable par ce qu’elle 
prévoit une accumulation infinie d'électrons autour 
d’un ion positif. On pourrait rejeter entièrement 
une formule établie de manière purement intuitive. 
Mais la singularité qui figure dans notre dernier 
résultat est inévitable : c’est une conséquence de la 
forme de l'énergie potentielle. Seul un traitement 
quantique pourrait nous en débarrasser. 

Essayons d’autre part d'améliorer l’approxi- 
mation du paragraphe précédent. Revenons aux 


four fat 


PRE CC 


di iiememet 


st teinte. > ane nés 
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équations V (1) et (2) que nous récrivons en les 
combinant : 


Del set +) 


d: 
= nf (15 — M13) (Eos + us) As dTs, (5) 


= 2X% + n fe Æ Es Æ M13) (Eos — us) Vis Ars (6) 


Ce système peut être résolu par approximations 
successives en considérant X,, (mais non pas X,) 
comme du premier ordre — ou, plus précisément, 
en négligeant 


e2]kT devant n(e?[kT) dr3. 
La solution d’ordre zéro est : 
Ey2 = Mo = 0. (7) 


Au premier ordre on retrouve les résultats déjà 
obtenus. Au second ordre il vient : 


il 4 
€12 me Éa + Mie Treo T2 = 0, (8) 


n) 1 il 
kT dt, (12 mou no E12 Se 9 ) 


= 2X 3 + n [ e + €13 + Nus) (E23— 123) À15 dTs. (9) 


On remplace dans les quantités du second ordre 
les € par leur expression du premier ordre. Il vient 
d’abord : 


32 + yo = (rofr)* er. (10) 
Ensuite on trouve que €,» — 112 à la même 
expression qu’au premier ordre. Ainsi, au troisième 
ordre près : 
PU Ne 2 
E12 Me er +; () CES 
(11) 


CE 

Ces résultats sont conformes aux formules empi- 
riques (3) et (4), ce qu’on peut considérer comme 
satisfaisant, 

Ïl serait assez difficile de passer à l’approximation 
du troisième ordre : mais comme le principe de 
superposition, sur lequel sont fondées nos équations 
de départ, n’est pas rigoureux, une telle entreprise 
n’est peut-être pas indispensable. 


VIII. Situation inhomogène. — Nous désirons 
maintenant savoir ce que deviennent les corré- 
lations lorsque le plasma n’est pas homogène. Pour 
ce problème plus difficile nous serons moins exi- 
geants et nous nous contenterons de l’approxi- 
mation linéaire du principe de superposition. On 
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revient donc aux équations VI (7), qui s’intègrent 
et donnent : 


62 e2 
KT eg = — 2 — [ve es des, (1) 
12 13 
ra 1 
[= Va Eps dre =f— Here (2) 
13 23 


On notera d’abord qu'aux grandes valeurs de 
la distance r,, le paramètre de corrélation e,, est 
nul. Il en résulte d’après (7) la condition : 


1 1 s 
_—— + [= V3 ©o3 dtr3 = 0 pour r,, grand. (3) 
le F13 


ce qui peut encore s’énoncer de la manière suivante: 
19 par la relation : 


1 + fr €12 dT2 = 0; (4) 
20 par le fait que 


CA 
: rat ? Vo ©1o ro (5) 


est orthogonal aux harmoniques sphériques centrés 
sur æ, d'ordre { — sauf pour { = 0 (il n’est pas 
facile de dire jusqu’à quel ordre en / cette condition 
doit être retenue). 

Ceci dit, si on admet que la densité globale varie 
lentement de point en point de manière que seul 
son gradient détermine les corrélations, l’expres- 
sion de €,, peut se mettre sous la forme : 


€ 


= __700—14/2rnmi +0), 


2 
12 r 


[X. Répartition macroscopique. — Nous 
sommes en état maintenant d'appliquer l’équa- 
tion VI (6) à la répartition des particules dans un 
champ de forces. A l’approximation du paragraphe 
précédent et compte tenu du peu d'influence des 
corrélations aux distances de l’ordre de r, on 
obtient les relations suivantes : 


o 
es kT —=n4x ue no FE, 


dx 3 dt 

d 2 ù () 
fr] 2 œ ni 

x. RENET 3 n— e = eV 


On en déduit l’équation d'équilibre entre pres- 
sions et forces : 


Ù : VE ‘ ) 

Le 2 = v9/2 | — F Érb=,e2 

D(var 4e PA n+ Fr +n (2) 
ce qui signifie qu'à la pression propre aux gaz 
parfaits s’ajoute une pression due aux interactions 
coulombiennes. C’est une pression scalaire qui se 


formule à 
3 Tiers 
Le VE d'A (3) 
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Étant donné que cette quantité est négative, il 
est d’ailleurs plus physique de dire que les actions 
coulombiennes créent une tension. Il est inutile d’y 
insister car, dans le domaine physique qui nous 
intéresse, cette tension est négligeable (#). Dési- 
gnons par W, l’énergie potentielle d’où dérive la 
force résultante. Les équations (1) s’intègrent : 


AT log RE = —W4 +50, 
; (a) 
RE 2 

KT ee LA 


Les constantes d'intégration C} et C- laissent à 
notre gré le choix de l’origine des énergies poten- 
tielles. 

Nous nous intéresserons aux situations où la 
densité électrique abstraite est faible devant la 
densité globale. Les équations précédentes con- 
duisent à l’approximation suivante : 


y G\ LA 2 À 
KT (log + ©) = W4 +Seta : 

Ÿ G\ICE 2 : 

KT (lg —©) = —W= +30 œ 
en continuant il vient : 
V 1 

KT 108 ———© = -(W} + W_) Lea, 6 
y a 5 (W+ ) (6) 

ae V + —W CE 
Lo} KT 2 —vlg =. (7) 


L'énergie potentielle électrique disparaît dans 

pe ire 

la combinaison W, + W_. SET 

l'unité, on en profite pour remplacer cette quantité 
par 


x est petit devant 


1e? 

log (1 + 3 er). 

(5) En astrophysique, température plus basse, densité 
plus élevée, la situation est un peu différente. 
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La répartition de la densité globale est donc 
donnée par l'équation : 


1-5 @)"-2ee 


Au laboratoire on ne saurait obtenir que de 
faibles variations de la densité globale. 
Supposons donc celle-ci uniforme. Séparons 


dans W, ou W_ l'énergie de gravitation de l'énergie | 


électrique : 


W__=Wo-— eV. 


Étant donné que : 
(10) 


AVI Tec, 
l'application du Laplacien aux deux membres 
de la relation (7) nous conduit au résultat suivant : 


Na Wot —Wo- 


D'ET RE Mi) 


Aot== + &? 0. 

Cette relation montre comment un champ gravi- 
fique sépare les charges électriques d’un système 
neutre. Elle donne aussi la loi de décroissance des 
charges dans un système purement électrostatique 
et elle permet d'estimer comment un plasma forme 
écran aux champs électriques qui lui sont appliqués. 


X. Conclusion. — Cette étude est en soi acadé- 


mique, étant donné le faible effet des corrélations. 


dans un plasma très chaud et modérément dense en 
équilibre thermodynamique. Toutefois elle pré- 
sente un certain intérêt si on veut bien la consi- 
dérer comme une introduction à l’étude. des 
plasmasihors d'équilibre. 


Manuscrit reçu le 29 mai 1958. 
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CONTRIBUTION DE LA PARTIE IMAGINAIRE DE L’AMPLITUDE DE L'EFFET DELBRUCK 


A LA DIFFUSION ÉLASTIQUE DES PHOTONS 


Par M. Pau KESSLER, 


Laboratoire de Physique Atomique et Moléculaire du Collège de France. 


Résumé. — Nous avons établi l'expression de la partie imaginaire de l’amplitude de l'effet 
Delbruck (diffusion élastique d’un photon par le champ électrique des noyaux) sous forme de 
deux intégrales quintuples, correspondant aux deux états possibles de polarisation du photon. 
Un calcul approché de ces intégrales a ensuite été effectué pour une énergie du photon de 2,62 MeV 
et quatre valeurs de l’angle de diffusion, à l’aide d’un ordinateur IBM 650. 


Abstract. — We have calculated an expression for the imaginary part of the amplitude of ‘the 
Delbruck effect (elastic scattering of a photon by the electric field of nuclei). This expression is 
given in the form of two five-fold integrals, corresponding to the two possible polarization states 
of the photon. An approximate calculation of these integrals has then been made for a photon 
energy of 2,62 MeV and four different values of the scattering angle, with the help of an IBM 650 
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computing machine. 


Introduction. — Des expériences récentes sur la 
diffusion élastique des rayons y par les noyaux [1] 
paraissent permettre une mise en évidence de 
effet Delbruck. Cet effet, prédit depuis 1933 [2] 
correspond à la diffusion des y par le champ élec- 
trique des noyaux ; ou de façon plus précise, à une 
diffusion photon par photon, où l’un est le photon 
réel et l’autre le photon virtuel associé au champ 
coulombien. 

Le calcul de ce processus du quatrième ordre, 
correspondant à l’ensemble des trois diagrammes 
de la figure 1, est extrêmement compliqué et n’a 


Prome 


encore été entrepris que pour deux cas particuliers : 
le cas de la diffusion en avant (angle de diffusion 
nul) [3], et le cas de la diffusion à très haute énergie 
sous des angles très faibles [4]. 
Cependant, il nous est apparu, à la lumière des 
résultats des expériences citées [1] qu’un calcul 
limité à la partie imaginaire de l’amplitude de 
l'effet Delbruck, pouvait éventuellement suffire à 
nous apporter les informations recherchées. En 
effet, dans ces expériences portant surides y 
de 2,62 MeV diffusés’élastiquement dans le plomb, 
les sections efficaces trouvées à différents angles se 


situent systématiquement au-dessus de la courbe 
théorique correspondant à l'addition des effets 
Thomson (diffusion par le noyau, limitée au second 
ordre) et Rayleigh (diffusion par les électrons péri- 
phériques). Or, si l’on se base sur les calculs faits 
pour la diffusion en avant, on peut s’attendre à ce 
que la partie réelle de l’amplitude de l’effet 
Delbruck donne plutôt une contribution négative 
à la section totale (par interférence destructive 
avec les effets Rayleigh et Thomson), la partie 
imaginaire au contraire une contribution positive. 
Ces considérations justifient l’intérêt que nous 
portons à la partie imaginaire. Comme nous allons 
le montrer, le calcul-de celle-ci nous est facilité par 
l'emploi d’une méthode qui associe l'effet Delbruck 
au phénomène de production de paires. 


Exposé de la méthode. — Notre procédé consiste 
à établir une formule analogue au « théorème 
optique » (relation de Bohr-Peierls-Placzek), qui 
met en rapport, d’une manière classique, la partie 
imaginaire de l’amplitude de diffusion vers l’avant 
d’une onde monochromatique dans un milieu dif- 
fuseur quelconque, avec la section efficace totale 
pour le phénomène d’absorption associé : 


{ ao, 0) }o=0 = 7 09b{w), 


où « est la fréquence de l’onde, 8 l’angle de dif- 
fusion ; notons que la formule n’est valable sous 
cette forme que pour of(o) € 06°%%(«)). 

Ici le processus d’absorption associé à la dif- 
fusion Delbruck est la création de paires électron- 
positron dans le champ coulombien. La formule 
ci-dessus a été effectivement appliquée pour obtenir 
la partie imaginaire de l’amplitude de l'effet 
Delbruck vers l'avant [3]. 
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Nous cherchons à généraliser le théorème pour un 
angle de diffusion quelconque. Pour cela, nous 
écrivons, utilisant la formule bien connue du calcul 
des perturbations au second ordre : 

-> > 
5 Ar) Hex(%) 


HKK) = 
: 8,8” k — E — E" 


p° dp dQ p’?dp a@/. 

Dans cette formule, A? et HP sont respecti- 
vement les éléments de matrice de l’hamiltonien 
d'interaction correspondant à l'effet Delbruck et à 
la production d’une paire (H?* correspond à l’anni- 
hilation de la paire). Les quantités (E, p, Q, s) et 
(E', p', Q, s') définissent l'énergie, le moment, 
l’angle solide et le spin de l’électron et du positron 
respectivement ; les indices 1 et f se rapportent à 
l’état initial et à l’état final du photon d’énergie k. 
Les AH? sont évidemment fonctions des diverses 
variables caractéristiques de l’électron et du posi- 
tron. 

Le terme du second membre représente une 
sommation sur tous les états intermédiaires tels 
que le photon incident se trouve absorbé, et que 
sont présents un électron et un positron d’énergies 
quelconques (fig. 2a). (Remarquons que dans ce 


; État. AA AR 
intermédee. e” et 


FIG. 2. 


formalisme l'énergie totale n’est pas conservée à 


l’état intermédiaire, par contre les relations énergie- 
moment sont conservées pour toutes les particules.) 

En toute rigueur, il faut ‘ajouter au second 
membre un deuxième terme, correspondant au 
processus où l’émission du photon final précède 
l'absorption du photon initial (fig. 2b). Pour plus 
de simplicité, nous n’écrivons pas ce deuxième 
terme ; nous nous contentons de remarquer qu’il 

> > 


se déduit du premier en permutant 4, et k; et en. 


changeant k en — k au dénominateur. | 
Nous allons montrer maintenant que la quantité 


1(E, E') = [EX ae(k) 44) ao ao 
8.8’ 


est purement réelle pour toute valeur de E et E. 
Pour la simplicité de l'écriture, nous grouperons 
"les coordonnées de spin et les coordonnées angu- 
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laires de l’électron et du positron sous un indice 
unique }, et nous écrirons de façon schématique : 


IE, E) = 2 H(, ÿ) HÂ(L 


Il suffit de montrer que : 1 = 1*. 
Pour cela, nous défimissons un plan # passant par 
le centre de l’atome diffuseur et parallèle aux 


> > + > | 
vecteurs k, + k; et k;, x k. L’atome, présentant 
une symétrie sphérique, a évidemment % pour plan 
de symétrie. Ceci permet d’écrire : 
Hi, D H%(, à) = HU 5) H*(,)) 

où i et f définissent le symétrie par rapport à ® de 
toutes les coordonnées (y compris la polarisation) 
du photon initial et du photon final, j la symétrie 
par rapport à ® de toutes les coordonnées (spins 
compris) de la paire électron-positron désignée 
par j. 

Tout ceci est représenté schématiquement par 
la figure 3 (en projection sur le plan de diffusion). 


Pics) 


Comme nous le montrons plus loin, il faut dis- 
tinguer deux cas : 4) la polarisation du photon 
initial et celle du photon final sont parallèles au 
plan de diffusion ; b) la polarisation initiale et 
finale est perpendiculaire au plan de diffusion, On 
remarque alors immédiatement (4 étant une fonc- 
tion linéaire de la polarisation) : 


HG, ) = + Hi); HGj) = + HG,)j), 


le signe + ou — correspondant respectivement à la 
polarisation perpendiculaire ou parallèle. D'où : 


Hi, j) AX(f, j) = A(f, j) H*(i, j) 
et 
I\E, E’) — : Hi, j) H%(f, j) = Z A(f, j) H*(i, ÿ) 
j 
= 2 H(f, j) HX(, j) = IE, E. 
Il en résulte que l’on peut partager A? en une 


partie purement réelle et une partie purement 
imaginaire, comme suit (à condition s’ajouter un 
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terme infiniment petit — 1e au dénominateur sous 
le signe somme) : 


HR) = #f | : ILE, E”) p° dp p? dp” 
re \/G, if Heat , P P P P 


HErT) G Ta 3 
T'ATES n [ak—£—#" ÎLE, E!) p° dp p?dp. 


Le terme correspondant à la figure 2b vient 
s'ajouter à la partie purement réelle. 
._ D'où l'expression suivante pour la partie ima- 
ginaire de l'amplitude de l'effet Delbruck : 
D > 
 aPa(ki, ke 


fau—r—r > 4e(k) He+() pi pp ap dQ do’. 
8,8 


ML 


Notons que nous pouvons donner une démons- 
| tration différente, basée sur l’unitarité de la 
matrice $, pour établir cette expression. 

En définissant la «matrice de transition » T par : 
S — 1 + 1T, l’unitarité de S a pour conséquence 
la relation : 
ALT TT 0 
D'où : 
1 


Im(T) = —+ (T— T*) = 5 TT* 


et si l’on considère un élément de matrice Tir : 
1 
Im(Ti) =; 2 Ti5(T*)jr. 


Si Ti représente la diffusion Delbruck, et si l’on 
se borne au quatrième ordre, il est évident que Ti; 
et (T*);, correspondent à la production et à l’annihi- 
lation d’une paire électron-positron. Compte tenu 
de : 


5 
Tis = 2Tr0(E: 2 Eÿ) HD ki, ki) ; 
— 
Ty = 2x0(E; — Ej) HPÛK) : 


(Ty = 2702 — E9 HP4(Rr), 


on retrouve bien la formule obtenue plus haut pour 
la partie imaginaire, 


PS 

Caleul de a? im. — L'élément de matrice HP) se 

compose de deux termes correspondant aux deux 
diagrammes de Feyman indiqués par la figure 4, 


Fic, 4. 
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Une remarque s’impose avant d’écrire ces termes : 
Un calcul exact nécessiterait l'emploi de fonctions 
d’onde coulombiennes pour décrire l’électron et le 
positron ; mais cela nous entrainerait à des compli- 
cations mathématiques très considérables. Nous 
utiliserons donc des ondes planes, ce qui revient à 
nous contenter de l’approximation de Born. 
À priori, ilest vrai que cette approximation n’est 
pas très légitime pour le problème considéré 
où Za« 1. Cependant, nous pouvons admettre 
que la méthode de Born nous fournira tout de 
même un bon ordre de grandeur ; il a été montré 
en effet qu'elle donne une approximation à 
environ 20 % près dans le calcul des sections effi- 
caces de production de paires dans la région d’éner- 
gie considérée [9]. 

Appliquant les règles bien connues de Feymann 
aux diagrammes de la figure 4, nous avons (en 
prenant À —=c —1): 


.Za? mm 1 


uk) = iZ D 
FUN) 


= u(Q$F + QP) »’, 


Où : 
m est la masse de l’électron ; 
u, v’ sont les spineurs représentant l’électron et 


le positron respectivement (u — u* y°); 


Mk —p—p'; 

NEA Y (Ki D) £ 
Qi % m? — hi pas at 
m + Y.(p — ki) 0. 
me — p— ki: |? 


Q? = Y.E& 
e, est le quadrivecteur polarisation du photon 
initial. Utilisant les relations : 
DPAEUnRE KE= 0", Ke = 0; 
d’autre part la propriété suivante des matrices y : 
(va) (vb) = — (vb) (v.a) + 2a.b, 
et enfin l’équation de Dirac : 
(Y.p+m)v =0,; al Yy.p +m) = 0, 
nous aurons : 
Qi = (2ki.p) 1 y°{(y.ki) (v.e) — 2p’.€i] 
Q = (2ki.p)—t v'{(v.ki) (v.e) + 2p.e]. 


Nous avons d’autre part : 


Zoo? m 4 De x 
P* Te D \()f et 
# (x) &TÈ EE pa Q?) U, 
où : 
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Qt = (2kr.p)—! y'[(Y.ki) (ver) — 2p’.€] 

Q? = (kr.p)t v°[(.kr) (v.e) + 2p.er]. 


En effectuant le produit HUE ) +?) et en 
sommant sur les spins de l’électron et du positron 
suivant un processus bien connu, on obtient ‘: 


—? > Z2 04 m2 1 
Z ax) ae*() “ARE Ep 
qi qt 
x Tr (Qt + QÙ) LEE 40 (08 + QD ve LEP. 
D'où : 
— Z? a 
PA ki, ke) Fe 16 = 
rec S(x — EE") p? dpp’? dp’ dQ d®, 
EE) (a) 


avec : 

X = Tr[ (Qi + Q) (v 
Tenant compte de : 

pdp = E GE, 


pm) v,(Q5 + Q7 Ÿ v{v.p + ml). 


œ = Mr ; p'dp = dE" 


on peut écrire : 


ae Z2CE ‘a ‘dE 
( qi ) ( sh 


Avant d'effectuer la trace X, il nous reste à consi- 
dérer de plus près la question de la polarisation du 
photon, et d’autre part à définir un système d’axes 
convenable. 

Pour le photon incident comme pour le photon 
diffusé, nous pouvons définir les deux états quan- 
tiques suivants : polarisation parallèle au plan de 
diffusion, polarisation perpendiculaire à ce plan. 
Des considérations élémentaires, basées sur le fait 
que l’opérateur de symétrie par rapport au plan de 
diffusion est une constante du mouvement, 
n ontrent qu'il ne peut y avoir transition d’un des 
deux états à l’autre dans la diffusion. Donc, nous 
avons à traiter les deux cas : a) la polarisation 
initiale et finale est parallèle au plan de diffusion ; 
b) la polarisation initiale et finale est perpendi- 
culaire à ce plan. 


dQ dQ’ X. 


DFA ® mpp' dE 7 
aPa(k, ki) —= 16 © °° | PCR 
« qi gt 
—> \2 ANS 
e Ca ) + ( qi ) (° ‘) (: 1) 
k? SRE RE 


n(ra+e)(ré+ rt) 
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Nous choisissons maintenant le systèmen qe 
suivants (fig. 5) : l’axe 1 est dirigé suivant s + ke 


IG 5: 


> > 
l'axe 3 suivant À; X kr, le plan (1, 2) est donc le 
plan de diffusion. Dans le cas de la polarisation 
parallèle (a), l’axe 2 constitue la bissectrice de 


Le _ 

Ei» 
+ 

culaire (b), Es == a) est porté par l'axe 3. Soit 0 le 

demi-angle de diffusion. 


Il nous sera utile de définir les quantités sui- 
vantes : 


«) : ; dans le cas de la polarisation perpendi- 


= k-LK;.p = E — p, cos 0 + p, sin 0 ; 
Ke KID EN te 1e cos 0 + Pa sin 0 ; 


sg = k1 K£.p = E — p, cos 0 — p, sin 0 ; 
st = k—ki.p = E’ — p, cos Ô — p, sin 0, 
et d’autre part les quantités : 


Ti = Ps COS 0 + p,sin 0; 


pa cos 0 + p' sin 0: 


: 
rf = Po COS 0 — p, sin 0 ; rt 


+ >> 
correspondent respectivement à p.æ, p'.æ, 

++ >> 

P.&, p'.æ dans le cas de la polarisation parallèle. , 


Le calcul de la trace nous conduit aux résultats 
suivants, après un travail assez laborieux de regrou- 
pement des termes : 

a) Polarisation parallèle : 
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Pa cos 0 — Pi sin 0, - 


pk pm A 6 


Mn AO msi 
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D 


Si se Si Sf 


+ 2 sin? 0 (2) (— 


i Sj/ \Sf 


) Lipa + 79% + (ps + pi)? — 


BR Past pari Pat Pan + Post 
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(E — E’}° sin? 6] 


4 LA 
3) (P3 Ps — m?) 


Piri rs) 


+ 2 sin o 


Si sf 


Sf si 


Don [( A) (ES PP 
SH eus : 


Sf 


= (E=*) É cos 0 + pr 
sf 


St si 


b) Polarisation perpendiculaire. 
ZEiRE 


> » 
aPa(K, ke) miss ul Ce 


gi 


V [Ca ne Cœ)° (:: 


k2 Si 


Si Sÿ 


mpp' dE 


PO 
or 


E’ eos 0 + Hi) 


St 


E" cos 0 + bi] | 


Si 


dQ dQ 


P3 


, Ps Ps Ps P3 
(#2 + EË) (res s+ 52) 


5 ++ 


Si Sf Si St 


] (pa + ps)? + (pa + ps)? — R? sin? 0) 


+ 2 sint 0 (© + =) (2 + + (Ps P3 + m?) 


S; sf 


S 1 
+ 2 sin Le ae) (Ep, — E’ P2) 


i Sf 


Ro [(2 “ (Ë cos 0 + p, E’cos 0 + &) 
Si si : St 


S$ 


St 


Rappelons que dans ces formules on a : 
E =k—E. 


On vérifie aisément, que si-dans l’une ou l’autre 


de ces expressions on fait 0 —0, et que l’on 


applique le théorème optique : 


oP(k) = _ ab CH, ki, #4) 


( a) ( cos 0 + p, 


E” cos 0 + m)] 


si Si 


1 


on obtient bien la formule de Bethe et Heitler pour 
la section efficace totale de production de paires. 

Le calcul des deux intégrales quintuples a été 
effectué pour £ — 2,62 MeV et pour quatre valeurs 
de 9, à l’aide d’un ordinateur IBM 650, en utilisant. 
la méthode dite « de Monte-Carlo ». Les résultats 
ainsi obtenus sont approchés à environ 10 %, près, 
ce qui nous à paru suffisant pour nos besoins. 

Ces résultats sont (ab, en unités Z? &? 7) : 


k 
(] 00 150 30° 150 
aPalpol. 1) 0,248 + 0,023 0,082 + 0,006 0,031 + 0,003 0,020 + 0,002 
aha(pol. 1) 0,239 + 0,031 0,058 + 0,006 0,011 + 0,003 0,000 + 0,001 


Conclusions. — Les résultats ci-dessus corres- 
pondent à des angles de diffusion de 00, 300, 60o, 
900, 

Notons que pour le cas de la diffusion en 
avant (00), où l'intégration a pu être effectuée de 
façon analytique, le résultat est exact est : 


FLE (pol. 1) — hn (pol. 1) = 0,265, 


_ 


en bon accord avec les valeurs approchées que nous 
obtenons. 

Les résultats étant régulièrement décroissants 
dans les deux cas (parallèle et perpendiculaire) à 
partir de l’angle de diffusion nul, il est facile d’inter- 
poler pour les angles intermédiaires. 

Confrontés avec l’expérience, ces résultats ne 
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comblent malheureusement pas notre attente. En 
effet, ils n’expliquent que dans la proportion de 10 
à 15 % l'écart entre la courbe expérimentale 
donnant les sections efficaces de diffusion élastique 
des y, et la courbe théorique correspondant à 
l'addition des effets Rayleigh et Thomson. Il faudra 
donc chercher l'explication principale de cet écart 
dans une voie différente qui peut être l’une des 
suivantes : 

19 Un calcul plus exact de l'effet Rayleigh 
à 2,62 MeV, les valeurs utilisées jusqu'à présent 
provenant d’une extrapolation sujette à caution. 
Un tel calcul est actuellement en cours d'exécution. 

20 Une contribution positive de la partie réelle 
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de l’amplitude de l'effet Delbruck (contrairement 
à ce que nous pensions @ priori). Notons à ce propos 
que le fait de connaître la partie imaginaire autorise 
l'espoir de pouvoir calculer la partie réelle grâce à 
l'emploi de formules de dispersion. 

Je tiens à remercier M. Philippe Eberhard qui 
m'a considérablement aidé à mettre au point la 
méthode de calcul que j’ai utilisée dans le présent 
travail. 

J’exprime également ma reconnaissance à 
MM. Feix et Nicourd et Mme Sajaloli, qui se sont 
chargés du calcul à la machine IBM 650. 


— Manuscrit reçu le 11 juin 1958.. 
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REVUE DE LIVRES 


Riçcaz (R.), Les hyperfréquences, circuits et propagation 
des ondes. 32 édition. (4 vol., 16 X 25 cm, 256 pages 
175 fig., Éditions Eyrolles, Paris, 4957, 2 000 F.) 
L’auteur prend soin, au début du livre, de rappeler, 

sans s’encombrer de démonstrations supposées connues, ce 
qu’il faut déjà savoir pour le comprendre : les équations 
de Maxwell, le théorème de Poynting, les résultats qui 
résument l’étude des lignes de transmission (diagramme du 
cercle). 

Les chapitres se succèdent logiquement et à l’intérieur 
de chacun on va du plus simple au plus compliqué. C’est 
ainsi que le chapitre IT a pour objet, avec la notion essen- 
tielle d’impédance d'onde, l’étude d’ondes tombant norma- 
lement d’abord, obliquement ensuite; sur la surface de sépa- 
ration de diélectriques purs ou sur la surface de conducteurs 
parfaits, puis l’étude dans les mêmes cas des ondes électro- 
magnétiques. (polarisation horizontale, verticale). 

De l’étude au chapitre [IT des ondes guidées dans des 
milieux idéaux (diélectriques purs, conducteurs parfaits, 
on passe aux cas pratiques où il y a des pertes (chap. IV) 
puis aux particularités provenant de l’introduction d’obs- 
tacles ou de discontinuités dans les guides. Les cha- 
pitres VIE, VIII et IX comprennent chacun une vingtaine 
de pages. L'auteur y montre comment, à.partir des réso- 
. nateurs classiques, on arrive aux volumes résonnants utilisés 

en hyperfréquences. Les systèmes rayonnants diffèrent 
aussi des systèmes classiques. On étudie (chap. VIII) les 

‘cornets, miroirs, fentes. Puis la propagation des hyper- 

fréquences dans l’espace (chap. IX). 

Trois paragraphes constituent le dernier chapitre. Ils 
ont pour titre : utilisation d’obstacles en ligne pour amé- 
liorer le taux d’ondes stationnaires, réalisation de filtres, 
emploi des ferrites en hyperfréquences. 

née photos illustrent, ‘à la fin du volume, certaines 
réalisations, 


Ce livre fait partie de la collection du C. N. E. T. et 
s’adresse aux étudiants ou ingénieurs, ayant déjà une 
bonne formation en électricité, et qui veulent aborder 
Pétude des hyperfréquences. Cet ouvrage est clair. L’auteur 
part de phénomènes simples et fait souvent appel à l’intui- 
tion. Il évite ainsi des développements mathématiques 
compliqués. 


E. DIEULESAINT. 


CARRIÈRE (P.), Méthodes théoriques d’étude des écou- 
lements supersoniques. (1 vol. broché, 18,5 X 27 cm, 
de 234 p. Publications scientifiques et techniques du 
ministère de l’Air, Paris, 1957. Prix : 3 000 F) 

Un excellent ouvrage, d’une grande clarté, illustré de 
nombreuses figures. Cet ouvrage destiné aux ingénieurs a 
pour but de les aider effectivement dans le calcul des 
tuyères, des écoulements et contient à cet effet des tables 
numériques : Tables des écoulements permanents avec et 
sans choc, pour diverses valeurs de y = CJc; Tables des 
écoulements supersoniques par la méthode de Busemann. 
L'ouvrage comprend : Chapitre 1€r : Théorie monodimen- 


sionnelle des écoulements permanents et théorie des ondes | 


de choc. Chapitre 22 : Théorie élémentaire des écoulements 
plans supersoniques. Chapitre 3° : Exposé de la méthode 
des caractéristiques et applications simples. Chapitre 4e : 
Application de la méthode des caractéristiques au caleul 
des écoulements'de révolution: 

Nous recommandons vivement la lecture de cet ouvrage, 
qui, orienté vers la pratique, n’en contient pas moins un 
exposé clair de la théorie des caractéristiques et est d’un 
niveau mathématique élevé. : 


J, WINTER, 
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SUR L'EFFET BLOCH-SIEGERT ET L'EFFET D'UN CHAMP ALTERNATIF PERTURBATEUR 
EN RÉSONANCE MAGNÉTIQUE ET QUADRUPOLAIRE 


Par François LURÇAT, 


Laboratoire d’Électronique et de Radioélectricité de la Faculté des Sciences de Paris (Fontenay-aux-Roses). 


Résumé. — On indique une méthode qui permet d'étudier le déplacement de la fréquence de 
résonance dû à l’amplitude finie du champ radiofréquence, ou à la présence d’un autre champ de 
fréquence différente, à partir des équations macroscopiques. Cette méthode est utilisée pour cal- 
culer ces déplacements dans le cas de la résonance quadrupolaire de noyaux de spin 1, 3/2 et 5/2. 


Abstract. — À method is given, wich allows the study of the resonance frequency shift due to 
finite radiofrequency field amplitude, or to another field with a different frequency — starting 


from macroscopic equations. 


This method is applied to the computation of these shifts, in the case 


of quadrupole resonance of nuclei with spins 1, 3/2 and 5/2. 


LA 


1. Introduction. —" Bloch et” Siegert” [1] ont 
montré qu’en résonance magnétique, pour des 
atomes de spin 1/2, la composante circulaire du 
champ magnétique ‘alternatif qui tourne en sens 
inverse de la précession de Larmor produit un 
déplacement de la fréquence de résonance. Celle-ci 
n’est donc pas exactement ©, = |y| 44, mais 


o% = wf1 + (2) (1) 


S0 


où H, est l'amplitude du champ radiofréquence 
linéaire (l'amplitude de chacune des composantes 
circulaires est donc /7,/2). Cet effet n’est pas déce- 
lable expérimentalement dans les conditions 
“usuelles. Mais Seiden [2] a montré qu’un champ 
perturbateur ,/2 tournant à la vitesse w,, voisine, 
mais cependant suffisamment éloignée de la fré- 
quence de résonance, produit un déplacement 
beaucoup plus grand, donné par la formule : 


ESS 1 (A) 
ARC tee (2) 
(on retrouve la formule (1) pour ©, = — «,, 


c’est-à-dire pour un champ tournant en sens inverse 
de la précession de Larmor). Ce même résultat a été 
ensuite obtenu par Ramsey [3]. F. Bloch [4] a 
indiqué qu'il se produit, non seulement un dépla- 
cement, mais encore un dédoublement de la raie de 
résonance, sous l’action d’un champ de fréquence 
voisine (mais distincte) dé la fréquence de réso- 
nance, et Anderson [5] a observé ce phénomène. 
L'effet Bloch-Siegert et l’effet d’un champ de fré- 
quence différente de la fréquence de résonance ont 
été observés récemment par plusieurs auteurs [6], 
[7]. Winter [8], [9] a observé et calculé cet effet 
dans le cas de transitions atomiques à plusieurs 
quanta. 

Au point de vue théorique, les auteurs cités ont 
calculé ces effets, soit par une méthode quan- 


tique [1], [31, [4], [8], soit par une méthode”géomé- 
trique [2] Nous donnons ici une méthode qui 
permet de calculer cet effet directement à partir des 
équations macroscopiques auxquelles obéit le phé- 
nomène de résonance (équations de Bloch dans le 
cas de la résonance magnétique). 


2. Méthode de calcul. Cas de la résonance 
magnétique. — Ecrivons les équations de Bloch [10] 
dans le système tournant à la vitesse angulaire 
quelconque «, (on déterminera ensuite «, de 
façon à ce qu’un champ tournant à cette vitesse 
produise la résonance) : 


My = u CoS wyt— v sin wt 


My = F (usin w, t + v cos œt) i 


et en présence d’un champ perturbateur de fré- 
quence ©, dont l’amplitude sera dans ce para- 
graphe désignée par 2/7, : 


Hay = 2H; cos w;t RENE, (4) 
On obtient : 
du 
Fa c— Bu— Dv — M, 
sin (@2 + 3) t + sin (os — w;)t 1] 
ARTS NE NON 
o (5) 
{ cos (@» ++ &3) t + cos (os — w;)t }] 
= e[u | sin (oo + ©) t + Sin (© — w,)t 
+ v cos (@2 + &3) t + cos (2 — w;)t | AM, — M)] 
1 
ÉRarUPT} PRE 
ME Ft Ts 
LES 1 Se [xl Ho — © 
| lvl Hi To lvl A 


Nous utiliserons pour intégrer ces équations la 
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«méthode de la moyenne » exposée par Bogolioubov 
et Mitropolski [11]. Cette méthode consiste à déve- 
lopper suivant les puissances de € la solution du 
système différentiel 


DE Mules 2 Da) Les AN VUS 4702 NO) 
où € est un petit paramètre, et où les X, sont de la 
forme : 


een ol eu NE ne IR 
V 
les v étant des fréquences quelconques. 
La première approximation de la solution de (6) 
est la solution de 
dax 


d CAGE ANNEES) (8) 


où X, désigne la moyenne dans le temps de À. Ce 
fait est bien connu en résonance magnétique, 
puisque pour les équations de Bloch usuelles, 
l’approximation (8) revient simplement à négliger 
la composante circulaire du champ qui tourne en 
sens inverse de la précession de Larmor. 

La deuxième approximation de la solution de (6) 
est: 


TR — Ex + € Xr(t, Er, …. En) (9) 


où 
* De. 
DAME ER CN EE a eR 
xl ? G1» En) v£0 IN kv(éx ë ) 


et où &, est la solution du système 


dé ——— ———— 
= SAN En... 6) 
! 2 T [a a ù es (4 
er DIN CE Es) 57 XE(tET een LU) 
q 4 


La correction e X, dans (9) n’a pas d'intérêt ici. 
Elle représente en effet de petites composantes des 
grandeurs w, v, M,, composantes qu’on n’observe 
pas par les procédés expérimentaux usuels parce 
qu’elles ont une fréquence différente de la fré- 
quence de résonance. ; 

Par contre le système (10) donne directement 
l'effet du champ perturbateur. Il s’écrit en effet, 
pour les équations (5) : 


du 


FPE e-— Bu—3v] 
D = up (11) 
= ef a(Me— M] 
où 
helene at en h (12) 


N°10 


Les équations (11) sont exactement les équations 
de Bloch ordinaires, sans champ alternatif, mais 
où Gest remplacé par 9’. Quand on introduira un 
champ magnétique tournant à la même vitesse w, 
que le système tournant, on trouvera donc que 
ce champ produit la résonance pour à = 0, c’est- 
à-dire pour : 


4. AN PANPRETSSS 
@9 + O9 — 


, 1 
& = we = oo +3 (r Hi)° ( )- 43) 
Si nous appelons maintenant l’amplitude du 
champ alternatif ,, comme dans le reste de cet 
article, la fréquence de résonance s’écrit : 


1 1 
O9 + 1 


‘ Las 4 2 
@% = @o + Z (YA) ( ne (14) 

Ce résultat est en accord avec la formule (2). Il 
est valable pour la résonance en présence d’un 
champ  perturbateur : FT, = core 
Hy = H, — 0, dont la fréquence est suffisamment 
éloignée de «w, pour qu’il n’excite pratiquement pas 
de transitions. Pratiquement, seul le second terme 
est important : 


= 0 + (TA (15) 

Les deux termes de la formule (14) proviennent 
des deux composantes circulaires du champ per- 
turbateur. Si maintenant on étudie la résonance 
sans champ perturbateur, avec un seul champ 
alternatif rectiligne, de fréquence voisine de la réso- 
nance, on peut voir directement que le premier 
terme disparaît, et on obtient avec le second seu- 
lement la formule (1). 


3. Résonance quadrupolaire de noyaux de spin 
1 — 1. — Écrivons les équations de la résonance 
quadrupolaire [12], transformées par 


% ñ : 
M, — Mÿ cos out — © 0% sin ot 
(16) 
À LE YR % 
Len Que = Me Sin ot + a Quz COS wat 


et en présence d’un champ magnétique alternatif 


He == PHONE NH (17) 


perpendiculaire à l’axe de révolution du gradient 
de champ électrique. On obtient : 


ge = — (oo — 0) 0% —27 [08 — 0h] 


[sin (©: + ©) t + sin (os — &,)t] 


RL 0%) = (o— 0) ME —24 1, 4 (08 — 0) 


[cos (w2 + ©) t + cos (os — w,)t] (18) 


PO ST PE 


= 19 [mi {sin (os + @j) 1 sin (0: — et} 
+ 1 Qù{ cos (@1 + &2) t + cos (o— 0,)1}|. 


- En écrivant l’approximation (10), nous obtenons 

les équations de la résonance quadrupolaire ([12], 
équations (27-I) avec À, = 0), où la fréquence de 
résonance &, est remplacée par 


: + —— ). (19) 


O9 — @1 


6 = O0 + : DA 
Bien entendu, comme en résonance magnétique, 
il y a dédoublement de la raie de résonance, et 
_ l’autre composante se trouve à la fréquence 
| 201 — : 

Si le champ alternatif n’est pas perpendiculaire à 
Paxe de révolution du gradient de champ élec- 
trique, mais fait avec cet axe un angle 6, y 7, doit 
être remplacé par y À, sin 0. 

Si on étudie la résonance sans champ pertur- 
bateur, avec un champ rectiligne de fréquence 
voisine de la fréquence de résonance, un caleul 
direct montre que le terme en 1 /(o,— «w,) dans (19) 
disparaît, et on obtient l’effet Bloch-Siegert pro- 
prement dit, analogue à (1) : 

Op — ce [1 +7 (=) | 


20 
me (20) 

Mais il faut observer que ce déplacement n’est 
pas dû, comme pour la résonance magnétique, à 
l’une des composantes circulaires du champ alter- 
natif rectiligne. 


4. Résonance quadrupolaire de noyaux de spin 
1 — 3/2. — Les équations de la résonance qua- 
drupolaire, transformées par la transformation 
adaptée à la fréquence «>, 


x = U COS @2 t — UV Sin ot + uw’ 


(21). 


ñ 5 
LE Que = u sin ©O9 t + U COS mot 


vie 


et en présence d’un champ 4, — H, cos «, t, 
He 0s'écrivent [12] : 


d H ae 
TS = — (@9 — ©) 0 + 37 [—w {sin (@2 + o)t 
+ sin (wo — &) t } + 2x cos wt] 
d H 
_ = (@9 — @) u + VS —w { cos (w2 + @)t 
+ cos (2 — 1) t } — 2y cos ot] 
dæ H 


= +5 lu {sin (oa + ©) t + sin (wo: — ©) t} (22) 


+ D{ cos (o, + w) t + cos (w, — w,)t}] 
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141 
dæ 1 s 
ane (@o — &2) y + 9 Y Ail(8u = — 4) 
{ cos (@2 + 1) t + cos (we — ®1) t}— 4u cos ot] 
d 1 ; 
D = —(o0— 0) 2 +7 Hu — 2) 


{ sin (©: + ©) t + sin (@2 — w:) 1} + 4o cos wi t]. 


En écrivant pour les trois premières équa- 
tions (22) l’approximation (10) nous obtenons les 
équations dela résonance quadrupolaire ([12], équa- 
tions (37) avec Æ, — 0) où la fréquence de réso- 
nance «, est remplacée par 


L 1 


, 3 
Gi 0 +5 A (— & 
0 1 


Ê (23) 


O9 — O1 


On peut faire les mêmes remarques qu’à propos 
de l’équation (19). 

Pour la résonance sans champ perturbateur, 
avec un champ de fréquence voisine de la fréquence 
de résonance, on obtient pour l’effet Bloch-Siegert 
proprement dit : 


2 
% = 9 [1 ++ (2) | 
Lorsqu'on écrit pour les trois premières équa- 
tions (22) l’approximation (10), on constate que les 
termes en x et y dans ces trois équations ne four- 
nissent aucune contribution. Les phénomènes de 
résonance quadrupolaire sont donc explicables, du 
moins à cette approximation, par les trois premières 
équations (22) sans termes en x et y, c’est-à-dire par 
les équations de Bloom, Hahn et Herzog [13]. C’est 
pourquoi l'effet Bloch-Siegert est € — 3 fois plus 
grand que pour la résonance magnétique, résultat 
qu’on obtiendrait en calculant cet effet directement 
à partir des équations de Bloom, Hahn et Herzog. 
Physiquement ceci s’explique très bien, puisque les 
équations de ces auteurs supposent que seuls inter- 
viennent, dans les phénomènes de résonance qua- 
drupolaire entre les niveaux + m, + (m — 1), les 
niveaux + m, + (m 1) entre lesquels se pro- 
duisent les transitions de résonance — ce qui est 
bien entendu rigoureusement vrai pour À — 3/2, 
puisqu'il n’y a pas d'autre niveau. 


(24) 


5. Résonance quadrupolaire de noyaux de spin 
I — 5/2. — Les équations de la résonance quadru- 
polaire, transformées par : 


My = u coS w3t — v sin w,t + u’ cos 2w,t 


— v’ sin 20,t +u” 
(25) 


ble 


ñ Û 
À Que = u sin &9t + D COS ot 


+ Qu’ sin 20, t + 20” cos 2œ, t 
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et en présence d’un champ A, —H, cos o;t, H; — H, — 0, s’écrivent [12] : 


du _ 
dog 


— 2y° coS ot + y} cos (202 + ©) t + cos (202 — w)t}] F 


+ 2x’ cos w,t + y {sin (202 +") t + sin (Co, — w,)1}] 


dæ H;, 


 - 1 Luf sin (o2 + @) t + sin (0, — w,)t}— w{sin (262 + ©) t + sin (20: — o1) 1} 


+ vf cos (os + ©) t + cos (2 — w;) 1} — v’ { cos (203 + &,) t + cos (26: — w;) t}] ; 


—5{ cos (wo + ©) t + cos (op — &;)t } + v' À cos (2w, + ©) t + cos (26 — ©;)1 | 


— 5u À sin (262 + ©.) 1 + sin (26: — &,)t } — 8v' | cos (@> + &,)t cos (os — w,)t | 
+ 5v | cos (26% + ©) t + cos (20, — w,)t } — 10h cos w;t] 


-— 5u À cos (269 + ©) 1 + cos (262 — w) 1} + 8v' | sin (@2 + ©) & + sin (we — ©) 1} 
— 5v | sin (262 + ©) t + sin (26: — w;) t} + 2a cos œit] 


(les autres équations ne sont pas utiles). 
Si nous étudions l’effet d’un champ perturbateur 
1 3 
= 5) a 5) , nous 
calculerons l’approximation (10) des trois premières 
équations (26). Nous obtenons ainsi les équa- 
tions (46) de [12], avec 7, — 0, et où o, est rem- 
placée par : 


sur la résonance de la transition | - 


1 1 
6 = 00 + (HD ( = + ———) 
O9 + O1 Op — 
(27) 
5 


= (y A? ( 


= 


1 1 
209 + ©: n 209 — _. ; 
. On voit que la fréquence de résonance w,, corres- 
: VOTE 1 3 x 
pondant à la transition (= 9 5 , peut être 


modifiée notablement, soit par un champ de fré- 
quence voisine de «9, soit par un champ de fré- 
quence voisine de l’autre fréquence de résonance, 
2,. Le premier effet est 8 fois plus grand que pour 
la résonance magnétique, ce qui est bien le résultat 
que donneraient les équations de Bloom e.a. 
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— (©9 — &2) 0 — 4y H,fw sin (wo: + &,) t + sin (os — ©) 1} + z À sin (262 + ©.) t + sin (20, — w;) 1} 


Res (9 — ©) u— 4y H,[w cos (w3 + &3) t + cos (w: — w:;) 1}— x! cos (262 + ©.) t + cos (20: — w;) 1} 


ee — (260, — 2w;) u’ 5 [5° À cos (202 +: @4) à + cos (26: — ow;) 1} + 8x | cos (@2 + &1) t + cos (2 — w:) 1} 4 
(26) 
= [Si sin (@» + ©) 1 + sin (os — w;) 1} + u' À sin (26: + ©) t + sin (26: — «w;) 1} 


= — (304 — 302) y — YH, [8u’ | sin (w2 + ©) t + sin (os — w,)t 


: | 
a (309 — 3) x Het H;[8u' | cos (wo + ©) t + cos (oo — ow,) 1 | 


= — (26, —202) y (5w” À sin (2w2 + &w:) t + sin (260: — &) t} — 8x | sin (@o + @) t + sin (2 — w;) 1} 


— 8y { cos (w + ©) t + cos (@ — 1) t}] 


— 8y {sin (o, + ©) 1 + sin (0, — w;)1}] 


(6 — 8). Par contre, le second effet ne peut évidem- 
nent pas être prévu à l’aide des équations de 
Bloom e.a. . 

Pour la résonance sans champ perturbateur, avec 
un champ de fréquence voisine de la fréquence de 
résonance &,, le terme 1/(w, — «,) dans (27) dis- 
paraît, et les trois autres termes donnent : | 


O6 — 09 [1 + 5 (21: 


Si nous étudions l’effet d’un champ perturbateur 


(28) 


De 3 5 
sur la résonance de la transition | + 5 . 5): nous 


calculerons l’approximation (10) des trois équa- 
tions (26) relatives à u', v' et w', et nous trou- 
verons une fréquence de résonance : 


| 5 «( 1 1 ) 
_ = RER A EUR «2» 
20% AE 8 (x A5) 209 + © 209 — @r 
| (29) 

1 À 1 1 
Fee (r-H3) Le + ©1  @p — ) 


La modification de la fréquence de réso- 


No 10 


| 
| 
| 
| 


&,] 


N°10 


nance transition 


3 
(+5 


de 265, est 5 fois plus grande que l’effet corres- 
dondant en résonance magnétique (€ — 5). 
Pour la résonance sans champ perturbateur, avec 


2&,, correspondant à la 


ne 3]: par un champ de fréquence voisine 


… un champ de fréquence voisine de la fréquence de 
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résonance 264, on obtient l'effet Bloch-Siegert pro- 
prement dit : 
, ATEN 
9 == tee 
26% 20, [1 nn j: 
Les phénomènes produits par un champ pertur- 

bateur, et l’effet Bloch-Siegert, se calculent de 


façon analogue pour les autres valeurs du spin 
nucléaire. 


(80) 


Manuscrit reçu le 7 mai 1958. 
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REVUE DES LIVRES 


Rogin (L.), Fonctions sphériques de Legendre et fonctions 
sphéroïdales. (1 vol., 16 X 25 cm, xxxv + 201 pages, 
Gauthier-Villars, Paris, 4957, 4 000 F.) 

Ce volume est le tome I d’un traité détaillé sur les fonc- 
tions sphériques et sphéroïdales, fonctions qui jouent un 
rôle de plus en plus considérable en Physique Mathéma- 
tique. 

Ce tome I contient la partie la plus classique du traité : 
Chap. I : « Introduction de l’équation différentielle associée 
de Legendre en Physique Mathématique. Les polynômes 
de Legendre et les fonctions de seconde espèce, d’indice 
entier, positif ou nul. » — Chap. II : « Les fonctions asso- 


ciées de Legendre, Py(u) et Qu(u), pour m et n entiers, 
n positif ou nul. » — Chap. III : « Les harmoniques sphé- 
riques de degré entier, positif, négatif ou nul. » 

Quatre brèves annexes donnent des expressions explicites 
des polynômes P,{u) et des fonctions Q,{u), pour n allant 
jusqu’à 7, ainsi que des fonctions associées Py(u) et Qu (u), 
pour mn et n allant jusqu’à 4. 

Toutefois il ne s’agit pas là d’un simple recueil de formules 
et de recettes entassées vaille que vaille à l’usage des 
physiciens et autres mathématiciens par nécessité. Bien au 
contraire : ce traité est véritablement le fruit d’une révi- 
sion minutieuse et perspicace de toutes les connaissances 
actuelles en la matière. L’auteur n’a pas hésité à donner 
toutes les démonstrations et surtout les conditions de 
validité, assurant ainsi à chacun la possibilité d'appliquer à 
bon escient chaque résultat. Chemin faisant, des lacunes 
ont été comblées par des résultats originaux et des démons- 
trations ont été renouvelées et simplifiées. En outre de 
nombreuses erreurs — de calcul ou simplement d’impres- 
sion —, demeurées sans rectifications dans les textes de 
peus ont été détectées et corrigées pour la première 

ois. 

Une copieuse « Introduction et vue d’ensemble », placée 
en tête de cet ouvrage et coiffant le traité tout entier 


(3 tomes), constitue un fil d'Ariane qui sera sans doute 
précieux à tous les utilisateurs des fonctions sphériques ou 
sphéroïdales. ‘ 

M. Jesse. 


ALLEN (J. S.), Le neutrino. 4 vol. 16 X 24 cm, 168 p., 
relié, Princeton University Press, 1958 (en anglais). 
Prix.:..$ 4,50: 

Il s’agit d’une monographie complète, avec indications de 
toutes les expériences faites pour établir l'existence et les 
propriétés de ce corpuscule doué d’énergie (sans masse au’ 
repos) et de spin seuls, et imaginé par Fermi pour sauver le 
principe de conservation de l’énergie au cours de l’émission 
des rayons bêta. Les schémas de montage et les courbes 
donnant {es résultats des mesures sont d’une grande clarté. 
Mais les explications théoriques sont très sommaires, et la 
lecture de l’ouvrage est difficile pour le physicien qui n’est 
pas spécialiste nucléaire. Il est vrai que la condensation 
d’un si grand nombre d’expériences dans un si court 
volume exclut la possibilité de dissertations théoriques. 
Chapitre I : Propriétés générales du neutrino (neutrino et 
antineutrino, observation indirecte et directe du neutrino). 
II : Masse au repos du neutrino. III : Reculs dus aux neu- 
trinos émis lors de la capture par un noyau d’un électron 
orbital, IV : Corrélation angulaire électron-neutrino dans 
la désintégration bêta (conservation de parité, invariance 
pour un renversement du temps, conjugaison de charge 
dans la désintégration bêta). V : Expériences sur la corré- 
lation angulaire électron-neutrino. VI : Double désinté- 
gration bêta. VII: Détection du neutrino libre (expé- 
riences d’absorption, moment magnétique, désintégration, 
bêta inverse). VIII : Réactions mésons-neutrinos (désinté- 
gration du méson pi, capture et spectre bêta du méson mu, 
interaction universelle de Fermi). 


J. Winter. 
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LES GROUPES CRISTALLOGRAPHIQUES MAGNÉTIQUES ET LEURS PROPRIÉTÉS 


Par Y. LE CORRE, 


Laboratoire de Physique Théorique, Collège de France, Paris. 


Résumé. — Tavger et Zaitsev [6] ont montré que les opérations de symétrie qui ramènent le 
milieu cristallin sur lui-même sont en nombre plus élevé pour un cristal antiferromagnétique que 
pour un cristal paramagnétique. Il faut en effet considérer les «antiopérateurs» qui sont le produit 
d’une opération géométrique d’espace par un retournement du.temps. Le nombrè des classes 
cristallines (groupes de symétrie d'orientation) n’est plus 32 mais 122 (Choubnikov [7])}. Parmi 
ces 122 classes, 31 sont ferromagnétiques, 59 sont nécessairement antiferromagnétiques (dont 
38 piézomagnétiques) et 32 sont soit diamagnétiques, soit paramagnétiques, soit antiferromagné- 
tiques (il existe alors deux sous-réseaux antiparallèles de spins décalés l’un par rapport à l’autre dans 
une antitranslation élémentaire). Divers tableaux décrivent les propriétés des 122 classes magné- 
tiques. 

Le nombre des groupes d’espace n’est plus 230 (Schônflies, Fedorov) mais 1 651 (Zamorzaev [8], 
Beloy, Neronova, Smirnova [9], [10]), 230 sont dia ou paramagnétiques, 272 sont ferromagnétiques, 


517 comportent deux sous-réseaux antiparalèles de spins décalés dans une antitranslation élémen- . 


taire, ils sont donc antiferromagnétiques ; les 632 groupes restants sont également antiferro- 
magnétiques. Les groupes d’espace de quelques structures magnétiques sont donnés. On en déduit 
certains renseignements sur les différents types de domaines, sur l’orientation exacte des spins, 
et sur l’existence du piézomagnétisme. ; 

Cette étude suggère certains traitements capables de favoriser un type de domaine aux dépens 
des autres. Ce serait le cas, par exemple des cristaux de MnF, ou NiAs franchissant leur point 
de Néel par température décroissante en présence d’un champ magnétique et d’une pression conve- 
nablement orientés. 

Certaines structures donnent naissance à un « ferromagnétisme faible ». Il en est ainsi, par 
exemple, lorsque deux spins antiparallèles sont localisés en des sites non équivalents de sorte que 
les moments magnétiques ne se contrebalancent plus. Il en est encore ainsi lorsque la symétrie 
cristalline n’impose plus l’antiparallélisme rigoureux des divers spins : ce serait le cas de Fe,0,, 
MnCO, et CoCO, (Dzyaloshinsky [18]). 


Abstract. — Tavger and Zaiïtsev [6] have shown that the operations of symmetry which return 
the crystalline medium upon itself are more numerous for an antiferromagnetic crystal than for a 
paramagnetic one. One must indeed consider the ‘‘ anti-operators ? which are the product of a 
geometrical spatial operation by a reversal of time. The number of crystalline classes (groups of 
orientation symmetry) is no longer 32 but 122 (Shubnikov [7]}} Among these 122 classes, 31 are 
ferromagnetic, 59 are necessarily antiferromagnetic (38 of these are piezomagnetic) and 32 are 
either diamagnetic, or paramagnetic, or antiferromagnetic (In this case, there exist two antiparallel 
spin sub-lattices which are related by an elementary antitranslation). Some tables describe the 
properties of the 122 magnetic classes. | 

The number of space-groups is no longer 230 (Schônflies, Fedorov), but 1 651 (Zamorzaev [8], 
Belov, Neronova, Smirnova [9], [10]) 230 are dia or paramagnetic, 272 are ferromagnetic, 517 are 
composed of two antiparallel spin lattices which are related by an elementary antitranslation, 
they are then antiferromagnetic ; the 632 remaining groups are also antiferromagnetic. 

This work suggests some treatments which are able to favour one type of domain at theexpense 
of the others. It should be the case, for instance, of crystalline MnF, or NiAs cooled down through 
their ne point in the simultaneous presence of a magnetic field and a mechanical strain properly 
oriented. 

Some structures induce a ‘‘ feeble ferromagnetism That happens, for instance, when two 
antiparallel spins are located in non equivalent sites so that the magnetic moments does not counter- 
balance any more. That happens also when the crystalline symmetry does not assign a rigorous 
RE Ho to the different spins : it should be the case of Fe,0;,, MnCO;, CoCO, (Dzyalo- 
shinsky [18]). 


LE 


1. Introduction. — L'étude des structures tiques ordonnées s'effectue généralement de la 


magnétiques n’est qu’à ses débuts. Seuls quelques 
cas particuliers simples sont résolus. Cependant, si 
les spins s’orientent à température suffisamment 
basse pour former un édifice triplement périodique, 
on doit s'attendre à découvrir des structures de 
plus en plus complexes chez tous les cristaux conte- 
nant des atomes paramagnétiques. 

La recherche théorique des structures magné- 


façon suivante : 
On suppose connue la structure paramagnétique 
du cristal au-dessus de son point de Curie ou de 


Néel. On peut supposer ensuite qu’il existe une 


direction grossière de ferro ou d’antiferromagné- 
tisme et on cherche empiriquement quelles sont les 
distributions possibles de spins sensiblement paral- 
lèles et antiparallèles. Ces distributions sont 


l 
2 
É 
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_ordonnées par complexité croissante. On peut alors 


chercher lesquelles de ces structures paraissent 
plus stables, moyennant une interaction donnée 
entre spins proches voisins. On cherche ensuite à 
confirmer l’une des structures proposées par les 
méthodes classiques, diffraction de neutrons, dif- 
fraction de rayons X (pour atteindre la défor- 


_ mation spontanée), mesure de susceptibilité, réso- 


nance ferro ou antiferromagnétique, énergie magné- 
tocristalline… 

La recherche des structures possibles a généra- 
lement été effectuée empiriquement (Néel [1], 


Anderson [2], Van Vleck[3]). Wollan et Koehler [4] 


ont tenté une étude systématique en composant 
un réseau cubique simple à l’aide de 7 types de 
cubes élémentaires dont les 8 sommets sont 
occupés par 4 spins (+) et 4 spins (—). Récem- 


_ ment enfin, Gersch et Koehler [5] ont développé 
_ un schéma et une nomenclature simple pour les 


réseaux cubiques et hexagonaux dans le cas où 
chaque spin est environné d’une configuration 
magnétique équivalente. 

Cependant une étude complète des structures 
magnétiques est possible si, en plus des propriétés 
de symétries d’espace on tient compte d’une nou- 
velle opération : le renversement du temps (cette 
opération sera notée 1”). 

Considérons le groupe à deux éléments (12 45et 
effectuons le produit direct des groupes cristallo- 
graphiques de symétrie d'orientation (classes) par 
le groupe (1, 1”). Tavger et Zaitsev [6] ont remarqué 
que l’on obtient ainsi 122 classes (1) qui rendent 
compte des propriétés magnétiques macro- 
scopiques. Ces 122 classes sont constituées : a) par 
les 32 classes usuelles ferromagnétiques, b) par les 
32 classes usuelles astreintes de plus à être symé- 
triques dans le temps et, c) par les 58 groupes 
mixtes de Choubnikov [7]. 

Notre intention est de tirer parti des rensei- 
gnements apportés par la connaissance des 
122 classes magnétiques. Ensuite nous chercherons 
à utiliser les groupes magnétiques d’espace obtenus 
en effectuant le produit direct des 230 groupes 
d’espace de Schôünflies-Fedorov par le groupe (1, 1”). 


: Les groupes magnétiques d’espace sont au nombre 


de 1 651, tabulés par Zamorzaev [8] puis par Belov, 
Neronova et Smirnova [9], [10] sous le nom de 
«groupes de Choubnikov ». 


IT. Les classes et les groupes de Choubnikov. — 
Le travail de l’école russe (2?) (Choubnikov [7], 
Zamorzaev [8], Belov, Neronova et Smirnova [9], 


_ [10]) est la généralisation à 3 dimensions des études 


(:) Tavger et Zaitsev ne donnent que 90 classes magné- 
tiques mais ils oublient les classes symétriques dans le 
temps (grises) qui peuvent être antiferro, para ou dia- 
magnétiques. 

(2?) Ces travaux n'étant malheureusement pas traduits, 
on pourra lire les articles de MacKay [19] ou Le Corre [20] 
sur cette question. 
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de Speiser [11] (une seule dimension) et d’Alexander 
et Herrmann [12], puis Weber [13] (deux dimen- 
sions). Ces auteurs généralisent la notion de 
groupe cristallographique en affectant une couleur, 
soit blanche, soit noire aux atomes du cristal. 

En plus des opérations de symétrie habituelle 
qui transportent un atome coloré sur un atome de 
même couleur il convient de considérer des « anti- 
opérateurs » qui transforment un atome blanc en 
noir et réciproquement. On peut donner plusieurs 
significations à l'opération du changement de 
couleur : changement de signe d’une fonction, 
transformation d’une grandeur imaginaire en sa 
complexe conjuguée, traversée d’une surface ou 
retournement en doigt de gant dans l’espace à trois 
dimensions... Nous l’assimilerons au renversement, 
du temps. Les « anti-opérateurs «seront notés avec 
le même symbolisme que l'opérateur ordinaire 
correspondant, à Ceci près que nous placerons un 
accent en haut et à droite du symbole. 

Aiïnsi l’opérateur 1’ a pour seul effet de changer 
la couleur de chaque atome du cristal. L’opé- 
rateur 3’ (dénommé antiaxe ternaire inverse) con- 
siste à effectuer dans n'importe quel ordre les 
opérations suivantes : rotation ternaire, inversion, 
changement de couleur : 


DS Se 4x 47: 


On constate aisément que 1’ commute avec tous 
les opérateurs de symétrie. 

Les opérateurs de base sont donnés dans le 
tableau I. 


TABLEAU I 


OPÉRATEURS FONDAMENTAUX DÉCRIVANT LES DIVERSES 
SYMÉTRIES CRISTALLOGRAPHIQUES DE! L'ESPACE TEMPS 


“ Directs 12 RSS SR 6 
Éléments de symétrie es pr e. hQRE 
Inverses AMEL 
Directs 7000 ER ER CE 
Antiéléments de symétrie >: LOL 
Inverses DONS EE MECE 
2= m 3—3+1 6 — 3/m 
Nota : = is 2 ee / 
D TE Je le 6—= 3/m’ 


Une convention intéressante consiste à consi- 
dérer une couleur intermédiaire entre le blanc et le 
noir : le gris. Un atome gris ne change pas de 
couleur dans une antiopération de symétrie. 

Considérons donc que l'opération 1’ correspond 
à un renversement du temps. Dans un cristal dia 
ou paramagnétique tous les atomes sont insen- 
sibles en moyenne au renversement du temps, donc 
gris. Au contraire un atome libre dont le moment 

magnétique a une valeur moyenne non nulle 
possède la symétrie co/m 2'/m' qui comporte un 
centre de symétrie. Le moment magnétique moyen 
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est dirigé suivant l’axe d’isotropie (co). Un mirage 
normal à l’axe d’isotropie laisse l’atome inchangé ; 
mais «in mirage parallèle à cet axe renverse le 
moment magnétique ; pour redresser le moment 
magnétique il suffit de changer le signe du temps. 
-L’atome est de même inchangé dans un antiaxe 
binaire normal à l’axe d’isotropie. Si l’atome est 
engagé dans un cristal, sa symétrie propre est 
réduite par celle de son entourage mais reste sous- 
groupe de com 2'/m'. : 

Ces sous-groupes font partie des groupes de 
symétrie d’orientation bicolores ou « classes de 
Choubnikov ». 

Les classes de Choubnikov sont au nombre 
de 122 reproduites dans le tableau II, confor- 
mément au symbolisme de Hermann-Mauguin. Ce 
symbolisme est particulièrement intéressant pour 
représenter les groupes de symétrie, il suffit d’effec- 
tuer une légère réduction pour retrouver le symbo- 
lisme des tables internationales de cristallographie. 


TABLEAU II 


LEs 122 CLASSES MAGNÉTIQUES 
(GROUPES DE SYMÉTRIE D'ORIENTATION) 


D : diamagnétique ou paramagnétique 
F : ferromagnétique 
AF : antiferromagnétique 
P : piézomagnétique. 
L’ordre du groupe est donné dans la dernière colonne 


c 42 D ou AF 2 
Lt ia P il 
At D ou AF a 

CG 1 F P ÿ 
4* AF 2 

( DNA D ou AF A 

Co 2 F 1 2 
(2 AN 

| m 41’ D ou AF k 

CG m F 12 2 
( m' +P) 0 2 

MER NUtA D ou AF 8 

2/m F 12 L 

Con 4427/m AF k 
2/m'’ AF [A 

27}! F 1 4 

ATEN ES UE D ou AF 8 

DS DAT AF IÈ k 
DDR, F P [A 
MAMAN ETS D ou AF 8 

C HTNTR TR AF P 4 
SN PTE ANT RE Je R A 
MANN NT F P 4 
2/m2/m2/m 1’ D où AF 16 
2/m2/m2/m AF F 8 

Da ( 2/m’2//m2'/m AF 8 
2'/m'2’/m'2/m F P 8 

| :2/ra°%2/ma2 nt AF 8 

RO L D ou AF 8 

Ca A F P k 
4° AF 4 


S4 


Can 


Dy 


Cs 


Ci 


D; 


Ce 


Csh 


Con 


Ds 


HZ ET 


ÆUELELEUR 
Lo] 


2 
Lk]m2/m2/m 1’ 
4/m2/m2/m 
4/m’2//m2’/m 
4'/]m2’/m'2/m 
4‘/m’2/m'2/m 
4/m2//m’2//m’ 
k/m'2/m'2/m'’ 


SAS 06 


[SE] 
je 
> 


na 


Co] ol ol 


R9 

>. 
É; 

* 


» 


4’ 


ialel Do co] cal ulcol]co|l 2 co Go -u oo co 
Lo] 
ES 
3 


67 

6) 47 

6/m 

6’/m 

6/m° 

6’/m'’ 

GAP S MITA 
(er EP 
GED 
G+"2"0r27 
GRAS 
6 mem 
6’ m’ m 

6: *°m°’'"m" 


AF 
D ou AF 


i D ou AF 


HU DUT 


golize Pet v se, 


UT 


Re 


_ 


— 


= 
Go O0 OÙ 00 Œ > O0 O0 © A Oo & 
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6 Mon de D ou AF 2% 
6 ms 2 AF P 42 

Dan 46 m’ 2 > AF Pre 
Genre AF 1 12 
Game2" F 15 12 
6/m 2/m 2/m 1° DouAF 48 
6/m 2/m 2/m P 24 
6/m’ 2’/m 2//m AF 24 

Dan { 6//m 2/m’ 2’/m AF 24% 
6//m’ 2//m’ 2/m AF Pr 
6/m 2‘/m’ 2’/m’ F P 24 
6/m’ 2/m° 2/m’ AF 24 

T PC 1 D ou AF 24 
2 3 


TRICLINIQUE 


: MONOCLINIQUE 


2 nes D ou AF 48 

qe | 2/m 3 AF PE 
Ans AF 24 
SAR D ou AF 48 

O &, 3 2 AF P 24 
HS C2 AF P 24% 

RSS ne 40 D ou AF 48 

Fd 4 3-1m AF P 24 
PARENT AF P 24 

&jm 3 92/m 1’ D ou AF 96 

&jm 3 2/m AF P 48 

Oh l4/]m’3° 2m AF 48 
jé 3 2//m’ AF P 48 

48 


Fic. 1 a. — Les 36 réseaux de Bravais bicolores (d’après Belov, Neronova et Smirnova [9], [10]). 
Suite de la figure pages suivantes. 


Le nombre des groupes d’espace s'élève à 1 651 
parmi lesquels il faut compter les 230 groupes ordi- 
_naires de Fedorov puis 230 groupes possédant en 
plus des précédents la symétrie dans le temps et 
auxquels appartiennent donc les cristaux dia et 
paramagnétiques, enfin 1 191 groupes mixtes. 


Parmi les 1 191 groupes mixtes on en compte 517 
qui possèdent, en plus des opérations de symétrie 
des 230 groupes de Fedorov, un certain nombre 
d’antitranslations élémentaires avec les antiopé- 
rateurs qui en résultent, Une antitranslation trans- 
forme un spin en son antiparallèle ; ces groupes 
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sont donc antiferromagnétiques, lä classe corres- 
pondante est grise. 

Ces antitranslations élémentaires sont déter- 
minées par 22 nouveaux réseaux de Bravais décrits 
figure 1 et nôtés avec les indices suivants : 


a, b, c : arête anticentrée, 
À, C : face anticentrée, 
I : corps anticentré, 


triclinique P, arête anticentrée 
s { orthorhombique et cubique : 
F, : arêtes et corps anticentrés. 


N° 10 


La maille élémentaire magnétique est alors un 
multiple pair de la maille paramagnétique (aux 
déformations spontanées près). 


Les 674 groupes restants sont obtenus en rem- | 


plaçant, dans un groupe de Fedorov, au moins un 
opérateur de symétrie (autre qu’une translation) 
par son antiopérateur correspondant. Ces groupes, 
de même que les 230 groupes qui ne comportent 
aucun antiopérateur, sont ferro ou antiferro- 
magnétiques: Pour ces 904 groupes la maille 
magnétique est un multiple pair ou impair de la 


Fi1G. 1 b. — Les 36 réseaux de Bravais bicolores (d’après Beloy, Neronoya et Smirnova [9], [10)]. 
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maille paramagnétique (aux déformations spon- 
tanées près). Les propriétés des classes qui sont 
données plus loin indiquent que 272 groupes sont 
ferromagnétiques et 632 antiferromagnétiques. 


On trouvera un tableau des 1 651 groupes de 
Choubnikov dans les références [9] et [10]. Nos 
notations seront celles de Belov, Neronova et 
Smirnova [10]. 


QUADRATIQUE 


HEXAGONAL ET 


RHOMBOEDRIQUE 


CUBIQUE 


Fi. 4 c. — Les 36 réseaux de Bravais bicolores (d’après Belov, Neronova et Smirnova [9], [10]). 


III. Propriétés des classes magnétiques. — Les 
propriétés macroscopiques d’un monocristal sont 
déterminées par sa classe. 

L'opération «renversement du temps » change le 
‘signe des grandeurs antisymétriques dans le temps 
comme le moment magnétique, le pouvoir pyro- 
magnétique, le pouvoir rotatoire magnétique spon- 


tané, la conductivité rotationnelle électrique ou 
thermique, les tenseurs piézomagnétiques et 
magnéto-optiques… 


Î. FERROMAGNÉTISME, PYROMAGNÉTISME, POU- 
VOIR ROTATOIRE MAGNÉTIQUE SPONTANÉ, CONDUC- 
TIVITÉ ROTATIONNELLE. — La symétrie intrinsèque 
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du ferromagnétisme et du pyromagnétisme est 
celle d’un vecteur axial antisymétrique dans le 
temps (co/m 2'/m'). 

Le pouvoir rotatoire magnétique est déterminé 
par les parties imaginaires des susceptibilités diélec- 
triques et magnétiques. Ces susceptibilités étant 
hermitiennes, leurs parties imaginaires ont la 
symétrie d’un tenseur de rang deux antisymétrique 
dans l’espace et dans le temps, c’est-à-dire celle du 
ferromagnétisme et du pyromagnétisme. C’est 
encore la symétrie de la conductivité rotationnelle 
électrique ou thermique (effets Hall ou Rigghi- 
Leduc spontanés). 

Pour qu'un cristal possède ces propriétés, il faut 
donc que sa symétrie propre soit sous-groupe 
de (co/m 2’/m’). On compte 31 telles classes 
(tableau IIT). 


TABLEAU III 


CLASSES FERROMAGNÉTIQUES, PYROMAGNÉTIQUES, 
DOUÉES DE POUVOIR ROTATOIRE MAGNÉTIQUE SPONTANÉ 
OU DE CONDUCTIVITÉ ÉLECTRIQUE 
ET THERMIQUE ROTATIONNELLE 


Triclinique 1 

Monoclinique Dr NT IN F2 002 2/0 

Orthorhombique 2 2'2’,m’m’2,m m’2’,2/m 2’/m’ 2'/m’ 

Quadratique LEE RO in 0 2 Enr An nr, 
&jm 2//m° 2m’. 

Rhomboédrique 


Hexagonal 
6/m 2'/m' 2'/m’ 


La direction de l’aimantation moyenne est sou- 
vent déterminée : parallèle aux axes 2, 3, 4 ou 6 ou 
normale au miroir. Les exceptions sont les classe 1, 
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1,2’, m’, 2//m', cependant pour les trois dernières, 
laimantation moyenne est parallèle à l’antimiroir 
ou normale à l’antiaxe binaire (tableau IV). 


TABLEAU IV 


COMPOSANTES DES TENSEURS DU SECOND ORDRE ANTISYMÉ- 
TRIQUES DANS L'ESPACE ET DANS LE TEMPS ET COMPO- 
SANTES DES PSEUDO-VECTEURS CORRESPONDANTS 


Ces compôsantes décrivent l’aimantation spontanée, le 
pouvoir pyromagnétique, la partie imaginaire des suscep- 
tibilités diélectriques ou magnétiques (et, par conséquent, 
le pouvoir rotatoire magnétique spontané), la conductivité 
rotationnelle, électrique et thermique. L’axe non primé de 
plus grand ordre est pris suivant Ox.. ; 


sé 5 12 13 1 
1, 1 PE CPR 2 
—13 —923 :. 3 
hu 13 fs | 
2) 2 


AN 


ue — 3 


: 12 
— 12 : 


2. PIÉZOMAGNÉTISME. — Le piézomagnétisme 
est déterminé par un pseudo-tenseur du troisième 
ordre IT, 8, tel que 


M, = DIE pe Ts 
œ Ts a, Py !6y 


Toutes les autres 
classes ferromagnétiques 


Conformément à la théorie de Laval, nous consi- 
dérerons le tenseur comme dissymétrique en B et y. 

Il est aisé de construire par inspection directe 
le tableau des constantes piézomagnétiques 
(tableau V). 


TABLEAU V 


COMPOSANTES DES TENSEURS PIÉZOMAGNÉTIQUES 


Le deuxième indice j est affecté au tenseur d’effort ou de déformation X,g'; il court de 1 à 9 conformément au tableau 


suivant : 


MES EP AE MT 


DO ZA SD) 


ai — 1122 33 23 31 12 32 13 21 


à 18. 


Si le tenseur des efforts est symétrique, le nombre des composantes utiles s’abaisse, pour un cristal triclinique, de 27 


Pour les classes non cubiques, l’axe de plus grand ordre est pris suivant Ox;, s’il existe de plus des axes 2 ou 2 (=m) 
perpendiculaires, l’un d’eux est pris suivant Ox;, s’il n’y a pas d’axe 2 ou 2 perpendiculaires mais des axes 2’ ou 2’ (— m’) 


perpendiculaires, l’un d’eux est encore pris suivant Ox:. 


Pour les classes cubiques Ox;, Ox, et Ox, sont portés par les axes quaternaires du cube superposable. 
Le symbole F rappelle quelles classes sont ferromagnétiques, é 


AL OA AT A0 16 217 RAS AREO 
210220230024 0029 026 27028020 
91 32 1139 34% 9352867137, 238039 


2/m 2/m 2/m 


1% 17 
25 28 
36 39 |; 


_ No10 
| 
l HA" mme, .2/m'm, 2/m2/m2/m 
14-41 
36 
: &, 4, am 
18 15 2504 
24956 AS, 14 
31—31 36 36 
F 4, 4; 4]m:-6,"6, 6/m 
14 AB Vois 96 
24 95 1514 
31 31 33 36 —36 
4 2m’, 422’, &m'm’ &/m 2m’ 2m’ 
FOR E Nm 06227; 6m/m", 6jix. 2’)m’ 2’/m 


On constate que 69 classes sont piézo- 
magnétiques mais, parmi celles-ci, si le tenseur des 
efforts appliqués est symétrique, il faut supprimer 


les trois classes cubiques 4 3 2, 4 3 m et 4/m 3 2/m. 

On remarque que toute classe non grise qui ne 
possède pas un anticentre de symétrie est piézo- 
magnétique. 


3. CHANGEMENT DE PHASE. — Si On abaisse la 
température d’un cristal au-dessous de son point 
de Curie ou de Néel, il se divise en domaines. Seul le 
domaine possède la triple périodicité. Le domaine 
appartient à une certaine classe. Le seul rensei- 
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F car 
11—11 1% 415—22—25 24-922 
—22 22 2h 2511 15—14—11 
19 1MSS 36 —36 
F ASTM NS 27/m 
15—22 24—22 
—22 22 24 15 
31: 31 33 
LÉO EM RMAT ER" re22% 
eme net m9 7m 
14 25 
25 14 
36 36 | 
| 
& 2m, 422, 4mm, 4£/m 2/m 2/m 
6. m 2, 622, 6mm, 6/m 2/m 2/m 
14 —25 
25 —14 
36 —36 
| 
6’ 2’ m, 6’ 226 2/24%6/m'm, 
Cm men; n 
11—11 
—11 —11 
23, :2/m. 3 
14 47 
44 47 
14 17 


gnement qui nous est apporté par le principe de 
symétrie de Curie est le-suivant : la classe du 
domaine est sous-groupe de la classe paramagné- 
tique du cristal. Nous ne pouvons énoncer aucune 
règle sur d'éventuelles transitions entre états non 
paramagnétiques ; en particulier, il ne faut pas 
penser que la symétrie des domaines de chaque 
phase décroît quand ces phases sont obtenues par 
abaissement de température. Les classes donnent 
des renseignements sur les assemblages de domaines 
qui ne sont autres que des macles par pseudo- 
symétrie. 

Le nombre des types distincts de domaines est 
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égal au rapport des ordres de la classe du cristal 
paramagnétique et de la classe du domaine. 

Ainsi, si la transition de l’hématite Fe,O,x 
à 950 0K se fait avec passage de la classe 3 2/m {1 
(ordre 24) à la classe 2/m (ordre 4), il faut s'attendre 
à trouver six types distincts de domaines. La tran- 
sition à 250 0K qui rehausse la symétrie à la 
classe 3 2/m (ordre 12) réduit le nombre des types 
de domaines à deux. Les ordres des 122 classes 
sont donnés dans le tableau IT. 


4. CROISSANCE PRÉFÉRENTIELLE D'UN TŸPE 
DE DOMAINE DANS UN MONOCRISTAL, — Consi- 


PR 4/m nec 


arrangement 


Ca : cv a V2 : 
cubique centré de 
première espèce cv2Va: 
cvaV? : arrangement 
c.f.c. de première 


espèce. 


Ie 4,/a cd 


c.c. de deuxième espèce 


arrangement cvaV2: arrangement 
c.fic. de deuxième 
arrangement espece. 


c.f.c. amélioré de pre- 
miere espece. 


F 4,/m bc PF 4/mecc 


No 10 
dérons un cristal subissant le changement de phase 
&]m 2m 92/m1-—>4//m 2/m 2/m-’. 


Il apparaît deux types de domaines. On passe 
d’un type à l’autre par rotation de 90° autour de 
l’ancien axe quaternaire. Les composantes du 
tenseur piézomagnétique sont respectivement : 


le 4,/a cd 


C 


F1G. 2, — Symétrie des structures proposées par Van Vleck, Anderson, Wollan et Koehler. 
Les atomes portent des spins dirigés suivant l’axe quaternaires. 
Selon le sens du spin, l'atome est coloré en blanc ou en noir, 


No 40 


et 


Si on applique un cisaillement T4, un type de 
domaine prend l’aimantation 


MS OU Te 


et l’autre type prend une aimantation opposée de 
sorte que l’ensemble du cristal n’est pas aimanté. 
Mais en appliquant simultanément un champ 
magnétique le long de l’axe quaternaire, on abaisse 
l'énergie potentielle d’un type de domaine aux 
dépens de l’autre. Si le champ magnétique et le 
cisaillement sont appliqués au cours du chan- 
gement de phase, on favorise le développement d’un 
type de domaine et, après suppression du champ 
magnétique, le cristal reste piézomagnétique. 

Il en serait de même pour un cisaillement 7, 
(ou T';) et un champ parallèle à un axe binaire. 

Le tableau V, p. 756, permet de déterminer le 
processus d'orientation adapté à toute transition. 


IV. Symétrie de quelques structures magné- 
tiques simples. — 1. STRUCTURES POSSÉDANT UNE 
DIRECTION D'ANTIFERROMAGNÉTISME. — Néel [1|, 

. Anderson [2], Van Vleck[31, Wollan [4], Gersch [5], 
Koehler [4], [5].…., ont étudié les propriétés de 
quelques structures magnétiques simples que 


Fic. 3. — Exemples de structures à ferromagnétisme faible. 
Les spins sont parallèles à l’axe quaternaire. Selon le 
sens du spin, l’atome est coloré en blanc ou en noir, 
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pourraient présenter un cristal cubique ou hexa- 
gonal. 

Nous donnons dans la figure 2 la symétrie des 
structures quadratiques proposées par ces auteurs, 
dans le cas où la direction unique d’antifer- 
romagnétisme est portée par l’axe quaternaire. 
Il serait aisé de poursuivre cette étude dans le 
cas où la direction d’antiferromagnétisme est 
différente. On remarque que la plupart des struc- 
tures proposées sont constituées de deux sous- 
réseaux anti-parallèles et appartiennent à des 
classes grises. 

Toutes les structures de la figure 2 possèdent un 
anticentre de symétrie, elles ne sont donc pas 
piézomagnétiques. Notons qu’il en serait encore 
ainsi si la direction d’antiferromagnétisme était 
quelconque. 


2. STRUCTURES FAIBLEMENT MAGNÉTIQUES. — 
On note sur la figure 3 (où les spins sont toujours 
supposés parallèles à l’axe quaternaire), que la 
structure quadratique de gauche appartient à la 
classe ferromagnétique 4/m 2'/m' 2’/m’. Cepen- 
dant il y a autant de spins qui pointent dans un 
sens que dans L'autre et les spins sont bien paral- 
lèles en moyenne à l’axe quaternaire (les atomes 
porteurs sont dans des plans de symétrie) mais la 
symétrie cristalline n’impose pas la nullité rigou- 
reuse du moment magnétique global. Il faut 
s'attendre à un ferromagnétisme faible causé par 
les différents modes d’interaction entre voisins. 

La deuxième structure portée sur la même figure 
dérive du type D.E ; elle appartient à la même 
classe cristalline. Ici, la symétrie cristalline n’im- 
pose plus l’antiparallélisme rigoureux des moments 
magnétiques moyens et entraîne l'existence d’un 
ferromagnétisme faible. 

Nous verrons ultérieurement que, dans Fe,0;, 
MnCO; et CoCO3, la symétrie cristalline n’impose 
plus l’antiparallélisme rigoureux des spins, il n’y a 
aucune raison pour que le moment magnétique 
global soit nul. 


3. STRUCTURES A PLUSIEURS DIRECTIONS D’ANTI- 
FERROMAGNÉTISME. — Il est facile d'imaginer des 
structures quadratiques possédant deux directions 
d’antiferromagnétisme orthogonales. Nous en don- 
nons quelques-unes sur la figure 4. 

Sur cette même figure nous donnons également 
une structure cubique où les quatre directions 
d’antiferromagnétisme sont portées par les grandes 
diagonales du cube. 

Ces cinq structures ont un anticentre et ne sont 
pas piézomagnétiques. 


V. Exemples de structures magnétiques. —— Les 
structures que nous allons étudier, à part MnBr, 
sont décrites dans l’article de mise au point de 
Nagamiya, Yoshida et Kubo [14]. 
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1. MxO, FEO, CoO, N1O, MnS, MNSE. — Ces 
cristaux ont à l’état paramagnétique la structure 
cubique à face centrée du chlorure de sodium 


(fig. 5). F 4]m 3 2/m 1’(O$), tous les ions para- 


Pe 4h bm 
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magnétiques d’un même plan [111] ont leurs spins 
parallèles. Les spins des divers plans [111] sont 
alternés. La structure magnétique dérive donc de 
l’un des nouveaux réseaux de Bravais. 


le 4mbm Pcé/n mm Icé/mmm PB, n3n 


Fic. 4. — Exemples de structures à plusieurs directions d’antiferromagnétisme. 


Les structures magnétiques diffèrent selon la 
direction d’antiferromagnétisme : 

a) Direction [0 0 1] de cube ou, plus généra- 
lement, toute direction perpendiculaire à un axe 
binaire [110], [101] ou [0111]. 

La structure magnétique est monoclinique : 
Ce 2/c. 

(L’axe binaire est l’ancien axe binaire du 
cube [110] et l’axe ce peut être choisi suivant 
l’ancien axe [001 /]). 


ONa @CI 


F1G. 5. — Structure du CINa: C’est également la structure 
de MnO, FeO, CoO, NiO, MnS, MnSe à l’état para- 
magnétique. 


La classe est  2/m 1” (ordre 8). 

La classe paramagnétique était  4/m 3 2/m 1’ 
(ordre 96). 

Il faut done attendre 12 types distincts de 
domaines. 

Il n’y a aucune raison pour que la direction 
d’antiferromagnétisme soit dirigée rigoureusement 
suivant un axe quaternaire du cube. 


b) Direction [111] du cube. 

L’axe ternaire est conservé, le groupe magnétique 
est rhomboédrique R;32/c, appartenant à la 
classe 32/m 1” (ordre 24). On doit attendre 
4 types de domaines selon que la direction d’anti- 
ferromagnétisme est suivant l’un des 4 axes ter- 
naires du cube. 


c) Suivant l’un des axe binaires du cube paral- 


lèles à (111), par exemple [011]. 

L’axe binaire est conservé, le réseau est à 
nouveau monoclinique C;. 

Le groupes est C.2/m et la classe 2/m 1’. 

Il faut s’attendre à trouver 12 types de domaines. 

d) Direction générale : le groupe est P,1, la 
classe 1, 1” (ordre 4). Il y aurait 24 types de 
domaines. 


2. MNF,, CoF,, FEF,, N1F,, MO. — Ces 
cristaux ont à l’état paramagnétique la structure 


O0eT: 


Fi6. 6. — Structure du rutile TiO.. 
C'est également la structure de MnF,, CoF,, FeF,, 
NiF,, MnO, à l’état paramagnétique. 


du rutile (fig. 6) DK — P 4,/m n m 1’ de Ja 
classe 4/m 2/m 2/m 1’ (ordre 32). 

— La structure antiferromagnétique la plus 
simple est obtenue avec une maille quadratique 
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voisine de la maille paramagnétique et contenant 
encore deux molécules. Les spins sont néces- 
sairement antiparallèles et dirigés suivant l’axe 
quaternaire €. Le groupe de symétrie devient 
P 4;/mnm', appartenant à la classe piézo- 
magnétique 4'/m 2/m 2’/m’ (ordre 16). Ce serait 
effectivement le structure de MnF,, FeF, et CoF; 
Erickson [15]). 


F1G. 7. — Structure de. NiAs. 
C’est également la structure de CrS, FeS, MnTe, FeTe, 
CrSb, MnAs, MnBi à l’état paramagnétique. 


— Dans NiF,, d’après les mêmes auteurs, la 
direction des spins ferait un angle de 100 avec 
l’axe c. On note que s’il existe une seule direction 
d’antiferromagnétisme, la structure peut présenter 
une maille orthorhombique contenant encore deux 
atomes par maille. : 

— La structure proposée pour MnO, par Erick- 
son conserve la symétrie quadratique mais la 
maille magnétique contient quatre mailles para- 
magnétiques : les atomes situés aux sommets et au 
centre de la maille quadratique forment deux édi- 
fices antiferromagnétiques (réseau de Bravais P.:) 


F1G. 8. — Structure de CuCl,, 2H,0. 
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reliés entre eux par un axe quaternaire hélicoïdal : 
c'est la structure Pc 4,/n de la classe 4/m 1’ 
(ordre 16). 

Les spins sont astreints à être normaux à l’axe 
quaternaire mais ne sont pas nécessairement dirigés 
vers les oxygènes voisins. 


3. NrAs, CRS, FES, MNTE, FETE, CRSB, MAS, 
MxXBr. — Ces cristaux ont à l’état paramagnétique 
la même structure (fig. 7). Le groupe est : 

Dé — C 64/m m c 1’ et la classe : 6/m 2/m 2/m 1’ 
(ordre 48). 

La maille paramagnétique contient deux molé- 
cules. 

Expérimentalement, la maille magnétique est 
mal connue et paraît sensible à la présence d’impu- 
retés. Des études aux rayons X ont révélé dans 
certains cas des déformations orthorhombiques. 

Si la maille magnétique restait hexagonale, voi- 


Fic. 9, — Structure de Fe,O3x et de Cr,:03. 


sine de la maille paramagnétique, les spins seraient 
nécessairement dirigés suivant l’axe senaire. Le 
couplage antiferromagnétique impliquerait la struc- 
ture P6',;/m'm’c appartenant à la classe piézo- 
magnétique 6'//m’2’/m'2/m. 

D’après Snow [16], CrSb aurait effectivement 
cette structure. 
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4. CuCLs, 2H,0. — La structure de CuCl, a été 
très étudiée (fig. 8). 

La structure paramagnétique est 

Dan —Pmnal 
(classe orthorhombique  2/m 2/m 2/m 1”). 

D’après Poulis et Hardeman [17] la maille 
magnétique est doublée selon l’axe c. Le groupe 
serait donc: P.ben (spins suivant a). 

En présence d’un champ magnétique suffisant, 
dirigé parallèlement à l’axe a, les spins s’orientent 
suivant l’axe b, la structure devenant  P, n mn. 
Les deux structures antiferromagnétiques appar- 
tiennent aux deux classes grises 2/m 2/m 2/m 1’. 


5. FE,O3x, CR203. — La maille est un rhom- 
boëdre contenant deux molécules (fig. 9). Le groupe 
cristallographique, à haute température, est 
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R 3 2/c 1” de la classe 32/m1' (ordre 24). 

À l'intérieur de la maille, les quatre atomes para- 
magnétiques sont sur l’axe ternaire du rhomboëèdre. 
L'étude thermodynamique de ces cristaux a été 


donnée récemment par Dzyaloshinsky [18] qui fait 


lui-même appel aux considérations de symétrie. 

Cette étude est possible en supposant : 19 que la 
maille magnétique est voisine de la maille para- 
magnétique ; 2° que les spins sont sensiblement 
parallèles ou antiparallèles. 

On distingue quatre types de couplages magné- 
tiques (fig. 10) ; les deux premiers sont caractérisés 
par un centre de symétrie et les deux derniers par 
un anticentre. 

Le type I est un couplage ferromagnétique, les 
autres sont antiferromagnétiques. 

Dans la légende qui accompagne la figure 10, on 


O le) . © 0: 
# 
[o) e e [e) 
[e] e Oo e 
O [e] e [2 
#5 au I a: 
Spins parallèles à l’axe R 3 2’/c’ R 3 2/c R 3/2/c’ R 3/2//c Rhomboédrique 
ternaire 
Cas, par exemple, des C2/e G2’Jc". C2’/c G2}c? Monoclinique 
spins presque parallèles (ferro-magnétique) (ferro-magnétique) 
à l’axe binaire ‘ 
Cas, par exemple, des G2/c’ G2/c C2/c’ G2//c Monoclinique 
spins presque normaux (ferro-magnétique) (ferro-magnétique) 
à l’axe binaire 
Cas général P1 P1 P<’ P1’ Triclinique 


(ferro-magnétique) (ferro-magnétique) 


F1G. 10. — Les différents types de couplages dans la structure Fe,0O;x. 
Les cercles représentent les atomes magnétiques 1, 2, 3, 4 alignés sur l’axe ternaire. On suppose que la maille magné- 


tique est identique à la maille paramagnétique. 


Les atomes blancs ont des spins sensiblement parallèles, les spins des atomes noirs sont sensiblement antiparallèle à ceux 


des atomes blancs. 


Le tableau indique la symétrie correspondant aux diverses directions moyennes des spins. 


spécifie le groupe de symétrie correspondant aux 
diverses orientations moyennes de spins. 

D’après Brockhouse, l’hématite Fe,0,x a un 
couplage du type II. A basse température, en 
dessous de 250 °K, les spins seraient parallèles à 
l'axe ternaire, conformément au groupe R 3 2/c. 
Mais de 250 0K au point de Néel (950 0K), les 
spins ne sont pas parallèles à l’axe ternaire. On 
constate alors que, quelle que soit leur orientation, 
la symétrie est compatible avec le ferromagnétisme 
et c’est probablement l’explication du ferroma- 
gnétisme faible de l’hématite. 

En effet, la symétrie n’impose plus au moment 


> 
magnétique global m d’être nul. La structure peut 


ÆS 
être monoclinique C 2’/c’ (ni normal à l'axe 


binaire) ou CG 2/c Ex parallèle à l’axe binaire] - 
dans le cas le plus général, elle serait triclinique P 1. 

Au contraire Cr,0, présente un couplage du 
type III, quelle que soit l’orientation des spins, la 
structure possède un anticentre de symétrie et est 
donc antiferromagnétique, dénuée de piézoma- 
gnétisme. 


6. MxCO;, CoCO;, FECO;. — La structure para- 
magnétique est celle de la calcite. Le groupe est le 


même que pour Fe,03: D$a — R 3 2/e 1’. La 
maille contient deux molécules (fig. 11). 


? 
l 
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S'il existe une direction approchée d’antiferro- 
magnétisme et si la maille magnétique est voisine 
de la maille paramagnétique, plusieurs cas sont à 
considérer : 

a) Direction d’antiferromagnétisme parallèle à 
l’axe ternaire. 

Le groupe est R32/c, les spins sont rigou- 
reusement parallèles à l’axe ternaire. 

Chaque domaine est piézomagnétique. Ce Free 
le cas de CO;Fe. 


eCa 
EC 
OO 


Fic. 11. — Structure de la calcite. 
C’est également la structure paramagnétique de 
MnCO;, CoCO;, FeCO; 


b) Direction approchée d’antiferromagnétisme 
parallèle à l’un des axes binaire. La structure 
devient monoclinique et le groupe C 2’/c’. Aucune 
condition de symétrie n’impose aux spins d’être 
rigoureusement antiparallèles : il y a ferromagné- 
tisme faible. 

Le moment magnétique global est normal à 
l’axe binaire. 

c) Direction approchée d’antiferromagnétisme 
normale à l’un des axes binaires. 

Le groupe est monoclinique C 2/c. Il y a ferro- 
magnétisme faible suivant l’axe binaire. 

d) Direction approchée d’antiferromagnétisme 
quelconque, le groupe devient P 1’. 

Il y a encore ferromagnétisme faible dans une 
direction sensiblement normale à la direction appro- 
chée d’antiferromagnétisme. 

Les propriétés de MnCO; et CoCO; appar- 
tiennent à l’un des trois derniers cas. 


7 MxBr. — Certaines structures magnétiques 
ont été proposées pour les halogénures de Mn, Fe, 
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Co, Ni. Récemment, Koehler, Wilkinson, Cable 
et Wollan [21] ont réussi à étudier un monocristal 
de MnBr,, rendu sensiblement monodomaine à 
l’aide d’un champ magnétique, depuis 1,81 °K 
jusqu’à 25 ©K. La structure paramagnétique est 
celle de Cdl,;, du groupe P 32/m 1 1" de la 
classe 3 2/m 1” (ordre 24). La structure para- 
magnétique proposée par ces auteurs est pseudo- 
orthorhombique, chaque domaine ayant une direc- 
tion d’antiferromagnétisme parallèle à l’un des 
axes binaires. En toute rigueur cette structure est 
monoclinique (fig. 12) et appartient au groupe 


Fire. 12. — Structure magnétique de MnBr, (d’après [21]). 
Les atomes de Br ne sont pas indiqués, la direction 
d’antiferromagnétisme est suivant b, parallèle aux axes 
binaires. Les atomes sont colorés en blanc ou en noir 
selon le sens du spin. 


Ce 2/m de la classe 2/m 1 (ordre 8). La maille 
simple contient quatre molécules. Ainsi que ces 
auteurs ont pu l’observer, le nombre des types de 
domaines est limité à 3. Les miroirs imposent à 
tous les spins d’être rigoureusement parallèles à 
l’axe binaire. 


VI. Conclusion. — Les groupes de Choubnikov 
sont particulièrement adaptés à l’étude des cris- 
taux ferro ou antiferromagnétiques. Toutes les 
propriétés physiques du domaine triplement pério- 
dique sont déterminées par les symétries d’espace 
et de temps. Les grandeurs macroscopiques (aiman- 
tation spontanée, piézomagnétisme, conductivité, 
pouvoir rotatoire, énergie magnétocristalline) sont 
imposées par le groupe de symétrie d'orientation 
(ou classe). 

Les grandeurs microscopiques (facteur de Landé, 
couplage spin-spin, interaction spin-réseau, fac- 
teurs de structures aux neutrons.) sont déter- 
minées par le groupe d’espace, lequel renseignesur la 
symétrie du point où réside l'atome para- 
magnétique. La théorie des groupes donne les 
représentations irréductibles des systèmes de spin 
et facilite, par là même, l'établissement des théories 
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atomiques de l’antiferromagnétisme dans les cris- 
taux les plus complexes. Les groupes magnétiques 
constituent un outil mathématique puissant 
capable de lever la part d’arbitraire qui règne dans 
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la détermination des structures magnétiques ordon- 
nées. 
Je remercie À. Herpin et P. Mériel de leurs 
remarques. 
Manuscrit reçu le 11 juillet 1958. 
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DIFFUSION PAR FLUORESCENCE NUCLÉAIRE DU RAYONNEMENT y DE ‘5$e 
PÉRIODES DES NIVEAUX DE 265 ET 402 keV DE ‘As 


Par H. LANGEVIN-JOLIOT et M. LANGEVIN, 


Laboratoire de Physique Nucléaire du Collège de France à Orsay. 


Résumé. — Aucun effet de diffusion résonnante n’a pu être observé dans le cas de la transition : 
directe de 402 keV ; par contre, le rayonnement de fluorescence nucléaire du niveau de 265 keV a 
été mis en évidence. Les périodes des niveaux de 402 keV et 265 keV ont été ainsi fixées respec- 
tivement ait, 1,4 107 sec et +, — (2,3 + 0,3) 10-11 sec, 

L'étude systématique de l’intensité du rayonnement diffusé de 265 keV en fonction de la pression 
de la source gazeuse permet d’estimer à (1,5 + 0,4) 10° sec la période du niveau de 402 keV, 

Les facteurs d’embranchements de capture électronique vers les différents niveaux de "As ont 
été réexaminés. Un faible embranchement supplémentaire (0,04 %) a été observé vers le niveau 


de 575 keV. 


Abstract. — Nuclear resonance fluorescence from the 265 keV excited state of 75A$ leads to a 


value of (2,3 + 0,3) 10—11sec for the half life of this level. 


No resonance fluorescence from the 


402 keV level was observed and the half life of this level is therefore greater than 1,4 10—11 sec. 
The variation of nuclear resonance intensity from the 265 keV level with the pressure of the 
gaseous source gives the value (1,5 + 0,4) 10—° sec for the half life of the 402 keV level. 
Electron capture branching ratios have been reexamined and a weak (0,04 %) new transition 


leads to the 575 keV level of 75As. 


À. Introduction. — L'étude de la diffusion des 
photons par fluorescence nucléaire permet de déter- 
miner les périodes des niveaux excités ; toutefois, 
les sections efficaces correspondant à ce processus 
sont, sauf dans les conditions requises pour la 
résonance, très inférieures à celles de la diffusion 
élastique Rayleigh. 

Sauf dans le cas où la période du niveau consi- 
déré est extrêmement brève (< 10-14 sec) la lar- 
geur l'est petite devant l’énergie cédée au noyau 
de recul lors de l'émission ou de l'absorption 
(quelques eV) : en conséquence la diffusion par 
fluorescence reste inobservable même si elle a lieu 
sur des noyaux de même nature que les noyaux 
émetteurs. Dans certaines conditions, il est cepen- 
dant possible d'utiliser le changement d'énergie du 
rayonnement y par effet Dôppler pour retrouver 
la résonance ou s’en approcher. 

On y est parvenu notamment en donnant à la 
source une vitesse mécanique élevée {1], ou en la 
portant à haute température [2]. Les conditions 
expérimentales nécessaires limitent l’usage de ces 
méthodes au cas de niveaux de faible énergie 
existant dans des noyaux lourds. 

L’élargissement Dôppler dû à la vitesse de recul 
du noyau à la suite d’une émission £ ou y abou- 
tissant au niveau considéré, doit également per- 
mettre, s’il est suffisant, l'observation de la dif- 
fusion résonnante. Malheureusement, dans les 
sources solides, le noyau de recul se ralentit par 
collision en un temps de l’ordre de 10-18 s; ce qui 
explique l’échec des tentatives faites en ce sens par 
différents auteurs [3]. Le premier résultat positif a 
été obtenu par Ilakowac avec une source liquide [4]. 


L'utilisation de sources gazeuses, enfin, dans les- 
quelles les noyaux se ralentissent en des temps 
beaucoup plus longs, a permis à Metzger d'observer 
pour la première fois le rayonnement de fluores- 
cence nucléaire, sans être pratiquement gêné par 
les autres diffusions élastiques [5], [6]. 

Il nous a paru intéressant d'appliquer cette 
méthode à l’étude des niveaux de As. 75As est en 
effet le seul isotope stable de l’arsenic, et 5Se, qui 
se désintègre par capture électronique vers le 
niveau de 402 keV de 75As, peut facilement être 
mis sous forme d’un composé gazeux. L'énergie dis- 
ponible pour la capture vers ce niveau est 
de 470 keV [7] donc supérieure à l’énergie des tran- 
sitions ÿ (voir le schéma fig. 7) ; on pouvait espérer 
que le rayonnement de 402 keV en particulier, 
caractérisé comme Æ, ou M, + E, aurait une 
période suffisamment courte. 

Nous n’avons pu observer aucun effet de dif- 
fusion résonnante dans le cas de la transition 
directe de 402 keV, par contre nous avons pu 
mettre en évidence le rayonnement diffusé sur le 
niveau de 265 keV, et en fixer la période. L’étude 
systématique de l'effet des collisions subies 
par As*, sur l'intensité du rayonnement diffusé 
de 265 keV a permis d’atteindre la période du 
niveau de 402 keV. 

Au moment où s’achevait ce travail, nous avons 
pris connaissance des résultats de F. Metzger sur 
la période du niveau de 265 keV [8]. 


B. Généralités : Relation entre la section effi- 
cace moyenne de diffusion par résonance et la 
période du niveau. — La section efficace de dif- 
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fusion des photons d'énergie Æ, sur un niveau 
d'énergie Æ, est donnée par la formule (1 2 


SN ann ess 
SE) = Br (Er — Ent + GI DE () 
F'et l, sont respectivement la largeur totale et la 
largeur partielle pour l’émission ÿ vers le niveau 
fondamental, et : 


ne 12 À ATX A 
Er = Fo +3 F8 MC RE 
Par ailleurs, l’énergie d’un rayonnement y émis 
par un noyau dont la vitesse a la composante v 
dans cette direction est en première approxi- 
mation : 
E, = Er + Po 


SU Ur) OÙ Ur = 


0: 

o(v) étant alors la section efficace de diffusion pour 
les photons émis par les noyaux de composante de 
vitesse », et P(v) le nombre relatif de ces noyaux, 


On a : 
v max 


P{(v) o(v) du 


RUE da v max 


et comme o(v) est pratiquement nul, sauf dans la 
région très étroite de la résonance on obtient : 


12 ‘18 C 
Omoy — rs 8 Tr Æ. (or) ; (3a) 
RS = Par AC 
Di ter ete P(vr). (3b) 


En réalité, les molécules qui renferment les 
noyaux émetteurs ou diffuseurs ne sont pas au 
repos. V, V1, V, étant respectivement les compo- 
santes dans la direction du rayonnement y, de la 
vitesse propre du noyau émetteur, et des vitesses 
d’agitation thermique on a : 


Mit 
st PR Mi T. 
os = [ff FPS TS 


M; Mo” 2? 
2kT3 a, 
X Ru FT, ° o(v + v, — v,) dv dv, du, (4) 


en admettant que la répartition des vitesses dans 
la source et dans le diffuseur est celle d’un gaz 
parfait. 

Si la largeur propre l' du niveau est petite 


V2kT,]M de 


la source ou À, du diffuseur, la section efficace 
moyenne prend la forme (5) (à condition que P(v) 
ne soit pas nul dans la région de la résonance) : 


devant la largeur Dôppler A; =" 


+ max d E e— TR 
— — 1 — C{Ai?+ A5!) 
FE a: PVR RE + A A 
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— Lorsque P(v) est constant ou varie linéai- 
rement sur l'étendue de la raie Dôüppler, on 
retrouve POUT Gmoy le Même résultat que précé- 
demment (3a). Metzger a récemment préconisé et 
utilisé pour la détermination de la période d’un 


niveau [6], [8], non pas l’étude directe de la section 


efficace de diffusion par résonance, mais celle du 
coefficient d'absorption propre de la raie de réso- 
nance. Dans ce cas, à condition que P(v) varie 
peu sur l'étendue de la raie Dôppler, l'effet ther- 
mique de la source n'intervient pas, non plus que 
la probabilité P(v3): Iatténuation de la raie 
(corrigèe des effets Compton et photoélectrique) 
est en première approximation : à 

6 


F (nd : nombre d’atomes/em?). 


1 
Ain 


EXPRESSION DE LA PROBABILITÉ P(Vr). — La 
probabilité pour un noyau de recul dont la vitesse V 
est distribuée isotopiquement, d’avoir une compo- 
sante de vitesse v dans une direction donnée est 
constante entre — V et + V. P(v) dv = dv/2V. 

a) Le niveau étudié est alimenté par une cap- 
ture électronique directe. La vitesse de recul du 
noyau, constante (aux effets de réarrangement du 
cortège électronique près) est Vx = &/MC ; a étant 
l’énergie disponible pour la capture 
Si 5 Se P(v»). = 1/2V3. 

b) Le niveau est alimenté par une capture élec- 
tronique suivie d’une émission y ; en l’absence de 
corrélation angulaire y — y et V' étant la vitesse 
communiquée au noyau par l’émission du rayon- 
nement y précédant le niveau : 


Vr >vR+ V"  P(vr) = 1/2Vr (7) 
et 
Ur — V' < Vy <UrR + V° 
P{vr) = (1/2V%) (Vx — vr + V')/2V7. 


C. Étude expérimentale de la diffusion réson- 


nante sur les niveaux de *’As. — PRÉPARATION DE 
LA SOURCE. — 7%$e (121 J) a été obtenu en irra- 
diant pendant six semaines du sélénium naturel 
dans un canal à haut flux de la pile P, de Saclay 
(1,8 101? n/cm?/s) ; aucune activité y parasite n’a 
pu être décelée dans un échantillon de ce sélénium 
irradié. 

Le composé gazeux choisi, SeH,, forme par 
capture électronique un radical *AsH, dont l’éner- 
gie de dissociation est supérieure à l’énergie de 
recul. SeH, a été préparé par passage d’un mélange 
de vapeur de sélénium actif et d'hydrogène sur de 
la pierre ponce chauffée à 4000, (Le rendement de 
cette opération est supérieur à 50 %.) Le gaz actif 
traverse ensuite un desséchant (P,0;) puis se con- 
dense dans un piège à air liquide. Le vide étant 
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établi dans l'appareil, on transfère par recon- 
densation une fraction donnée de SeH, formé dans 
une ampoule de verre de 7 cm° qu’il suffit alors de 
sceller sous vide. 8 sources d'activités réparties 
entre 670 pG et 9 uC ont été ainsi préparées, la 
pression de SeH, étant alors de 89 cm de mercure 
dans la première ampoule et de 1,2 cm dans la 
dernière. 


— DisPOSiTIF EXPÉRIMENTAL. — Ce dispositif 
est représenté schématiquement figure 1. Le cristal 


SOURCE 


Error 


de INa — TI de 1” 1/2 de diamètre et de 1” d’épais- 
seur, est protégé du rayonnement y émis direc- 
tement par la source, par un écran de plomb tron- 
conique. Le diffuseur circulaire d’arsenic peut être 
rapidement remplacé par un diffuseur de sélénium 
de masse superficielle équivalente, ce qui permet 
d'éliminer tous les effets autre que la diffusion 
résonnante. Un écran de plomb de 1 mm d’épais- 
seur entourant le cristal absorbe une grande partie 


des photons diffusés par effet Compton, ce qui : 


permet avec les sources fortes d'étudier le pic 
photoélectrique dû aux photons de fluorescence 
sans être gêné par l’empilement des impulsions. 
Le spectre d’impulsions est analysé par un sélec- 
teur d'amplitude à 50 canaux. 

Les diffuseurs sont constitués par de l’arsenic ou 
du sélénium en poudre remplissant un boitier 
cylindrique d'aluminium de 2/10 mm d’épaisseur, 
assemblé à l’araldite ; les diffuseurs minces et épais 
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sont chargés respectivement avec 280 g et 1 080 g 
d’arsenic. 


— RECHERCHE DU RAYONNEMENT DIFFUSÉ PAR 
RÉSONANCE SUR LE NIVEAU DE 402 KEV.— Aucune 
différence significative n’a été constatée dans la 
région correspondant à un pic photoélectrique 
de 402 keV, entre les spectres d’impulsions obtenus 
avec les diffuseurs épais de As ou de Se, et ceci 
quelle que soit la source étudiée. Avec un facteur 
d’embranchement de 0,14 pour le rayonnement 
de 402 keV, les limites inférieures de la période 
partielle de la transition y directe, et de la période 
du niveau sont respectivement de 1071 $ et 
1,4 10-15. Avec une source de 20mC, Mertzger a 
fixé la limite inférieure de la période du niveau 
à 10710 s [8] : on peut remarquer que pour amé- 
liorer cette limite par la même méthode, il serait 
nécessaire à la fois d'augmenter l’activité de la 
source et d'opérer sous faible pression ; l’influence 
des collisions devient en effet très importante pour 
des périodes de l’ordre ou supérieures à 107125 et 
particulièrement lorsque la vitesse de résonance est 
voisine de la vitesse de recul (ici vr/Vx — 0,87). 


— RAYONNEMENT DIFFUSÉ PAR RÉSONANCE SUR 
LES NIVEAUX DE 205 ET 280 KE V. — Nous avons 
représenté figure 2a les spectres d’impulsions 
obtenus avec les diffuseurs minces d’arsenic et de 
sélénium, et la source n° 4, La raie observée avec le 
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diffuseur d’arsenic a été comparée avec celles pro- 
duites par le rayonnement direct de sources de "5e 
et de 2%Hg (fig. 2b). On peut en conclure que 90 0/0 


2 


x: Ÿ Direct 


o : Ÿ Résonance 


500 


300 Kev 


200 250 


Fic. 2b. 


sinon la totalité de l’effet observé doit être attribué 
à la diffusion par résonance du rayonnement 
de 265 keV. Dans la suite, la raie de résonance 
d'intensité VX sera entièrement attribuée au 
niveau de 265 keV. 

38 0}, des photons de 280 keV [9] sont associés 
au niveau isomérique de 0,017 s, et ne peuvent 
donner lieu à aucun effet. Compte tenu du rapport 
d'intensité des rayonnements de 280 et 265 keV 
(voir schéma fig. 7) on obtient t, (280) > 44, (265). 

Les résultats obtenus avec des ampoules conte- 
nant une quantité de plus en plus faible de SeH, 
font apparaître une augmentation très importante 
du rapport du nombre de photons diffusés par 
résonance au nombre de photons émis par la 
source ; un maximum parait être atteint pour les 
pressions les plus faibles. 

Nous avons représenté figure 3 les valeurs expéri- 
mentales de VN,/N, en fonction de la pression de 
remplissage des ampoules. La précision statistique 
devient mauvaise pour les mesures 7 et 8 corres- 
pondant à 2,1 et 1,2 cm de mercure. Toutefois 
il est probable que la diminution observée avec Ja 
source [8] est réelle et doit être attribuée à la 
décomposition partielle de SeH,, favorisée par la 
présence d’impuretés sur les parois ; une telle 
décomposition ne se produit en effet que très len- 


N°10 


tement lorsque la quantité de SeH, est plus élevée : 
ainsi la proportion de SeH, décomposé au bout 
d’un mois est inférieure à 15 0}, pour la source n°6, 
d’après l'intensité de la raie de résonance. Par 
ailleurs nous nous sommes assurés que la présence 
sur le verre des ampoules d’un dépôt de sélénium 
correspondant à la destruction de 10 04 de SeH, 
initial serait facilement visible, sauf pour les 
ampoules 7 et surtout 8. Enfin les valeurs portées 
dans la figure 3 ont été obtenues avec des sources 
fraîchement préparées. 


D. Période du niveau de 265 keV. — L’embran- 
chement de la capture électronique directe vers le 
niveau de 265 keV a été trouvé inférieur à 5 
(voir l’appendice). Dans ces conditions nous avons 
interprété les résultats en considérant que dans 
tous les cas le niveau de 265 keV était atteint par 
l'intermédiaire de celui de 402 keV. 

Nous avons utilisé pour la détermination de la 
section efficace de diffusion résonnante la valeur 
extrapolée NAN, = (1,25 & 0,15) 107$; la période 
du niveau de 265 keV alors obtenue est trop courte 
pour que le ralentissement des noyaux de recul 
par collisions sur les molécules du gaz intervienne, 
même dans l’ampoule (1). Par contre nous avons 
pu rendre compte correctement de la variation 
de NA/N, observée, en tenant compte de la période 
du niveau de 402 keV. 

Le nombre total de photons de 265 keV émis 
par les différentes sources a été déterminé en 
plaçant celles-ci à une grande distance du cristal, 
non protégé cette fois par le plomb. On observe la 
raie photoélectrique complexe des rayonnements 
de 280 et 265 keV, dont le rapport d'intensité a été 
pris égal à 0,48. 

Pour chaque élément du diffuseur de hauteur dh, 
le nombre de photons de.265 keV reçus varie selon 
l’angle d'incidence, de même que la masse super- 
ficielle effective ; il faut tenir compte également de 
Pabsorption du rayonnement incident dans le dif- 
fuseur (en moyenne 8 2}, pour les effets photoélec- 
triques, Compton et pour l’absorption par réso- 


nance), et de la variation du rendement photo- 


électrique du compteur selon la position du point 
d'émission des photons diffusés. 

Ce dernier effet a été étudié à l’aide d’une source 
de *#Hg, placée au-dessus du cristal dans une 
position standard, puis en différents points des 
parois du diffuseur ; une correction (5 °/) est néces- 
saire pour appliquer les résultats obtenus au cas du 
rayonnement diffusé de 265 keV. 

La position du diffuseur a été choisie de telle 
sorte que l’angle de diffusion moyen soit voisin 
de 1250. Dans ces conditions la section efficace 
moyenne est pratiquement la même que l’on 
admette une distribution isotrope du rayonnement 
diffusé ou une distribution en (1 + A4, P, (cos? 6)). 

Le niveau de 265 keV se désexcite dans 96 04 des 
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cas par émission y directe (4 0}, représentant la 
conversion interne et l’embranchement vers le 
niveau de 200 keV). 

En admettant pour les spins du niveau fonda- 
mental et du niveau excité la valeur 3/2, on obtient 
pour la période partielle et pour la période du 
niveau : 


A0 30 ES Ut — (2,8 0,210 Le 


la précision est principalement limitée par la 
mesure de la raie de résonance. 


Discussion. — Les périodes déduites par Metzger 
de la mesure de la section efficace de diffusion 
résonnante d’une part, et par la méthode d’absorp- 
tion propre d’autre part, sont respectivement 
D 'OMOEseL 165 1071ES. 

La première valeur est en très bon accord avec 
celle que nous avons obtenue ; le désaccord avec la 
seconde, d’ailleurs à la limite des précisions statis- 
tiques, ne devrait pas à notre avis être attribué à 
un effet résiduel des collisions dans le gaz. Une 
valeur limite de 1.8 1076 pour NZ/N, (fig. 3) nous 
paraît en effet tout à fait inadmissible, même en 
tenant compte de la décomposition possible 
de SeH, dans les ampoules 8 et 7; la variation 


de NA/N, avec la pression du gaz s’écarterait de 
plus fortement de la courbe théorique, même en 
ajustant au mieux la valeur de la période du 
niveau de 402 keV. 

Les valeurs de t, déduites de la section efficace 
directe, font intervenir l'énergie de la transition de 
capture électronique, déterminée par l’étude de la 
réaction >As(p, n)°%5e [7]. La vitesse V; de recul 
a été calculée à partir d’une masse de 77, les 
énergies de dissociation de SeH, ou *AsH, étant 
très supérieures à l’énergie de recul [10]. 

Enfin ces valeurs de t dépendent du rapport des 
largeurs l,/l, mais il est exclu que le niveau 
de 265 keV possède un embranchement supplé- 
mentaire de plus de quelques pour cent. 

On peut remarquer par contre que le coefficient 
d'absorption propre de la raie de résonance est 
proportionnel à l,/A, ; la largeur Dôppler A, due à 
Pagitation thermique dans le diffuseur est calculée 
pour un gaz parfait à une température effective 
légèrement différente de la température réelle [11]. 


E. Effets des collisions sur 5, (265 keV); 
période du niveau de 402 keV. — Le libre parcours 
moyen des molécules SeH, est 1/n, x d?si n, est le 
nombre de molécules par em% à la pression P, 
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et d le diamètre de ces molécules supposées sphé- 
riques. Avec une valeur de 2,6 10% cm pour ce 
diamètre et une vitesse moyenne de 2 105 cm}s 
pour les noyaux de recul, il apparaît que même 
dans l’ampoule 1 la probabilité de choc 
NS AT UE 

reste petite devant la constante de désintégration 
du niveau de 265 keV. La diminution de WN,/N, 
observée doit par conséquent être associée à l’exis- 
tence d’une période relativement longue du niveau 
de 402 keV. 

L'interprétation quantitative de la variation 
de NA,/N, en fonction de la pression (c’est-à-dire du 
libre parcours moyen L,) permet d’atteindre la 
valeur de cette période. 

Les noyaux de recul lors de la désexcitation du 
niveau de 402 keV sont animés non plus d’une 
vitesse constante V,, mais de vitesse V4, V;, V», 
Vs, etc. suivant que les radicaux auxquels ils 
appartiennent ont subi 0, 1 2, 3 chocs. Naturel- 
lement les vitesses V, — V, — V; peuvent prendre 
toutes les valeurs entre 0 et V3. 

Le nombre n{V.,) dV, de noyaux de recul appar- 
tenant à des molécules ayant subi un choc et un 
seul le ramenant à une vitesse V, avant la désexci- 
tation du niveau de 402 keV, dépend de la distri- 
bution des vitesses après le choc, des probabilités 
de collisions avant et après le choc, ainsi que de la 
constante de désintégration À. n(V,) et n(Vs) ... 
s’obtiennent en sommant les contributions pour les 
différentes vitesses possibles après le ou les chocs 
intermédiaires. En admettant pour la distribution 
des vitesses après un choc l'expression 


PIVOT er ar 
(collisions entre sphères rigides identiques), l, étant 
le libre parcours moyen et , le nombre de *AsH, 
formés, on a en définitive 

N=N x Ny [(Vx 3 1 Ny) 
pour le nombre de noyaux se désexcitant avant la 
première collision et : 


do M en Dray 
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Ve + \b Vr + |? 
Re ni 


les termes entre crochets se rapportant succes- 
sivement à 21(V), no V), n(V). 

Enfin les expressions (7) dans lesquelles on rem- 
place V; par la vitesse réelle V permettent de 
calculer, compte tenu du spectre de vitesse pré- 
cédant le rapport P(v})/P;(v3), soit 


Nr IN Nr IN) 
P;(v2) et NÂ se rapportent au cas où les collisions 
sont négligeables). 
Nous avons représenté figure 4a, en fonction du 
paramètre a — N/Vz, les contributions respec- 
tives des molécules ayant subi 0, 1, 2, 3 chocs, 
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+ 


rapportées à P;(v»), ainsi que P(v»)/P;(v») ; la 
contribution du quatrième choc, estimée pour 
PVR) 
PA) 


0,5 


0 0,5 1 1,5 
FiG. 4a. 


a —= 0,05 et 0,1 est négligeable au delà. Nous avons 
pris V'/V4 égal à 0,29 et v,/Vz égal à 0,565 d’après 
les énergies de la capture, du premier et du second 
rayonnement y. La figure 4b représente les mêmes 
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contributions, mais pour une valeur de @ 
donnée (0,05) lorsque v varie autour de ?3. 

Nous avons estimé, à partir de la dérivée seconde 
de P(v), la correction sur la valeur calculée 
de Nz/N;m due à l'effet de la température de la 
source et du diffuseur. Lorsque a est suffisamment 
petit pour que les molécules ayant subies plus 
d’une collision aient une contribution importante 
(voir par exemple figure 4b) cette correction peut 
ne pas être négligeable ; ainsi elle atteint 1500 
pour a — 0,05. Nous avons tenu compte de cet 
effet pour la comparaison avec l'expérience. 
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RÉSULTATS ; PÉRIODE DU NIVEAU 402 KEV, — 
Le paramètre a s'exprime explicitement en fonc- 
tion de {, période du niveau (en seconde) et de la 
pression P (en em de mercure) par la formule 


la masse spécifique de SeH, dans les conditions 

normales, la masse moléculaire et la section de 

collision ont été pris respectivement égales à 

3,63 mg/cem°, 81 g et 2,1 10715 cm2. 

Pour différentes valeurs possibles de Nm et 
différentes valeurs de {, on calcule ainsi une série 
de courbes de V,/N, en fonction de P. On peut 
constater sur la figure 3 que la courbe tracée 
pour MumiN, = 1,25 1076 et & — 1,5 1079 s est 
en bon accord avec les résultats expérimentaux 
sauf pour la pression la plus faible. Pour les deux 
autres courbes tracées à titre d'exemple et recou- 
pant la première pour les pressions moyennes, il y a 
décalage systématique de tous les points expéri- 
mentaux, dans un sens pour les pressions plus 
faibles, dans l’autre pour les pressions plus fortes. 
Des différents essais il est possible de conclure : 
t = (1,5 Æ 0,4) 10795. 

Discussion. —  L’existence d’un embran- 
chement de capture directe vers le niveau 
de 265 keV conduirait à une valeur plus élevée de 
la période t ; cet embranchement a été trouvé infé- 
rieur à 5 0/o. 

Les différents isotopes du sélénium sont présents 
dans le gaz actif et les collisions ne se produisent 
donc pas entre des molécules de masses stric- 
tement égales ; de plus ces molécules ne sont pas 
absolument sphériques, bien qu’on puisse estimer 
que la déformation due à la présence des atomes 
d'hydrogène n’occupe qu’une fraction faible de 
leur surface. L’énergie mise en jeu lors d’une choc 

est insuffisante pour dissocier le radical *AsH, ; 
bien que ce radical libre soit beaucoup plus actif 
qu’une molécule normale, la probabilité de réaction 
lors d’un petit nombre de collisions reste négli- 
geable. 

La principale cause d’erreur systématique sus- 


Exev Texp (par s) E; 
575 [2,8 1010 (F,)] 
402 6,5 107 2,1 1014 
304 3,6 
280 [1 X 10° (E,)] — < 810° 
265 0 134010 — [</10 (E:;)] 
200 [(E2) 2 108] 
136 2,5 108 81012. 
121 9,7407 4,7 1012 
97 2,7 407 
23 36 
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ceptible de modifier la valeur obtenue pour la 
période, provient du choix de la section de collision. 
Le rayon de la sphère de choc de SeH, ou *AsH, 
a été pris égal à 1,3 107$ cm, soit légèrement 
supérieur au rayon de l'atome de sélénium 
(1,17 10-8 cm) [12] et au rayon de covalence de 
la liaison avec l'hydrogène, 1,17 108 cm également 
[13]. 

Une limite supérieure de 2.10 $s pour la période 
du niveau de 402 keV a été fixée par Day [14] par 
la méthode des coïncidences différées. Metzger a 
par ailleurs observé comme nous-mêmes l'effet des 
collisions sur la section efficace de diffusion par 
résonance du rayonnement de 265 keV. Toutefois, 
disposant de sources contenant à la fois SeH, et H;, 
il a fixé seulement une limite inférieure pour la 
période (3 10-19s), 


F. Conclusions. — Les résultats que nous avons 
obtenus permettent de fixer d’une part la proba- 
bilité de transition du rayonnement y de 265 keV 
et la limite supérieure de celle du rayonnement 
de 280 keV, d’autre part les probabilités des diffé- 
rentes transitions y par lesquelles se désexcite le 
niveau de 402 keV. Nous avons pris pour les em- 
branchements les valeurs indiquées sur le schéma 
de la figure 7. 

Par ailleurs la période du niveau isomérique de 
304 keV a été déterminée par Vegors et Axel [15] 
puis par Schardt [9] ; les probabilités partielles £, 
pour les niveaux de 200, 280, 575 keV peuvent être 
déduites des valeurs de B(E,) obtenues par l’exci- 
tation coulombienne [16]. 

Nous avons rassemblé dans le tableau I l’en- 
semble des résultats ainsi que les probabilités théo- 
riques de transition (transition d’un proton 
impair) [17]. 

Le caractère M, du rayonnement de 265 keV 
semble bien établi, notamment par la distribution 
angulaire des photons diffusés par résonance, et par 
la nouvelle valeur du coefficient de conversion 
interne [8]. Le rapport 7, /T M, égal à 0,06 est 
près de 10 fois plus élevé que pour les autres tran- 
sitions M, observées dans les noyaux à nombre 
impair de protons [18]. 


M, E, M, F8 
3109 
2,21012 5,3 108 2,2 107 
5 105 180 
9 107 
521011 6,6 107 
1,6 107 
851019 2,35 105 105 
51010 1,3 105 5,5 104 
2,51010 4 X 105 4,4 104 
1,8 105 30 2 10—5 
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Le rayonnement de 280 keV est caractérisé 
comme un mélange M, + Æ; entre des ni- 
veaux 5/2 et 3/2—, les mesures de corrélation 
angulaire de [19], [20], [24] donnant pour le rap- 
port Æ,/M; de à 0,55 (si l’état du niveau 
de 402 keV est 5/2 +). La limite supérieure obtenue 
pour Tex ainsi que ni Fur de Tex É2s0nt compar 
tibles avec de tels rapports. 

Les caractères Æ, et M, respectivement ont 
seuls pu être attribués aux rayonnements de 304 
et 23 keV issus du niveau isomérique [9]. 

L'attribution d’un spin 5/2 et de la parité posi- 
tive au niveau de 402 keV soulève certaines diffi- 
cultés. 

Cette attribution s'appuie d’une part sur la 
valeur de log ft de la transition de capture élec- 
tronique, indiquant une transition permise, d’autre 
part sur les expériences de corrélations angu- 
laires [19], [20] et de corrélation associée à la 
polarisation [20] pour les cascades 121-280 keV et 
136-265 keV. Enfin la période partielle du rayon- 
nement de 97 keV est en accord convenable avec 
celle d’une transition Æ, (5/2* — 9/2+). 

Cependant les valeurs expérimentales [20] des 
coefficients de conversion des rayonnements de 121, 
136 et 402 keV sont plus de deux fois supérieurs 
aux valeurs attendues pour des transitions Æ;, et 


Deus 
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les périodes sont respectivement 6.104, 3.104 et 
3.106 fois trop longues. 

En attribuant au niveau de 402 keV une parité 
négative, les périodes des rayonnements de 121 et 
136 keV pourraient s’interpréter par un mélange 
M, + E, encore que les rapports £,/M, possibles 
soient difficiles à concilier avec le résultat des expé- 
riences de corrélation angulaire et les coefficients de 
conversion. Par contre, la période du rayonnement 
de 402 keV est supérieure à celle d’une transi- 
tion Æ, et celle du rayonnement de 97 keV 600 fois 
plus courte que celle d’une transition Y,. Or, sauf 
pour les transitions Æ,, les périodes expérimen- 
tales connues ne sont jamais inférieures aux 
périodes théoriques. 

On peut se rendre compte que le choix d’un 
état (3/2 —) ou (7/2 +) pour le niveau de 402 keV 
présente des inconvénients encore plus graves. Dans 
ces conditions il semble bien que le niveau ait un 
spin 5/2 et la parité +. On peut espérer confirmer 
encore une telle attribution en mcsurant à nouveau 
les coefficients de conversion des rayonnements 
de 121, 136 et 402 keV (les mesures existantes [19] 
ne sont pas sûres en valeur absolue). D’après les 
probabilités de transition obtenues, les proportions 
de M, dans les mélanges Æ£, + M, sont proba- 
blement négligeables dans le cas des rayonnements 


5000 
| 401Kev 
1S6Kev Spectre experimental 
| —— pics photoëlectriques corrigés 
4000 (401 kev direct et effet compton) 
3000 
2000 
280Kev 
| 304Kev 
1000 & 
Gi 
(0) 


200 


Fic. 


400 Kev 


No 10 


de 121 et 136 keV, mais non dans celui de 402 keV. 
Une autre cause susceptible de compliquer l’inter- 

 prétation des valeurs expérimentales des coef- 
ficients de conversion est l'influence sur ces coef- 
ficients Æ, des dimensions et de la structure du 
noyau [21]. 

L'existence de facteurs d’interdictions aussi 

élevés que 107 est bien connu dans le cas de tran- 
sitions Æ, observées dans les noyaux déformés. 
Pour des noyaux de À moyen on ne s'attend pas, 
d’une manière générale à l’existence de tran- 
sitions Æ, associés à des niveaux de faible énergie. 
Aussi, dans le cas de As est-il étonnant qu’un 
niveau 5/2*+ se présente pour une énergie d’exci- 
tation de 402 keV. Le modèle de Nilsson [22] 
pourrait expliquer pour un noyau déformé pos- 
sédant 33 protons une succession de niveaux 9/2*, 
1127, 9/27. 

I] n'existe pas à notre connaissance d’autres cas 
de transitions Æ, dans les noyaux moyens, dont le 
caractère puisse être considéré comme suffisam- 
ment bien établi. 


Appendice. — Schéma de désintégration de ’°Se. 
Compléments. — 7’Se a été l’objet de plusieurs 
études récentes [19], [23], [24], [9], [20]. Les inten- 
sités relatives des principaux rayonnements 7 
paraissent bien établies mais non les valeurs des 
coefficients de conversion interne. Ainsi, d’après 
les mesures préliminaires de Metzger, le coefficient 
de conversion À pour la transition de 265 keV 
serait voisin de 0,0064 au lieu de 0,016 obtenu par 
Schardt et Welker [19] ; toutes les valeurs données 
par ces auteurs seraient alors environ 2,5 fois trop 
élevées. 

Les études de corrélation angulaire pour la cas- 
cade 136-265 keV, conduisent à une anisotropie 
nulle ou très faible ; pour la cascade 121-280, le 
coefficient À, obtenu est — 0,36 [19], — 0,41 [24], 
— 0,435 [20]. Par contre, alors que Schardt [9] ne 
signale pas d’embranchement de capturc directe 
vers les niveaux de 265 et 280 keV, Van den Bold 
et coll. [20] évaluent ces embranchements à 24 0 
au total. 


LIMITE SUPÉRIEURE DE L'EMBRANCHEMENT DE 
CAPTURE VERS LE NIVEAU DE 205 KEV. — Nous 
avons repris l’étude de Se par la technique du 
cristal-puits, précédemment utilisée par Lu et coll. 
[23] pour mettre en évidence les principales cas- 
cades, et par Schardt pour l’étude du niveau 
isomérique. 

Nous avons représenté figure 5 le spectre d’im- 
pulsions obtenu avec une source de Se placée au 
centre d’un cristal de INa — T1 de 1” 3/4 de dia- 
mètre et de 2” de hauteur. (Le puits est creusé dans 
l’axe jusqu’à une profondeur de 1” 1/4.) La raie 
très intense observée correspond pour une part au 
pic photoélectrique des photons de 402 keV, mais 
surtout aux superpositions dues à la désexcitation 
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du niveau pour des rayonnements y en coïncidence. 
Les autres raies doivent être successivement attri- 
buées : 

— au rayonnement y de 97 keV aboutissant au 
niveau isomérique de 304 keV ({ — 0,018 5) ; 

— aux rayonnements de 121 et 136 keV — cette 
raie complexe est observée du fait de la désexci- 
tation des niveaux de 280 et 265 keV par conversion 
interne, ou par l'émission de photons qui ne sont 
pas détectés dans le cristal ; 

— trois processus peuvent rentrer en jeu pour 
rendre compte de l’existence d’une raie photo- 
électrique au voisinage de 280 keV : a) la détection 
des photons de 304 keV et surtout de 280 keV pro- 
venant du niveau isomérique ; b) celle des photons 
de 265 et 280 keV qui se trouvent en cocincidence, 
soit avec les électrons de conversion, soit avec des 
photons non absorbés correspondant ‘aux tran- 
sitions de 136 et 121 keV ; c) celle de photons 
de 265 ou 280 keV associés à des embranchements 
de capture directe. 

Nous avons calculé lefficacité du cristal et 
mesuré expérimentalement le rapport du pic photo- 
électrique à la surface totale, en fonction de l’éner- 
gie. Compte tenu du rapport d'intensité des pho- 
tons de 402 keV à ceux de 265 et 280 keV, il est 
alors possible de déduire du spectre d’impulsion 
total la contribution de la transition directe. 
L’intensité du « pic somme » ainsi obtenue est 
relativement peu sensible à la valeur de l’embran- 
chement de capture directe admis «a priort, et à 
celles des coefficients de conversion ; on peut 
remarquer que le rapport de cette intensité à celle 
de la raie complexe 121-136 keV est en bon accord 
avec la valeur calculée dans le cas du schéma de 
Schardt. 

La surface des raies associées aux processus (4) 
et (b) a été déduite, la première de’ celle de la raie 
de 97 keV, la seconde de celle du « pie somme ». En 
admettant les coefficients de conversion de Schardt 
et Welker, la somme de ces deux effets se trouve 
être nettement supérieure à la raie expérimentale. 
Avec des coefficients de l’ordre de 2,5 fois plus 
petits (voir plus haut) et compte tenu de l’impré- 
cision sur la détermination des surfaces, le pro- 
cessus (c) pourrait représenter au maximum 400/o 
de l'effet observé. Dans ces conditions la somme 
des embranchements de capture vers les niveaux 
de 265 et 280 keV est inférieure à 5 0/o. 

Par étude de la diffusion inélastique des neutrons 
sur l’arsenic Lind et Day [25] trouvent également 
que lé rayonnement y de désexcitation du niveau 
de 402 keV est indistinguable de celui d’une source 
de 75e. 


EXxISTENCE D'UN EMBRANCHEMENT DE 0,04 0/6 
VERS LE NIVEAU DE 579 KEV DE As. — Temmer 
et Heydenbourg [16] ont mis en évidence par exci- 
tation coulombienne des niveaux de 7°As à 200 et 
280 keV mais aussi à 573 et 811 keV. 
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Nous avons repris l'étude du rayonnement y 


de *5$e en interposant devant le cristal d'INa — TI, 


un écran de plomb de 2 cm, absorbant les photons 
de faible énergie. Nous avons observé dans le 
spectre d’impulsions (fig. 6) une raie produite par 


1577 Kav 


40kev filtré par 
2cm de Pb 


S00 


un rayonnement y de faible intensité et dont 
l’énergie a été trouvée égale à (577 F 2) keV. 
Cette valeur est en bon accord avec l’énergie du 
niveau observé par excitation coulombienne. De 
plus, le mode de préparation de la source 
(SeH, gazeux) exclut la présence d’impuretés 
radioactives autres que As. Pendant deux mois 
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&ueune variation du rapport d'intensité de la raie 
de 577 keV à celle de 402 keV n’a été observée. Ce 
rapport permet d'estimer à 0,040), l’embran- 
chement de capture vers le niveau de 577 keV, en 
admettant qu'il n’y a pas d’autres transitions que 
la transition directe vers le niveau fondamental 
(Temmer et Heydenbourg ne signalent pas d’autres 
modes de désexcitation.) , 


75 
É VA 127J. Se 
Kev 2 
700! 
200 € -470 
7 
PT 
702 577 004% 
+ 
400 % | 402 Kev 
| 97 a 4136 57% | 
| 4 | 
300 72 304 Kev 
LA 280Kev 
VA 265 Kev 
200 VA 199 Kev 
100! 
1 75 
0 2 AS STABLE 


Free 


On obtient pour log (ft) de la capture vers le 
niveau de 575 keV la valeur de 9, correspondant 
à une transition une fois interdite. 


Le schéma de désintégration de 7°Se finalement 


adopté a été représenté figure 7. 
Manuscrit reçu le 11 juillet 1958. 
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CONDUCTIBILITÉ THERMIQUE DES DIÉLECTRIQUES A BASSE TEMPÉRATURE 


Par P. RÉSIBOIS, 
Université Libre de Bruxelles (Belgique). 


Résumé. — Nous présentons une objection à l’argument de Peierls sur le rôle spécial joué par 
les « Umklapprozess » dans la conductibilité thermique à basse température : par définition même, 
il semble que la quantité J — E fn; soit seulement définie à un vecteur du réseau réciproque 
près ; nous ne pouvons donc pas donner « & priori» un rôle particulier à ces U. P. 

Nous discutons l’expression du flux thermique et nous montrons que la condition nécessaire 
pour obtenir une conductibilité thermique finie est l'existence de chocs à-trois phonons de fré- 
quences w, w’ et w” telles que w’ et w” restent finis quand © — 0 ; les-U. P. forment une classe 


particulière de tels processus. 


Abstract. — We raise an objection against Peierls’ argument concerning the special role played : 
by « Umklapprozesse » in thermal conduction at low temperature : by the very definition of wave 
vectors, it seems that the quantity J — E %f n; is only defined « modulo » a reciprocal lattice 
vector ; we may then no more give « a priori » a particular importance to U. P. Discussing the 
expression for the thermal flow, we show that the condition to obtain a finite thermal conduc- 
tivity is the existence of three phonons collisions of frequencies w, w’ and w”, in which w’ and w” 
remain finite when © — 0 ; U. P. are a particular kind of such processes. 


1. Introduction. — La conductibilité thermique 
des diélectriques dépend d’une manière complexe 
de la température. En effet, tant que nous nous 
limitons au modèle du cristal « harmonique » idéal, 
les quanta de vibration, ou phonons, se propagent 
librement sans atténuation. Pour rendre compte de 
la résistance thermique, il faut considérer des effets 
d'ordre supérieur comme : 10 l’anharmonicité du 
potentiel d'interaction entre les atomes du cristal, 
ce qui produit des « chocs » entre phonons ; 20 la 
présence d’impuretés (y compris les isotopes) ou 
d’imperfections qui diffusent ces phonons ; 39 les 
dimensions finies du cristal. 

En fait, pour des cristaux assez purs de dimen- 
sions normales, le premier effet est prépondérant 
tant que la température n’est pas trop basse. Nous 
nous limitons ici à la région de température située 


dans les limites approximatives 10 < - < 20 (où 0 


est la température caractéristique de Debye de 
l’échantillon), pour laquelle la quantification est 
déjà importante tandis que les facteurs 20 et 30 ne 
jouent pas encore un grand rôle. 

La formulation du problème due à Peierls [1] 
consiste à considérer les phonons comme des cor- 
puscules plus ou moins localisables (ce qu’on peut 
justifier en formant des paquets d’onde de pho- 
nons) et à écrire une « équation de Boltzmann » 
d'équilibre stationnaire entre la variation du 
nombre de phonons d’espèce X (k = (f, s) ; Î : vec- 
teur d’onde ; s: polarisation) due au flux sous 

A 
l’action du gradient de température, soit . 
et celle due aux chocs entre phonons, soit a - 

* [chocs 


1 
flux 


ANe 


dt 


flux dt 


= 0, (1) 


chocs 


Le premier terme de cette expression se calcule 
facilement [1]; à des termes d'ordre supérieur 
en 7 /0x : 


dx 
dt 


+1 ME (ae |) 


Lt 
flux LR 


où «4 est la fréquence du phonon k, 
DE 
 — 1) donne la distribution 


d'équilibre, 

(ve E donne la composante, selon la direc- 
tion x du gradient de température, de la vitesse de: 
groupe des phonons k. 

Le terme de choc peut se calculer également si 
l’on admet qu’il est le même que dans un système 
à température uniforme (hypothèse qui semble 
correcte dans le cas des gradients faibles) ; il est 
alors donné par l’équation de Peierls, quand on se 
limite aux termes cubiques dans le potentiel 
d'interaction : 


a 
dt flux 


= — 2 LEE R") © o or À NE(NE + 1) (NE +1) 
k'k" 


dgr 


= Né(NE +1) 


flux 


(gx — gr — gr) D(ox — ox — or) 
+ 2UNE +1) (NE + 1) Né(ge + gx — 8x) 
Sox + ox — ex) }. (3) 
Dans cette expression, nous avons introduit la 
fonction g, pour mesurer l'écart à l'équilibre : 
Nr — NN = ne = NUÛNe +1) &; (4) 
L{kk' k") est une constante de couplage qui est 


| 
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nulle sauf pour les processus qui satisfont la rela- 


tion de conservation des vécteurs d’onde à un 


vecteur du réseau réciproque, soit K, près : 


SE MS EE (5) 


le signe devant f’ est le même que celui qui inter- 


_ vient dans la conservation de l'énergie : 


off, s) F o(f’, 5°) — off”, s”) = 0 (6) 
assurée grâce aux fonctions delta de Dirac intro- 
duites dans (3). 

On appelle les processus pour lesquels K — 0 
! des processus normaux (P. N.) et ceux avec K 0 


des Umklapprozesse (U. P.). 


2. Rôle « a priori » des Umklapprozesse. — 
L’équation de stationnarité obtenue grâce à (1), 
(2) et (3) est très difficile à résoudre. Aussi Peierls, 
et après lui d’autres auteurs [2, 3], ont toujours, 
pour obtenir malgré tout une solution qualitative 
du problème, étudié la grandeur suivante : : 


2 h fe ri (7) 
(ni : nombre de phonons k à l’équilibre stationnaire) 


définie à partir des vecteurs d’onde f. 

Rappelons le raisonnement de Peierls : J, reste 
invariant au cours d’une collision du type normal ; 
seuls les U. P. sont susceptibles de le modifier ; 
si on réalise une situation initiale avec J, 0, 


les U. P. seront indispensables pour amener le 


système à l'équilibre thermodynamique, pour 
lequel on doit avoir J, — 0. Ayant établi ainsi le 
rôle fondamental des U. P., on en déduit faci- 
lement que la conductivité thermique doit être de 
la forme : 


x = Tr'exp. es — 2): (8) 

* Il semble cependant qu'on puisse élever la cri- 
tique suivante à cet argument : on sait que les 
vecteurs d’onde f (nous traitons le cas à une dimen- 
sion mais la généralisation est immédiate) s’intro- 
duisent lorsqu'on effectue la transformation cano- 
nique qui permet de passer des coordonnées indi- 
viduelles des atomes, soient (x,, p,) aux coor- 
données normales (g, ÿ;); on pose à cet effet la 
relation : 


jé MU > exp.ufan.qj(t) (9) 


c'est-à-dire que f intervient uniquement dans 
Pargument d’une exponentielle imaginaire. Il en 
résulte qu’on ne change rien quand on ajoute à fun 
multiple entier de 27/a. Autrement dit, aucune 
grandeur ayant une signification physique ne peut 
dépendre de la valeur exacte de f ; elle ne peut 
dépendre que de f. (modulo: 27 Ja). 

C’est seulement ensuite, puisqu'il n’y a que N 


CONDUCTIBILITÉ THERMIQUE 


DES DIÉLECTRIQUES 797 
coordonnées normales indépendantes, qu'il s'avère 
commode de définir les valeurs possibles de f dans 


a 
11 paraît clair qu'on ne peut pas attribuer de sens 
physique quelconque à une grandeur en se basant 
sur cette convention commode. Si l’on veut dis- 
cuter la grandeur J,, on est obligé de la définir de 
la manière suivante : 


la première zone de Brillouin (5 aus 5). 


Pa = Era (fe 217) Dares 


Dans ces conditions, s’il est bien évident qu’un 
P. N. conserve J,invariant, un U.P.nele modifie pas 
non plus, à cause de l’indétermination modulo27 /a. 
Dès lors les U. P. ne jouent plus, «4 priori, aucun 
rôle fondamental dans l’amortissement du flux 
thermique. 


3. Conductibilité thermique dans un modèle 
simple. — Comme on ne voit plus très bien quelle 
signification physique intéressante on peut attri- 
buer à J,, il est bien préférable de considérer le 
flux thermique lui-même, soit : 


. | 


comme la vitesse de groupe n’est pas la même 
pour les trois phonons qui entrent en collision, 
®, n’a, malgré les relations de conservation (5) 
et (6), aucune raison de rester invariant, ni dans 
un P. N,., ni dans un U. P. 

Nous allons rechercher la condition nécessaire 
pour obtenir une conductibilité thermique finie et 
nous montrerons, d’une manière approchée, qu’on 
doit s’attendre également à une relation du type (8). 

Supposons que l’équation d'évolution (5) puisse 
être remplacée approximativement par une expres- 
Sion : 


Dy = D À op ni À (10) 


dNx dnx 7x 


dé dt TL (A) 
où +4 est le temps de relaxation des phonons k. 
On a, grâce à (1), (2) et (11) : 
D PET [bale ce MAR 1 (12) 


RTS 


Dans l'expression du flux thermique (10), on 
passe d’une somme sur f à une intégrale en rem- 
plaçant le spectre de dispersion exact par un spectre 
de Debye, ©, = C;lf| (C, : vitesse du son dans la 
branche de polarisation s) ; on obtient alors faci- 
lement pour le coefficient de conductibilité ther- 
mique [2] : 


à ( 


13 
exp.x —1)? (13) 


où on à posé, pour simplifier, 


ñ © 


— — 14 
Ca TT (14) 
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Pour avoir une évaluation du temps de rela- 
xation r, il suffit de poser g’ = g” — 0 dans (3); 
ceci fournit un coefficient de conductibilité ther- 
mique «brownien», c’est-à-dire qu’on suppose que 
chaque mode normal se comporte comme s’il se 
trouvait dans un milieu complètement à l'équilibre 
thermodynamique. Nous discuterons cette approxi- 
mation plus loin. Il vient ainsi, passant de la 
somme sur f’ à une intégrale : 


(e2 — 1) (ex — 1) 


w’.&”.et lo + © — ®°] 
+2 rene A pt) US) 


Dans cette formule, nous avons supposé que le 
coefficient L(kk' k’) était approximativement indé- 
pendant de k, 4’, k”. De plus, pour f et f” donnés, 
f''et w(f”) sont entièrement déterminés par (5). 

À cause de la fonction à, l'intégration s’effectue 
sur la surface S de l’espace des f” qui satisfait les 
conditions de conservation (5) et (6). 

L'expression (15) croît rapidement avec © : en 
- effet dès que x est suffisamment grand les facteurs 
exponentiels se simplifient approximativement à 
cause de la conservation d'énergie (6) ; il reste un 
intégrant qui croît avec © et qu’il faut sommer sur 
une surface qui, elle aussi, augmente avec ©. Il en 
résulte que le temps de relaxation lui-même décroit 
très vite avec la fréquence ; dans un calcul grossier 
comme celui-ci, nous pourrons nous limiter à cal- 
culer + pour les petites valeurs de x, qui seules 
donneront des contributions importantes à (13). 

Or pour de faibles fréquences «, on peut dis- 
tinguer deux types de solutions qui satisfont les 
conditions de conservations (5) et (6) : 

19 Type I : w'et ©" — 0 quand © — 0. 

20 Type II : w’ et w” restent finis quand w — 0 
(pour la condition de résonance w” = © + «w' 
seulement). 

-On peut décomposer +=! en deux parties, corres- 
pondant à cette classification : 


4 1 1 


I TII 


(16) 


où trj'et zx: sont toujours donnés par (15) mais 
où les surfaces respectives sur lesquelles on intègre, 
soient S, et Sr, satisfont la première à la condi- 
tion 10, la seconde à la condition 2°. 

Pour « petit, le premier type de solution donné, 
puisque $, æ &?[4]: 

Lo TJ, af SOL Nes) 

(17) 


rm @T. 


On voit immédiatement que si ce type de solu- 
tion existait seul, on aurait, portant ce résultat 
dans (13), 


7 + Hé exp.x dx 
o (exp.z—1)? 
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qui diverge pour les basses fréquences. Ceci signifie 
qu'on aurait une conductivité thermique infinie 
parce que les phonons de grande longueur d’onde 
ne seraient pas suffisamment amortis. . 

Par contre, si les fréquences w’ et ©” sur les- 
quelles doivent porter l'intégration sont toutes 
finies, on peut poser Aw’ = how” —0/b(b > 1 mais 
fini). Il.en résulte pour xx! une expression de la 
forme : 


4 4 0 ) 
où /(0/T) est une fonction qui varie comme une 

certaine puissance de 0/7 ; en effet dans ce cas-ci, 
la surface sur laquelle se fait l'intégration ne 
dépend en général plus de « dans la limite 
& —> 0 [4]. Portant ce résultat dans (13), l’inté- 
grale ne diverge plus et on obtient un résultat de la 
forme. 


us À 


x Tnexp 0/bT. (19) 
4. Rôle de la loi de dispersion. — La formule (19) 
montre que la conductibilité thermique dépend 
exponentiellement de la température dans la région 
10 < : ce 20) , pour autant qu'il 
existe des processus du type IT ; dans le cas con- 
traire, notre modèle prévoit l’absence de toute 
résistance thermique. 

Le problème se pose donc de vérifier si les lois de 
dispersion usuelles fournissent effectivement de 
telles solutions. 

Nous discuterons brièvement les trois cas sui- 
vants : 

19 Approximation des milieux continus. 

20 Approximation de Brillouin-Debye. 

30 Lois de dispersion exacte, tenant compte de 
l’anisotropie cristalline. 


considérée 


19 Mizreux coNTINUS. — Un milieu continu 
peut être considéré comme le cas limite d’un milieu 
discret dont le nombre de degrés de libertés tend 
vers l'infini. On peut analyser ses vibrations en 
modes normaux (c’est-à-dire en série de Fourier, 
avec un nombre infini de termes) et on détermine 
les longueurs d’onde possibles en supposant les 
conditions aux limites périodiques usuelles : . 

_‘2T ne 


= — 00). (20) 


RO PRO TE 


La vitesse de propagation des ondes, pour une 
polarisation donnée, est indépendante de la lon- 
gueur d’onde ; on a donc la relation de dispersion 
simple : 

of) = Cf]. (21) 

En outre puisque toutes les longueurs d'onde, si 
petites soient-elles, ont un sens physique parfai- 
tement défini, la condition de conservation des 
vecteurs d’onde doit être exactement satisfaite, 


f 
| 


étudiés ont la même direction ; 
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c’est-à-dire que (5) doit être remplacée par la for- 


mule plus restrictive ; 


fr f—1 = 0. (57) 
Pour rechercher qualitativement les fréquences 
qui satisfont les conditions de résonance (5') et (6), 


nous utilisons la construction graphique de 


Peierls [1] dans le cas où les trois modes normaux 
rappelons briè- 
vement cette méthode : dans un plan («, f), on 
indique sur la branche s le point P correspondant 
au mode initial (f, s) (fig. 1) ; à partir de ce point 


0 QP' re 


Fic. 4. — Milieu continu : 
Processus normal : f = OP’; f = PQ’ ; f” = 


OQ”. 


comme origine on retrace la branche s’ (en pointillé) 
correspondant au second mode ; son intersection Q 


avec les lignes originales donne les solutions cher- 


chées, comme on le vérifie facilement. Dans le cas 
des milieux continus qui nous occupe nous rejetons 
toutes les solutions qui correspondent à un simple 
glissement d’une branche sur elle-même ; il y en a 


. une infinité mais elles ne peuvent physiquement 


correspondre à aucun processus réel car il suffirait 
un écart infinitésimal à la loi de dispersion (21) 
pour qu’elles disparaissent toutes. Dès lors, puisque 
les vitesses C, sont toutes du même ordre de 
grandeur, on voit sur la figure que quand © — 0 
(c’est-à-dire quand P tend vers l’origine), les seules 
solutions possibles sont du type I. 

L’approximation des milieux continus conduit 
donc à une conductivité thermique infinie. 


20 APPROXIMATION DE BRILLOUIN-DEBYE. 
Comme il est difficile de calculer exactement la loi 
de dispersion d’un milieu discontinu, on utilise 
souvent ce modèle qui repose sur les simplifi- 
cations suivantes : 

a) on néglige la dispersion et on suppose le 
milieu isotrope, c’est-à-dire qu’on conserve la rela- 
tion simple (21) ; 

b) on introduit un « cut-off » qui respecte Je 


nombre total N de vibrations propres de chaque 


branche en remplaçant la surface qui limite la 
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première zone de Brillouin de l’espace des f par une 
sphère de rayon fn». 

Nous pouvons généraliser légèrement ce modèle 
en conservant en outré la notion de vecteur du 
réseau réciproque : nous remplaçons simplement 
chaque vecteur K par un rayon-vecteur de lon- 
gueur 2fpm et dont la direction est celle de K. 

Les conditions de résonance (5) et (6) s’écrivent 
alors : 


2 DE 2{D:m réciproque (U. 12) 
Î FT Î’— 1 — à 7 


(P. N.) 
Ce Gt Cr 0 


Répétons la construction que nous avons décrite 
au paragraphe précédent en choisissant pour sim- 
plifier une direction dans laquelle se trouve un 
vecteur 2fn mrée. ; NOUS Voyons (fig. 2) que pour 


om 0 


Q P° f fo,m 


Fic. 2, — Milieu de Brillouin-Debye : 
Exemple de collision : (s = 1) + (s = 2) — 
Processus normal : f = OP’ ; jf’ — Aie d 
Umklapprocess : f = OP’; jf = LRO 


& —> 0, en plus de solutions du type I, nous obte- 
nons aussi des solutions du type II correspondant à 
des U. P., du moins si le phonon appartient à une 
branche transversale. Le modèle de Brillouin- 
Debye conduit donc à une conductibilité thermique 
finie pour les branches transversales, grâce à 
l'existence de U. P. ; par contre, pour la branche 
longitudinale, la résistance thermique est nulle, du 
moins tant qu’on se limite aux effets d’anharmo- 
nicité. 


39 Lor DE DISPERSION EXACTE. — Analysant les 
propriétés des cristaux réels, Herring [4] a pu 
montrer que l’anisotropie cristalline introduisait 
nécessairement, pour une large catégorie de sys- 
tèmes cristallographiques, des processus N du 
type IT ; en plus, les U. P. existent évidemment 
toujours. L'existence de tels P. N. du type IT est 
liée à des dégénérescences des fréquences, c’est-à- 
dire qu’il y a des surfaces ou des lignes de l’espace 
des f pour lesquelles les courbes w, — w,(f) corres- 
pondant à des branches différentes se recoupent 
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C'est ce que nous avons représenté très schéma- 
liquement ci-contre dans la construction de Peierls 
à une dimension ; on voit facilement que 
quand «© —> 0, on peut trouver ©’ et &” finis. Il 


s’ oO: RU 


Fic. 3. — Milieu anisotrope : 
P. N. du type II : OQ’ et OR. 
U. P. du type II : OS”. 


serait très difficile de calculer explicitement laquelle 
des deux sortes de solutions, P. N. ou U. P., est la 
plus importante mais nous sommes cependant 
assurés que pour chaque branche la condition 
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nécessaire pour que x soit fini peut être satisfaite, 
du moins pour la grande majorité des cristaux. 


En conclusion, il apparaît que ce sont les pho- 
nons à basse fréquence qui sont prépondérants 
dans le transport de l’énergie thermique ; la condi- 
tion générale à remplir pour amortir suffisamment 
ces phonons est qu'il existe des processus dans 
lesquels un phonon de basse fréquence « entre en 
collision avec un phonon «’ pour donner un nou- 
veau quantum «”, de telle sorte que quand © — 0, 
w'.et ©” restent finis. Les U. P. forment un cas 
particulier de tels processus. 

Remarquons encore que puisque les chocs impor- 
tants se font avec des phonons de haute énergie 
qui sont beaucoup plus près de l’équilibre que ceux 
à basse fréquence, l'hypothèse brownienne que 
nous avons faite plus haut s'avère satisfaisante. 


Je tiens à exprimer ma reconnaissance à M. le 
Professeur Prigogine pour les conseils qu'il ma 
prodigués durant la réalisation de cette recherche. 
Je remercie aussi M. Philippot pour d’intéres- 
santes discussions. 


Manuscrit reçu le 18 juin 1958, 


BIBLIOGRAPHIE 


[1] Peters, Ann. Physik, 1929, 8, 1055. 
Quantum Theory of Solids, Oxford, 1955. 
[2] KLEMENS, Proc. Roy. Soc., À, 1951, 208, 108. 


[3] Laesrrien et ScnLômAnn, Nachr. ak. Wiss. Gott., 
1954, 4, 71. 
[41 HERRING, Phys. Rev., 1954, 95, 954. 


LE JOURNAL DE PHYSIQUE ET LE RADIUM 


TOME 19, OCTOBRE 1958, PAGE 781, 


ÉNERGIE DE VIBRATION-ROTATION DES MOLÉCULES POLYATOMIQUES. 
TABLES DES COFFFICIENTS DE L'HAMILTONIEN TRANSFORMÉ DU SECOND ORDRE. 


Par M. L. GRENIER-BESSON, G. AMAT et S. MAES, 


Laboratoire d’Infrarouge. 
Laboratoire de Chimie Physique de la Faculté des Sciences de Paris. 


Résumé. — On rappélle le principe du calcul, au second ordre d’approximation, de l'énergie de 
vibration-rotation d’une molécule polyatomique et on donne les tables des coefficients des hamil- 
toniens transformés du second ordre k;, hs, h£!. Ces tables permettent une étude systématique plus 
commode des effets d’ordre deux et des corrections d’ordre supérieur. 


Abstract. — After recalling the principle of the second order computation of the rotation- 
vibration energy of a polyatomic molecule, the authors give tables for the coefficients of the second 


order transformed Hamiltonians A, hs, h3!. 


These tables enable second order effects and 


higher order corrections to be studied more easily. 


L'énergie de vibration-rotation Er d’une molé- 
cule est donnée par les valeurs propres de l’hamil- 
tonien de vibration-rotation H dont l'expression a 
été établie par Darling et Dennison [1] (voir Appen- 
dice). 

Cet opérateur étant trop complexe pour que ses 
valeurs propres puissent être obtenues de façon 
rigoureuse, on effectue un calcul de perturbation, 
après avoir développé Æ en série par rapport aux 
- coordonnées normales : 


SH, HU. (1) 


Shaffer, Nielsen et Thomas [2] ont proposé 
d’effectuer ce calcul de perturbation par la méthode 
de Van Vleck [3], [4], [5], [6], c’est-à-dire de sou- 
mettre l’hamiltonien à une transformation de 
contact : 


HO TAT = LR LR +... (2) 


choisie de telle sorte que A -+ À soit diagonal par 
rapport aux nombres quantiques vibrationnels v,, 
Vo, +.) Ds, --., nombres quantiques principaux 
relatifs aux diverses vibrations normales. Les élé- 
ments matriciels de À; qui ne sont pas diagonaux 
par rapport à tous les nombres quantiques v,, 
Do, -.., D ... n'interviennent pas dans un calcul 
de Eyr au second ordre d’approximation ; on 
pourra donc ne considérer dans un tel calcul, que les 
éléments matriciels de À + ki + k;, diagonaux 
par rapport aux v,. 

D'autre part, l'opérateur Æ et par conséquent les 
opérateurs 6, k; et À; sont diagonaux par rapport 
auxnombres quantiquesrotationnels J'et M. Ilreste 
donc, pour obtenir l'énergie Eye, au second ordre 
d’approximation, à diagonaliser À + k + h; par 
rapport aux autres nombres quantiques définis- 
sant les fonctions de base, à savoir : 1° les nombres 
quantiques vibrationnels (1) /, et m, ; 20 le nombre 
quantique rotationnel Æ. 


(1) A des vibrations harmoniques non dégénérées, dou- 
blement dégénérées, triplement dégénérées correspondent 
des fonctions d’onde définies respectivement par un (2), 
deux (v,, 1) ou trois (vs, {s, ms) nombres quantiques. 


Les valeurs Eye seront obtenues en résolvant des 
équations du type (?) : 


det { SAT D de mec 
Me Eki deb PJ M0 pe ms 


OR 0 ee Om me SU URI) 2010 (9) 


où Ô désigne le symbole de Kronecker. 

Les opérateurs k et ki sont donnés par les for- 
mules (4) et (5). (Les notations sont celles de 
H. H. Nielsen [7], [9]; la définition des divers 
symboles utilisés est rappelée dans l’appendice.) 


AI ET LE NS (ee \ 
DRM Ne 1/2 | Pso 2 
or) 2 > Le 2 2 Âs »2 As ss) (4) 
k: Ag > > > CES [se Pa 5 
; Aie Csoso” (dso Pso/ Qso’ Pso) PRIS (5) 
æ 8 oc, Lou 


O<o 


La forme analytique de l’opérateur , dont 
l’étude constitue l’objet principal de cet article, 
sera donnée dans le prochain paragraphe. Les pro- 
priétés des opérateurs h$ et k; peuvent être schéma- 
tisées comme suit : 


OPÉRATEUR 4  OPÉRATEUR h; 


Molécules à symétrie 


axiale diagonal diagonal 
Molécules linéaires diagonal nul 
Molécules à symétrie 
sphérique (°) diagonal non diagonal 
en K et m" 
Molécules asymétriques non diagonal nul 
en K 


I résulte de ces propriétés et de l’équation géné- 


(2) Des éléments matriciels non diagonaux par rapport 
aux nombres quantiques v, interviennent dans le calcul 
lorsque des résonances se produisent entre des niveaux 
possédant accidentellement des énergies voisines : réso- 
nances anharmoniques (résonances de Fermi, de Darling- 
Dennison), résonances de Coriolis [8], [9], [10]. 

(8) Par « molécule à symétrie sphérique », nous enten- 
dons une molécule dont l’ellipsoïde d’inertie est une sphère. 


56 
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rale (3) que l'énergie £yr peut être développée sous 
la forme : 


Eyr = Eo + E1 + E2 + .…. (6) 


et que les termes successifs du développement (6) 
sont, pour les divers types de molécules, donnés 
par les équations (7), (8), (9) et (10) (*). 


a) Molécules à symétrie axiale : 


(7a) 
Bo UD ML ve RON SORTE 
(7b) 
Bi = (LR, M, NUL LUE MC EME EN UE ee) 
(7c) 
det.{ (7, Æ, M, 210 le 0e DEC M 0 ER] 
+ Beer. dv ER 0: 
b) Molécules linéaires : 

(8a) 
Eos (J, M, 22208 lc RE PPT AI EE 
E, = 0 (8b) 
(8c) 
det LM, "de ls ce NE I M RUN CS) 

en 3 al 0. 

c) Molécules à symétrie sphérique : 

Es = UT RME ns Nm ee 
ROUE Me Eine Eee) Ra a 

dt ATAR, MO RS Lee 

| LA Je de M, CALE Us, Le, Ms; EX: 3) 

Le See 2: 0mymete Es | UOb) 

det! (J,K, M, :2.08 de; lbips ee 


FAR TS ATEN ue Lions) 
— dxx ... div, Omym'y -:: (E1 + E>)} = 0. (90) 


d) Molécules asymétriques : 


dot PAS KM ame te IEC EN NEA) 

— xx Es} —0 (104) 
Eu 10 3 (10b) 
det{ (RM, Ds. h ET IR UNE EEE 


— rx/(Eo + Es) } — 0. (10c) 


Le calcul de l'énergie de vibration-rotation au 
second ordre d’approximation a été effectué par 
Nielsen [7], [8], [9]. 


(*) En écrivant les équations (7), (8), (9) (10), on a tenu 
compte des faits suivants : 1) le nombre quantique Æ 
n'intervient pas dans le cas des molécules linéaires ; 2) les 
nombres quantiques m, n’ont de signification que dans le 
cas des molécules à symétrie sphérique, ces molécules étant 
seules susceptibles de posséder des vibrations triplement 
dégénérées : 3) le nombre quantique /, n’a pas de sens dans 
le cas des molécules asymétriques qui ne possèdent que 
des vibrations non dégénérées. 
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L'énergie d'ordre zéro Æ, apparaît comme Ja 
somme des énergies d’un rotateur rigide et d’un 
ensemble d’oscillateurs harmoniques. Les correc- 
tions d’ordre un et deux correspondent aux effets 
suivants : 


— Ordre un  : Interaction de Coriolis. 
Résonance de Fermi. 


Résonance de Coriolis. 


— Ordre deux : Termes  anharmoniques de 
l’énergie de vibration. 

 Distorsion centrifuge. 

Väriation des moments d’iner- 
tie en fonction des nombres 
quantiques vibrationnels. 

Résonances et dédoublements 
du type L. 

Résonance de Darling-Denni- 
son. 

Résonance de Coriolis. 


0 


Hamiitonien transformé d’ordre deux. — Les 
corrections du second ordre à l’énergie de vibra- 
tion-rotation sont liées à la valeur des éléments 
matriciels de l'opérateur A. Dans le cadre d’une 


_série d’articles consacrés aux énergies de vibration- 


rotation d'ordre élevé [11], [12], [13], Amat, 
Goldsmith et Nielsen ont proposé pour h; une for- 
mule très générale particulièrement adaptée à une 
étude systématique des effets du second ordre : 


hs, = ZE Y Pa Pe Py Po 
+2 ZE Y® pa Pa Pe Py 
+È Z (EYE pa pu + € Yab ga 1) Pa Pe 
a< + 
+2 2 S Ÿ%% pa Pv Pe Pa (u) 


a<b<c 


1 
A 2 Z @ Vas (da 9 Pe + Pe Ga 4) Pa 
a<b < 


Por 
bcd 


a 
a<b;c<d 
Eee 


abcd 
a<b<ce<d 


1 
2) Va (da 9b Pe Pa + Pe Pa Ga 4) 


(2) For Ja 4b de Ga- 


Dans cette expression les indices a, b, c, sont 

utilisés pour simplifier l'écriture, à la place de so, 
/ / L/4 (4 

Dans la référence [12], A., G. et N. ont également 
défini un opérateur k; et un opérateur k; donnés 
par des formules analogues à l’équation (11) et 
utilisés respectivement dans le calcul des correc- 
tions d’ordre 3 et 4 à l’énergie de vibration-rotation. 
Le coefficient de l’opérateur 


(1/2) (ab ... Pe Pa +. + Pe Pa -:. Ga gb...) Pe PR... 


| 
| 
% 
| 
Î 
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TABLE I 
my 15 ah 
s Br os Lis 
il AIS, — X) | 
(2 y ya roi Àa M a I. 
Ex 8 r1/2 Te 18 19 2 (tas 77 jet OS a%$) a (u Es x) 
(Sn 5) 
RS D (at cet) œes 
AES A 1 1 cs am am Am LEZ 
noi) 
HEC 
APTE  Aun a 
@ PES Ponte STE (1 + dom + dbm) 
Pat > (Ts de SE Uno ) 1 + Ds ae Ssmsa 
2A(1 + dab) m I, 16 MN Nice 
(sy F8) (CE) 
___ Pa)? Car (Z —;.) SF 
Aa — À) 14 Up y (28 
a aë 
a Y 4 ñ to 8 48 ag af + ab = 
POP Be ee MU + RL ‘loge To ; 
7 en? © Kabm 2% 1 Bab, 
o > Je Fe (1 Sam + Dom) a BA 4 
ñ Âa dr: Xon D Ày “es Xm 
5e 41, 180% Ay) 1/4 (1 ES ab) 2 (Cam a à Com (a) ee re . (1 “JE dsysm) A ee BE . (1 + Ssn5m) 
CEE So) (8, n À 3p) 
RAS Res 7 ss 
2()a — Xy) (a M)/4 1, U I, 6252.) 
à yabc HS A “2 : Kimi mi Go x -. eu de dy À g Oyni) (1 F Ôs,8) 


mn) tab 


Ka Ve: KA > au? At 4(ag se aa + aÿ te) 
(2) € fr TS Xp) /4(4 se dab) 


: ; 
LC TE NE £ en UP 
re DL : abm HS ED) (eee Le _ : ) 
FR œ m ke 07 
TC À 
ie ARTS im Com Kame Sue (1 AE ae + Î dem) (1 TE Var) (1 Ge Dsnse) 


(Q) ñ N 
Le Cam Kbme = sue (l Te beat dem) (1 tt dy) (LEE ÿ 


te 
— RC > FT Kabm PBab,e(1 #6 Van sf Spm) 
m œ 
3m S$ 
(rs) 


=>, DAC (0 de) 
LU 1) 741 + ab) (1 + ed) 


(2) 


= F2 RANCE dam + Obm) (1 + Sem + Dam) 
m 


2 Ÿ'abea = he kavea + F Ch À À kgu kmns Paille Des + Dj) (1 + Dms + Onj). 
nn) “abod) ; 
k< 
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TABLE II 
aByST _ 1 AS 
© 12 m An Lo Jp Ty 18 1 
séyTa ag) Aa!/ (3m — Aa) | 
mi = Sr A + | 
sy EE 
1 6 | 
Se ENEN PE PU A 7 (na NE Al sE 2 Es 9 Bas 
Sn ‘a 
no ) LR 
8 Tab mr C »à Kabm adm ae ES, 1 
@ SRE Le le m4 am or 
1+ Sm$b QE Bsnsa 
Pa UE ss (Cam Cno + mb Com io 2 
3A(1 + dub) m Lol Xa — ln Àÿ — Xm 
Sn AS) (A5) 
Cha An)! ce Eu Ce 
Abu de Ve il | 
£. A 
2) * ab 1/4 Aa A + CR 
21% Ip(%a 6) (1 + av) Q T5 
m chu? Kabrm # ; 1 GBab.m 
DE (À Ban + Bon) (a Ga 
Hans 2 (tn té + Un Cm) do. >) (: +98) + (7) (a + ni 
61e 20e AJ + Bo) me 0 PR 2 / Li — am) \ 262 
* (C7) (En Æ5) 


RER en en re 

ED — 0) De AT, y Tor 

œ 'T'abc 2m ce * An 1/4 Baja 
Cr Tam F2 : kimÿ mi (ee 7,14 (1 + Gi + Bms) (s A, “) 


fun) =(abc) 
(<m) 


SM GN( + dE) 


& 
( = 
MT EU tle re AA + Dub) 
DRE CT Le x ore (1 + Sam + dom) (ur Ets ee pe 
3 .m Te G UP TC RALRE TE Os à es 
mn FSe 


2x c S Km! 


ne 3(1 + db) m 1e 


ô 
ar SN ST SES (: es 


8 RATS 
R cn Rome 2 (4 Le Ôpe + dem) (1 + Sym) | 1 <: Pen 


Ap173 < 2 
TC 2Cme 
TO Ù SAT a Kabm GBab,el1 as Sam 2 Spn) 
æ œ 
( La So 
om Fe 


Ta _ 5 Me) Vt (Ce Cha + Via Lie) 
PUR Tale A) (EE dus] (1 + des) 


2 h 
ne TE 2 Kavm Kcam cam + dam + 5pm) A Sem + Sam) 


: 2x ch 
@T aboa = hC Kabea + = F2 2 kyts Æmns Ball + des + Du) (1 + mi +: Du). 
Lee tr à Ç 
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TABLE III 
LI Lee 1 ap rs 
M nl TT 
1 af\ Na 4(3Xn me) 
(4 net (4 
él 7 16%2 1, Te 7h > (sur + Lim dm) An (Ka — Ar) 
(Sn a) 
1 BY 1 1 Cr 
rs DEUEN ITU 2 Ed Rule FAT LE Dev 
8m a 
Sn a 
£ ê 
aBJJab ÆC Kobm adm am 
SU AXE Es Fe Ie NUE (1 | Dam A pm) 
dsm5b dsnse | 
1 . DE SRS 
(a A5)? (4 de ce Cnb En) “ 3 3 
kan ù JET Ka — ÀÂm bp —X 
ne + do 
Àp)1/4 1 dl 
____ Ua ) = | - I #4) 
3404 — 5) 1, My AN CET 
ñ 8 8 D a ap? 2 ap aëe 
aB [] | ES A NEA A% . @ | 
nn tee El SH Ts 
me T _. D abn Am Kabn au 


1 Bov,m 

RU, PTE (1 Se am LE Sm) Ps sl 
ñ 

814 Ig(a ho) 7/41 


: \ s À De Dsas 
Z (tèm Com + Com Cam) Le “ 2) (1 + Le (er) (: # 1 
n dub) nm Àa Km À EX À 3 - 
(sn #8,) (Sn A8) 
di ’ 
Mille ru). ae (Te is) 
6(4 = À) (X 1) / JT; 2 1 VASE 
be 
a Ur? Te à 2 kymi An es = 4 sn ie dv + mi) (: ns Pa) 
RAI 5 XjU 3 
Cp) tab 
au, 3 ME 20" Ce + dE Ce) 
PA 6 Lo mAl 


Et Jab) 
TC D Ce Xel/4 Kavm 


Bob 
m I (1 + dom ne Om) ( ar 
OR : 


TC Am L# 
F 2(1 an dab) im Te 


Cèm Et a ere 2 (1 Ne dac QE 0) (1 10 dam) (: airs 


) 


OR 8 
+ Lim Koné” eu + de + dom) (1 + Jo) (: + “ «)! 
DM NET-2C en à 
ER — Kabm Baboll + dam + db) 
STE E 
Eu) 


ges 3 Dead (EE Ch + it) 
PO © Ta A) (À + Duo) À + Bed) 


ÉE - 0 Kabrn keam egm(i + am 5 dy n) (IREE dem =) Sam) 
m 

$ 2 
Üavea — he kavoa + + è 


> : Ki Km GB, s(L TE dk te dj) (1 + Ôn + ni) 
(ktmn)= (abcd) 
(&<L):(m <n)_ 
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: dans ke est désigné 


qui est ni Dar ore 
dans . 


par: 76-1dans ot ié', BAUER 


Tables I, II, III. — Nous donnons dans les 
tables I, IT et III la valeur explicite des coefficients: 
DA date 0 Ten 
et 
28-05 Vo: 


Les formules données dans ces trois tables sont 
valables d’une façon tout à fait générale, pour tous 
les types de molécules et quelles que soient les 
valeurs particulières des indices &, 8, y, ... et 
AP brC 

La signification de la plupart des symboles figu- 
rant dans les tables est rappelée dans l’appendice. 
Nous compléterons la définition des notations uti- 
lisées par les remarques suivantes : 

19 Au lieu de caractériser les diverses coor- 
données normales par des indices so, s'a', 

"0", ..., nous préférerons souvent les caractériser 
(ainsi que cela a déjà été fait dans l’équation (11)) 
par des indices a, b, c, 

Chacune des lettres a, D, . représente alors en 
fait deux indices que nous écrirons Sa Ou, Sb Oh, --. 

Un symbole de Kronecker à,4 est égal à l’unité 
lorsqu'on a simultanément 


S5 = 1Sh MELUN ON 


Les grandeurs affectées d’indices a, b, c, ... sont 


de deux types : 
a) certaines dépendent à la fois des indices 54, 5%, 


5, LL'et des indices oo cu last A4 


Cab, Kabe; Kabca) 
b) les autres dépendent seulement des indices s,, 
Sp Se see (as love ave) ; NOUS les écrirons indif- 
féremment 


so ds OÙ A; abc; PBab,c- 


GB 
a Sb Se Sa Sb, #e 


20 ss et PB, désignent deux fonctions de À, 
À, À définies par les formules ci-dessous : 


Aave — 2: 2] 1/2 1/2 1/2 - 1/2 1/2 
Ja + À: + Àc En DV Ha 0: En me 
1 1 
HR NP, NRA, MB | 
WE 1 il 
a RE Lr ÉNÉENE CRE 
1 j 
+ UCI 72 dE VERS XP L CRE | 


qui peuvent aussi s’écrire sous la forme suivante : 
Fe Aave = 20 6 Xe)? Dabe (Sa À So F Sc) 


(Sa Su) 


a aab = 


PHYSIQUE > No 10 
/ Babe ES REX — 4 — 5) D (Sa À So À Sc) 
294 — À 
Borb — ne (55 "Se 
1/2 
Baba Æ x (Sa  Sb) 
1 
À CCE Taxe 
avec Dove : î 
Ne 
He Qu FEVER 7) (Ur x J'AI) 
1 
Xe - 
CHR Re 5 Ou? + M — x?) 


39 L'écriture 


‘à et 2 
rar Rte) ne (abcd) 
(<m) k<lm<EN 


indique une sommation sur toutes les permutations 
distinctes des trois indices lmn ou des quatre 
indices klmn, compte tenu des conditions imposées 
à ces indices. 

Sinous posons, par exemple : 


1/4 
ja 


œ J 
(2ylman Te — I, r Kimi Em Con vu 


X (1 He di #È ni) (1 25 ds, Sn 


la table I donne : 


AVR 


mn 
(Imn)= (abc) 
1<m) 


œylm,n 
(214 


ce qui signifie : 


@ fes = (yabe + Dyaeb + ia 
a<b<C) 

® Fa = (yat + (jyab.a 

à ns = (yat + (Hybba 

es (que (ya, 


4° Pour simplifier l'écriture, le symbole Z, a ebe 
utilisé à la place de 79. 

50 Comme nous le Dee dans l’appendice, 
chacun des indices «, f, y, ... peut désigner l’un 
quelconque des axes 2yz, axes principaux d'inertie 
de la molécule. De plus, les indices primés y” et y” 
sont définis à partir de l’indice y de telle sorte que 
la suite Yy'Y” soit identique à ZyZ OU à yZX ou 


1, 


à zxy. (Il en sera de même pour les indices xxx 


86e"). 
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60 La quantité Ya est donnée par la formule : 


ù 1/2 C4 } , 
n xBYaa — Par ( CH ) | au ) 
FOI ENTI AIRES PE EN PES 
D ) (x) 
To Ter Iy xp JL, Te IL, /4=6 Tale 


ne NL 


ay’ 
LR Tao | ( ) ( | 
_ T6 ne. Mlgrls/g er Le le I, pr 


dans laquelle le premier terme est nul si «  Ê 
le second terme est nul si «” =£ B, etc. Cette for- 
mule générale prend une forme plus simple dès que 
l’on tient compte des valeurs relatives des indices «, 


= YY F 4 7 
4 Yo Cie Ga CG =) +% (z ) 
OEM Tel: 

Si a ÆBEÉY 


2 “ANGERS SFlLa=sÆrz 
PRES LS ATOME Si CAE SAC ER 


NOTE ue ( à RP de 
241015 y 15 ë LE. 


= 0 SI 


Explication des coefficients. %-:; Y:::. — Pour 
calculer pratiquement les diverses corrections du 
second ordre à l'énergie de vibration-rotation, il est 
commode d'écrire les formules du tableau I sous 
une forme plus explicite obtenue de la façon sui- 
vante : 

— pour chaque coefficient %:;; Yé-::, on dis- 
tingue plusieurs cas correspondant aux diverses 
relations d'égalité pouvant exister entre les 
indices Sa Ga, Sy OC, Se Ge, ASSOCIÉÉS AUX 4, b, C, .…. 

Pour le coefficient #$ Y®, par exemple, on écrira 
trois formules donnant respectivement : 


D D «sr. 
ee) (Sa 4) 
Sa 7=0b 
— On développe, dans les formules, les 


sommes, *»,*», > en tenant compte de ce que 
m mn Hklmn j ; À 
pour une molécule d’un type donné, certains coef- 


ficients Æymn Sont nécessairement nuls : 

19 Pour une molécule à symétrie axiale ou une 
molécule linéaire (dont nous désignerons les 
vibrations non dégénérées par s, s', s”, ... et les 
vibrations doublement dégénérées par £, to, di, te, 
6, ...) on peut dresser le tableau suivant : 


(5) « = 8 Æ y £ 0 signifie que, dans le calcul de % **vza 
par exemple, il faut prendre 3 = 3 (et à = x, à” = y). 
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TABLEAU A 
COEFFICIENTS COEFFICIENTS POUVANT 
NÉCESSAIREMENT ÊTRE DIFFÉRENTS 
NULS DE ZÉRO (6) 

Kaaa kttt, k totata Kiss 

Ksst, Ksst, ksss’ 
Kaav Ktit Ktitte kits Kttas 

titil'a totat'1 

ktsse ktss’ Rss’s” 
Kave knts kits kits 

ktit kit kits 

ktit ka 


Le tableau À peut être résumé ainsi : 

— lorsque un ou trois des indices !, m, n, corres- 
pondent à des vibrations doublement dégénérées, 
Krmn eSt nécessairement nul, 

— lorsque zéro ou deux des indices /, m, n, 
correspondent à des vibrations doublement dégé- 
nérées, Æmn peut être différent de zéro (toutefois, 
kits A 0). 

20 Une molécule asymétrique ne possède pas de 
vibrations doublement dégénérées, les seuls coef- 
ficients km susceptibles d’être différents de zéro 
sont donc 


kss, ksss’, Kssrer. 

39 Le cas des molécules à symétrie sphérique est 
plus complexe par suite de l’existence de vibrations 
triplement dégénérées et nous ne l’envisagerons pas 
ici. Il importe d’ailleurs de noter que la plupart des 
molécules à symétrie sphérique sont du type XY, 
tétraédral. Une formulation générale présente donc 
pour ces molécules un intérêt beaucoup plus 
restreint que pour les molécules linéaires, asymé- 


triques ou à symétrie axiale. 


Table IV. — Les coefficients %::;Y%:-, explicités 
comme il vient d’être dit, sont donnés dans la pre- 
mière colonne de la Table IV. On a admis que les 
seuls Æmn Susceptibles d’être différents de zéro 
étaient ceux qui figurent dans la dernière colonne 
du tableau À. Il en résulte que la Table IV est 
valable pour toutes les molécules à l’exclusion des 
molécules à symétrie sphérique, pour lesquelles des 
termes supplémentaires devraient être introduits 
dans les formules. 

Dans les deuxième et troisième colonnes de la 
Table IV figurent les facteurs par lesquels il faut 
multiplier chaque terme du cœfficient %:; Yé::: 


pour obtenir 8-47: et 8; U%:: respectivement. 


La suite de la table prend un développement 
trop important pour être publiée dans le cadre 
de cet article. Élle a été entièrement établie, 


(S) Le fait que les coefficients Æyris, ttss, kut,s peuvent 
être différents de zéro a été souligné récemment par 
F, L. Keller, A. H. Nielsen et W, H. Shaffer [14]. 
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TABLE IV 
(2) £ (2) T (2) U | 
eBy8 ÿ 1 ae a 2 1 
PEUT TEL en UE 3 3 À 
1/4 £ 
ru? [ci «By Aa (3m — da) 2 1 
AA MORE PERRET NDS US MGR MES NET 2 1 
() x 81, Te 1 _ (on m + Cam | NA RATE 3 5 
(on #5 
PATE 1 Le ST RQ \ 
RL 
+ 4e Ip (ing) am M su m tL s ; 
Om Fa 9 3 
+ noyau | | 0 
1/2 
ñ 8 Ya 
72 af ÿas Ds D œ 
2" (2) TE Ig = Cam Cam Fe 4 Fo 
(u#%) 2 { 
x : 2 
k chu? a$ a8 Àa 3 
+ k + D k : be 
La Te | er pos Fa AP A 
m (2 
SAUT ; 
ñ j 
n2. ë re D'ENCS En + Cm con) Li 9 1 
80 Sp 8 Te Ta m À — Xm = s 
Be (s,#s ) n) 2 
ñ 8 + 1 1 
2.8 yab _ D æœ GB œ ) ( XVE ( ] 
Los) I, TG We (Gén Cbm + Com Cam) (ha À5) TA ZAR He PRE 2 1 
Ris 3 DA 
Sn F 35 
A (a Xp) H4 \ 
FR D (tacle > (ete) 
c a 
(cxe,) (22) 1 1 
m [) a l 2 3 
ES A6) Co Lr — | 
a — À Tell, (HER ES) F 
x chi/2 32/4 ë/4 \ 
ST Er ag” 15 < + 2havb 4 = ? 
CRC) Xp (4 — Àp) ÀG (224 — da) 
2 1 
Te 23 Kkabm a Lu A) UE ae Dobn & 3 
m 
(2) } 
Sn ES à ] 


Au cours du récent séjour de l’un d’entre nous 
(G. A.) à Ohio State University, les questions 
traitées dans le présent article ont fait l’objet de 
fructueux entretiens avec le PT H. H. Nielsen ; il 
tient à l’en remercier très vivement, 


Appendice : Rappel des notations utilisées. 


1. OPÉRATEURS QUANTIQUES ASSOCIÉS AUX 
COORDONNÉES ET AUX MOMENTS : 


P, avec « = %, y, 2: composante du moment 


angulaire total suivant la direction d'un axe prin- 
cipal d'inertie «. 

Qsa : coordonnée normale. L'indice © caractérise 
les diverses composantes d’une même vibration dé- 
générée de fréquence co, = 1/2/2r. 

(6 =1,0o —1,2ouc = 1, 2, 3 suivant que la 
vibration est non dégénérée, doublement dégénérée 
ou triplement dégénérée). 

pr : Moment conjugué de Q,, : 

Pas = — iA D[dQs 
coordonnée normale sans dimensions 
eo = (M [?)48 Qsa 


es : 


No 10 


moment conjugué de 4 


Pso : 
Pas = — if D [dGs0. 


+ 2. PARAMÈTRES DÉFINISSANT LA CONFIGURATION 
D'ÉQUILIBRE DE LA MOLÉCULE : 

m,; : masse du ième noyau. 

a?, y?, 22 : coordonnées de la position d’équilibre 
. du ième noyau par rapport aux axes principaux 
d'inertie de la molécule. 

(a«ÿ, 89, ... désignent l’une quelconque des trois 
coordonnées précédentes). 

I@, : moment d'inertie à l’équilibre par rapport 
à l’axe principal «(x = %, 7,2) 


19 = > mi(Be + Yi ass 


3. COEFFICIENTS és Cars me, Aus Asset : 

Ts, Yi, & : COMpOsantes, suivant les axes prin- 
cipaux d'inertie, du «vecteur déplacement » du ème 
noyau à partir de sa position d'équilibre («;, B;, ... 
désignent l’une quelconque des trois composantes 
précédentes). 

2 : élément de matrice définissant la transfor- 
. mation linéaire permettant de passer des coor- 
données normales aux déplacements cartésiens pon- 
dérés : 

ml Ge » Léo Qso 


or UD, A7, A, coefficients définis, à partir 
des précédents, au moyen des formules suivantes : 


Con Ê= D (ee or ra De! 4] - (avec AY xyz, 
î 


YZX OÙ 3ty) 
HU AD (blem), (ep Æy 
à < 
GE = ZT mi (ef Do + 69 Ho), 
(aÆ8) î 


cel af œ 8 
A 050 — Asas’o — À soso” CS'a's"a" 
C) 
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4. FONCTION POTENTIELLE : 

V : Fonction potentielle, développée en série 
par rapport aux coordonnées normales sans dimen- 
SiOn : 

Roos os Ksos/o/# ss" 0" : coefficients, exprimés en 
cmt, des termes cubiques et quartiques dans le 
développement de V : 

ñ 1/2 2 
pi 9 > As ‘so + oi 2 2 ksos’o's"0" Qso s'0’ Qs’o 
s0<s/0°<s°0 
+ DE D >» » Ksos’0/s’0"s/"0!" Qso Gs'a’ Asa! so” + …. ): 


so so’ 8/0” NE a 
so<s'o/ <s”o” <s/’0”’ 


5. HAMILTONIEN DE VIBRATION-ROTATION (FOR- 
MULE DE DARLING-DENNISON) : 


1 
2 mix (Pa — Pa) Las 7 V#(Pe — pl ut + 


1 
+ QUAD pie ue pes uit + W 
= ke] 


Pa avec « —=%, y, £: composante du moment 
angulaire interne suivant la direction d’un axe prin- 
cipal d'inertie «. 


FÉOIb> 


so s’o” 


À 1/4 a 
cv Ssos’o’ so Ps'o’ 
s 


u : déterminant de la matrice { deg À. 
us : éléments de la matrice inverse d> 


1e TER RUE AU | 
re 1 Ta S To 


dont les éléments sont définis de la façon suivante : 


# ) à / 
Tux — 1% + >> dsg Oss + D » Asos'o’ Osa Qs'o 
sa sa so 


Añastat = 2 (ho stat + ha la), (x # B 4 1 Lg = —2 0% Qu — DE Asso Oo ss 
x Æû) sa se s’a 
LU EE HOT ( 3 
“(a AE Manuscrit reçu le 7 février 1958, 
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LETTRES À LA RÉDACTION. 


ACTION D'UN CHAMP MAGNÉTIQUE 
SUR LA SENSIBILITÉ PHOTOÉLECTRIQUE 
D’UNE CELLULE AU SÉLÉNIUM 


Par G. BLer et J. VINGENT, 
C. N. R. $., 66, rue Saint-Sébastien, Marseille. 


Une cellule à couche d’arrêt au sélénium, soumise à 
un champ magnétique, subit, lorsqu'elle est placée 
parallèlement aux lignes de force de ce champ, une 
diminution de son courant photoélectrique, indépen- 
dante du sens du champ. 

Nous avons étudié ce phénomène en fonction des 
paramètres suivants : 1° intensité lumineuse et lon- 
gueur d’onde du rayonnement incident, 2° intensité du 
champ magnétique. 


1. Dispositif expérimental. — La figure { représente 
schématiquement notre dispositif. 


Filtres F 


RICA. 


La lampe L a une tension stabilisée au moyen d’un 
montage en tampon. Les filtres F sont, soit des filtres 
interférentiels, soit des filtres Wratten Kodak, auxquels 
on ajoute jusqu’à 0,654 un filtre anti-calorique. C, 
et C, sont deux cellules de même type (cellules montées 
sur laiton) que le prisme P permet d’éclairer par le 
même rayonnement incident. L’électro-aimant E 
donne un champ magnétique A atteignant 
20 000 œrsteds. 


2. Principe des mesures. — La méthode utilisée est 
une méthode de zéro. Le courant photoélectrique débité 
par la cellule C, compense exactement, en l’absence de 
champ, le courant de la cellule C, : l'intensité lumi- 
neuse reçue par chaque cellule est réglée pour être la 
plus voisine, donc r est très petit. 

On mesure les variations Aï du courant photo- 


électrique à l’aide d’un galvanomètre G qui, à la dis- 
tance utilisée, a une sensibilité : 
_G = 0,47.10—% A/mm. 


Ce montage est tel que, même si la résistance r, du 
galvanomètre est assez élevée (ici r, — 100 Q), la cel- 
lule C, se comporte comme si elle était en court-circuit 
et le courant photoélectrique mesuré au moyen d’un 
micro-ampèremètre est le courant de court-circuit vrai 


Pre; 


En effet, pour une variation At du courant photo- 
électrique la différence de potentiel V aux bornes de la 
cellule C; qui était nulle en l’absence de champ 
deviendra : 


V=rhi. 
ceci pour un courant débité z (ici i © 10 A), C; débi- 
tera donc sur une résistance apparente r” 
A Pin CAQA 


T = — Torre 
L CE 


Dans les conditions de nos mesures, la valeur maxi- 


mum de Atfi est 6.10%, donc dans le cas le plus. 


défavorable : 
r = 100.6.10—% = 0,6 ohms. 


La cellule étant en court-circuit, sa courbe de 
réponse peut être considérée comme linéaire dans un 
très grand domaine d’éclairement et le niveau d’éclai- 
rement pourra être repéré au moyen du courant 
photoélectrique débité. 


Résultats de mesures. — A. MESURE A CHAMP 
CONSTANT. — 1) En faisant varier l’éclairement, on 


constate que Aï/i n’en dépend pas lorsque la compo- 


sition spectrale du rayonnement incident reste la 
même. 

2) En faisant varier la longueur d’onde du rayon- 
nement incident : 

Pour ces mesures, le faisceau lumineux sortant d’un 


N° 10 


monochromateur donnait une intensité lumineuse trop 
faible pour que Aifi soit mesurable et surtout iden- 
tifiable par rapport aux fluctuations. 

Dans le visible, les filtres interférentiels donnaient 
des bandes relativement étroites de longueur d’onde 
mais malgré le filtre anticalorique, laissaient passer 
une certaine proportion de rouge et de proche inira- 
rouge, variant de quelques % à quelque °/ du violet 

. à l’orange, Ceci nous a obligé à introduire un facteur 
de correction, le phénomène étant environ cinq fois 
plus grand aux environs de 0,74 u. 

Au delà de 7 000 
filtres Wratten avaient des bandes passantes beaucoup 
plus larges calculées en tenant compte de la répartition 

spectrale énergétique de la lampe, du coefficient de 
transmission des filtres et de la sensibilité de la cellule. 

Les résultats obtenus sont donnés (fig. 3) pour un 
champ Æ — 10 000 cœrsteds. 


5.10 | Az/ 


210° 


Grandeur maximum 
de l'écart autour de la 
valeur moyenne . 


10° 


ER XA 


8000 


4000 6000 


Fic. 3. — Variation de Az/i en fonction de la longueur 
d’onde pour A = 10 000 œrsteds. 


Étude de la courbe (fig. 3).— On obtient de 4000 À 
à 6000 À une courbe lentement variable avec la 
longueur d'onde à concavité tournée vers le bas. Au 
delà de 6 000 À, Ai/i croit très rapidement de 10-% 
à 5.10-3 et la courbe présente un maximum accusé 
aux environs de 7 400 À. Ce maximum dépassé, 
Aifi diminue lentement. Il est intéressant de noter la 
similitude entre cette courbe et celles obtenues par 
G. Blet dans l’étude des phénomènes de fatigue ou de 
réactivation aux basses températures [3]. 


B. MESURES AVEC CHAMP MAGNÉTIQUE VARIABLE 
(fig. 4). — Si on compare entre elles les courbes obte- 
nues pour différentes longueurs d’onde en les ramenant 
toutes à passer par un même point, on constate qu’à la 
précision des mesures, ces courbes se superposent. 
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D'une façon générale Aï/i croît moins vite que le 
champ appliqué et d'autant moins vite que le champ 
est plus élevé. 


25 | . 


EAN r 4 oersteds 


15000 


2 

3000 000 

Fig. 4. — Variation de Ax/i en fonction du champ magné- 
tique :4 À = 5 275 ; 5 825 ; 5 980.—2 > A — 5 9275 ; 

5 425. — 37 À — 4 575 ; 4 895 ; 5 075. — 4 — À — 

6 375. — 5 À — 8 600. — 6—>AÀ = 6 750. — 7 — À — 

8 000. — 8 — À — 7 000 ; 7 800. — 9 —>X — 7 400 en À. 


La superposition des courbes permet de conclure 
que les 2 variables 77 et À sont indépendantes. Aï/: est 
une fonction à deux variables séparées, Æ et à, et peut 
se mettre sous la forme : 


Aifi = f(H).80), 


Interprétation des résultats. — Les résultats: que 
nous avons obtenus montrent que le courant de court- 
circuit à — Ne (peu différent du courant de court- 
circuit réel) diminue avec le champ magnétique, il y a 
donc diminution du nombre ÆV de photo-électrons 
produits par lumen. Cette diminution est sélective par 
rapport à la longueur d'onde, elle doit par suite corres- 
pondre à une variation du rendement quantique p : 
nombre de photo-électrons arrachés par quantum inci- 


n n ; +. 
dent, p— 5 R v, pe nombre d'électrons émis par watt. 


Pour compléter les renseignements obtenus, nous 
avons étudié la cellule en circuit ouvert, à l’aide d’un 
électromètre ; nous avons constaté que la tension à 
vide, v — N'g(v) r, diminue avec le champ magnétique. 


ÜN'= NS ; 
$ — surface de la cellule ; 
N — nombre de photo-électrons par cm? ; 
r = résistance de la couche d’arrêt ; 


g(v) = n{N fraction d'électrons atteignant l’élec- 
trode antérieure. 
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Cette diminution peut être due à la diminution de 
l’un ou l’autre des facteurs N our. 

La résistance r, à la précision de nos mesures, ne varie 
pas, le nombre d’électrons de la bande de conduction 
n’est pas modifié alors qu’il y a variation du nombre 
d'électrons atteignant l’électrode antérieure de la 
cellule. Cette variation serait donc imputable à une 
modification du rendement quantique sous l’influence 
du champ magnétique, et cela plus particulièrement 
aux niveaux d'énergie responsables de la sensibilité 
spectrale dans le domaine des longueurs d’onde centré 
sur 0,74 pu. ; 

Cette sensibilité particulière aux champ magnétiques 
confirmerait l’épithète de métastable attribué par 
G. Blet à ces niveaux à la suite de ses recherches effec- 
tuées aux basses températures. à 


Conelusions. — L’effet que nous venons d’étudier 
peut être représenté par deux courbes, l’une donnant 
la variation de Aï/1 en fonction de Æ pour une longueur 
d’onde déterminée, l’autre donnant la variation de Aï/i 
avec la longueur d’onde pour un champ donné. 

Cet effet diffère de l’effet photo-magnéto-électrique 
longitudinal qui se produit pour les corps massifs et 
qui a déjà été étudié sur des monocristaux tels que le 
germanium et l’oxyde de cuivre. Dans la cellule, l'effet 
se produit dans la couche d’arrêt située entre la surface 
antérieure de la couche de sélénium et la surface 
métallique transparente, d’autre part le sélénium étant 
formé d’un agglomérat de microcristaux est isotrope 
alors que tous les cristaux étaient anisotropes. 


Lettre reçue le 31 mai 1958. 
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CONTRIBUTION 
A L'ÉTUDE DES PHOTOMULTIPLICATEURS 


Par P. Cacnon et A. SARAZIN, 
Institut d'Études Nucléaires de l’Université d'Alger, 


Nous avons mis au point un appareil générateur 
d’impulsions lumineuses très brèves destiné à l’étude 
des photomultiplicateurs (1). Ces impulsions sont 
obtenues par simple transformation d’impulsions 
d'électrons par un écran luminescent. 


(:) Ce travail a été fait dans le cadre du contrat n° 1687 
du Commissariat à l'Énergie Atomique, 
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Le matériel mis en œuvre est le suivant : 


Le générateur d’impulsions lumineuses brèves utilise 
le spot obtenu sur un écran spécial de tube cathodique. 
à très faible rémanence par un faisceau électronique « 


pulsé commandé par des impulsions électriques de 


l’ordre d’une millimicroseconde de durée appliquées 
sur une grille de commande. Le spot obtenu est focalisé 
sur la surface de la photocathode du photomulti-” 


licateur en étude par l'intermédiaire d’un système 
P P ÿS 


optique comportant un diaphragme permettant de 
faire varier l'intensité de l’image. La surface éclairée 
est de l’ordre de 3.10-4 cm2. Le photomultiplicateur | 
est monté sur un chariot permettant de former l’image 
du spot en un poiñt quelconque de la photocathode | 
sans modifier les impulsions lumineuses dont les carac- « 


téristiques sont les suivantes : temps de montée 
de 10 % à 90 Y : 2 mus; durée à mi-hauteur : 12mys; 
variations d'amplitude inférieure à 1% ; variations de 
position inférieure à 0,2 mus ; longueur d’onde moyenne 
du rayonnement émis # 4000 À; fréquence de récur- 
rence des impulsions : 50 par seconde. 

Le signal électrique est observé à l’aide d’un oscil- 
lographe déclenché extérieurement par un signal syn- 
chrone de celui produisant l’allumage du flyimg-spot 
générateur d’impulsions lumineuses. 

La mesure des temps se fait en repérant la position 
du milieu de la montée de l’impulsion par rapport au 
début du balayage. 

Avec l’oscillographe Tektronix 517 A, cette mesure 
peut être faite à 0,2 mm près, la précision dans la 
mesure du temps est donc de 0,2 mus. 


1. Relevé de la carte des retards d’un photomul- 
tiplicateur en fonction du point d’impact sur la photo- 
cathode. — On appelle retard la différence de temps de 
transit entre un point donné de la photocathode et le 
centre. La mesure peut être faite à 0,5 mus près. La 
planche n° 1 représente les retards mesurés sur un 
tube 51 AVP. 


2. Mesure du temps de transit global dans un photo- 
multiplicateur. — Les retards étant différents suivant 
le point d'impact sur la photocathode cetre mesure est 
effectuée en n’illuminant que le centre de cette der- 
nière. Nous avons mesuré pour diverses valeurs de la 
tension d’alimentation U, les retards correspondant 
aux photomultiplicateurs types 51 AVP, 6 819 ot 


6 097 B et vérifié qu’ils variaient comme 1A/U. 


3. Mesure des fluctuations du temps de transit. — 
Les fluctuations du temps de transit sont directement 
observables à l’oscillographe, elles correspondent à des 
variations désordonnées de position de l'impulsion 
obtenue. Pour le photomultiplicateur type 51 AVP 
utilisé dans ses conditions normales, ces fluctuations 
sont inférieures à 0,5 mys. Elles deviennent plus impor- 
tantes pour des flux lumineux faibles ou pour des ten- 
sions accélératrices importantes, mais dans ces condi- 


tions l’observation est fortement perturbée par le bruit 


propre du photomultiplicateur. Les résultats montrent 
que les fluctuations sont fonction : du type de photo- 
multiplicateur, du potentiel accélérateur, du point 
d'impact sur la photocathode, du flux de photons reçus. 


| 


Fe 
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4. Mesure de l’allongement du temps de montée de 
l'impulsion électrique en fonction du flux lumineux 
incident. —, Nous avons constaté un allongement 


_ notable du temps de montée pour des valeurs crois- 


santes de l’intensité lumineuse incidente. Des essais 
systématiques ont été entrepris sur différents tubes et 


en divers endroits de leur photocathode. Pour le photo- 
* multiplicateur 51 AVP, nous avons obtenu un allon- 


gement atteignant 7 mus pour des flux lumineux 
variant dans un rapport de 16. 


re PT Ee SERRE 
Frc. 1, 


tiplicateur lorsqu'on déplace le point lumineux sur la 
photocathode. 

La figure 2 représente la carte des sensibilités du 
type 51 AVP, la fenêtre carrée de-collection ayant ses 
bords parallèles à ceux de la planche. Les différentes 
courbes d’atténuation sont chiffrées en décibels par 
rapport au point le plus sensible. : 


6. Conelusion. — Les résultats obtenus avec divers 
types de photomultiplicateurs tant en ce qui concerne 
les variations du temps de transit ou de sensibilité de la 
photocathode en fonction du point éclairé que les déter- 
minations absolues de temps de transit ou de leurs 
fluctuations sont utiles pour étudier la réponse du 
photomultiplicateur à une impulsion lumineuse inci- 
dente et les déformations qu'il apporte. 

Cette étude est actuellement poursuivie de façon à 
étudier séparément le comportement de l’ensemble 
photocathode-première dynode collectrice et le mul- 
tiplicateur lui-même. | 


Lettre reçue le 24 juin 1958. 
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Cet allongement du temps de montée serait compa- 
tible avec l'hypothèse d’après laquelle l’impulsion de 
sortie résulterait d’une superposition d’impulsions 
décalées dans le temps dont la dispersion pourrait peut- 
être dépendre d’une certaine défocalisation fonction 
de l’intensité. 


5. Carte des sensibilités. — Ces mesures ont été 
effectuées en compensant par un atténuateur étalonné 
les variations de sensibilité globale du photomul- 


FrG22: 


POLARISATION DE LA LUMIÈRE 
ENTRE 100 y ET 300 y 
PAR LES RÉSEAUX « ÉCHELETTES » 


Par MM. Armand Hapni, k 
Edmond Décamps et Pierre DELORME. 


Dans un spectromètre à réseau pour l’infrarouge 
lointain dont le principe a été donné par ailleurs [1], 
nous avons utilisé comme filtre un réseau « echelette » 
travaillant dans l’ordre zéro. Sa surface, carrée, permet 
d’orienter les traits de deux façons : 

19 parallèlement à ceux du réseau disperseur, 

20 perpendiculairement. 

Dans le premier cas, le récepteur ne reçoit presque 
plus d’énergie ; dans le second cas, la perte causée par 
le réseau filtre est très faible. 

On en déduit, non seulement que la lumière réfléchie 
par le réseau filtre est fortement polarisée, en bon 
accord avec des expériences effectuées avec des ondes 
de 3 em qui montrent que le vecteur vibration est 
perpendiculaire aux traits [2], mais encore que la 
lumière dispersée par le réseau principal dans le 
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premier ordre est aussi polarisée, mais ici parallèlement 
aux traits du réseau. 

Les théories concernant les réseaux « echelette » [3] 
ne tiennent pas compte de ce partage de la lumière 
dans les spectres d’ordre 0 et 1 du seul fait de la polari- 
sation et. ne semblent pas exactes dans le cas étudié. 
Par contre, pour les spectres-d’ordre élevé, nous avons 
vérifié avec un polaroïde que la lumière visible dis- 
persée dans le blaze n’était pas polarisée. L'effet de 
polarisation que nous signalons ne se produirait donc 
plus lorsque la longueur d’onde devient nettement infé- 
rieure à la largeur « a » des facettes. mais il est net dès 
que À > af2. (Ici, a = 250, af2 =1254) 

Ce sont précisément là les conditions dans lesquelles 
sont généralement employés les réseaux pour l’infra- 
rouge lointain. Il semble donc indispensable de tenir 
compte de cet effet qui présente un inconvénient : 
perte de près de la moitié de l’énergie, mais qui peut 
constituer un avantage dans l’étude des cristaux et de 
certaines matières plastiques. 


Lettre reçue le 18 juillet 1958. 
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EFFETS TRANSITOIRES 
DES DÉFORMATIONS MÉCANIQUES 
SUR LA CONDUCTIBILITÉ ÉLECTRIQUE 
DU FLUORURE DE LITHIUM 
ET DES HALOGÉNURES ALCALINS 


Par H. Curien et Z. MiHaILovic, 


Laboratoire de Minéralogie-Cristallographie 
de la Faculté des Sciences, Paris. 


Dans un précédent travail (Curien, Mihaiïlovic, 1957) 
nous avons montré comment la résistivité et l’aptitude 
à la polarisation électrique du fluorure de lithium 
étaient influencées par l’écrouissage. Nous avions pro- 
posé une explication faisant intervenir les ecrans 
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F1G. 1. — Une pointe d’acier peut venir piquer latéralement 
le cristal en 1 ou en 2. 
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chargés des dislocations. Nous nous proposons de 
décrire maintenant les phénomènes transitoires obser- 
vés dans les cristaux d’halogénures alcalins pendant la 
déformation. Des expériences de ce genre avaient déjà 
été réalisées par Gyulai et Hartly (1928), Stepanow 
(1933), Fischbach et Nowick (1955). Ces auteurs opé- 


raient à température ordinaire ou au plus égale | 


à 100 0C, températures auxquelles les déformations 
restent petites. Nous avons, au contraire, fait nos 
mesures à 600 0C. 

Un cristal de forme cubique de 1/2 cm de côté, 
obtenu par clivage, est placé dans un four maintenu 
entre deux électrodes d'acier inoxydable. Ce four est 
percé sur le côté d’uné ouverture par laquelle on peut 
introduire une pointe (aiguille de phonographe) qui 
vient s’enfoncer dans le cristal au voisinage de l’une 
ou l’autre des électrodes repérées par les chiffres 1 et 2 


(fig. 1). 


Observations sans tension extérieure. — Dans une 
première série d’expériences, l’échantillon cristallin 
placé dans le four est simplement court-cireuité à 
travers un galvanomètre. En l’absence de toute défor- 
mation, on constate l’existence d’un courant thermo- 
électrique très important. Ce courant correspond, dans 
le cas de LiF à 6000, à une différence de potentiel de 
l’ordre de 0,1 V, alors que les gradients de température 
dans le four, vérifiés à l’aide de couples, sont tels que 
les différences de température entre les divers points du 
cristal n’excèdent pas quelques degrés. Cet énorme 


pouvoir thermoélectrique dépend d’ailleurs de la modi- 


fication chimique des électrodes par les ions du cristal 
qui viennent s’y déposer du fait de la conductibilité 
ionique. Le pouvoir thermoélectrique se modifie en 
effet avec le temps. On lui redonne sa valeur initiale 
par élimination par clivage d’une mince couche au 
voisinage des électrodes. 

Si on pique le cristal en 1 ou en 2, on observe les 
phénomènes reproduits sur la figure 2 : au moment où 
on enfonce la pointe dans le cristal apparaît un brusque 
pic de courant dont la hauteur est de l’ordre de 10-$ A 
qui disparaît après quelques minutes. Le sens du cou- 
rant est tel qu’il sort par l’électrode la plus voisine du 
lieu de pénétration de la pointe. La piqüre laisse 
d’ailleurs dans le cristal une cavité de 2/10 mm de dia- 
mêtre et 4/10 de mm de profondeur environ. 


Observations avec tension extérieure appliquée. — , 


Reprenons maintenant les mêmes expériences, mais en 
plaçant dans le circuit extérieur une tension de l’ordre 
de 6 à 12 volts. Avec un échantillon de LiF, on obtient 
alors un courant permanent de l’ordre de 10-4 A. 
L’enfoncement de la pointe, en 1 ou en 2, fait diminuer 
ce courant comme le montre la figure 3. Dans tous les 
cas, quel que soit l’endroit de la piqüre ou le sens de la 
tension appliquée, c’est toujours une diminution tem- 


poraire du courant de conductibilité qui est observée. ? 


A 600 ©, le phénomène est plus-spectaculaire avec KCI 
(point de fusion 776 °C) qu'avec LiF (point de fusion 
870 °C). Pour un même effort sur la pointe, le cratère 
obtenu est d’ailleurs beaucoup plus important dans KCI 
(8/10 mm de diamètre et 1,8 mm de profondeur). 
D’autre part si, après pénétration, on laisse la pointe en 
place dans le cristal, le courant transitoire diminue 
lentement (temps de demi décroissance dans KCI 
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{minute 


116. 2. — Au courant thermoélectrique viennent se superposer des pics de courant provoqués par la pénétration de la 
pointe au voisinage de l’électrode 1 ou de l’électrode 2. (T = 6250, sans polarisation.) 


I 


1minute 


men) temps, 


T=625°C 
? -U=+6 Votts 


4 2 


F1c. 3. — Variation transitoire du courant de conduc- 
tibilité de LiF à 600° sous l’effet d’une piqüre. 


à 600 0C : = 25 min). Si on retire la pointe (portion b 
de la courbe de la figure 4) la décroissance est beaucoup 
plu s rapide (+ 1 min). | 


Interprétation. — Nous ne voulons ici qu’esquisser 
une interprétation que nous développerons dans une 
publication ultérieure détaillée. 

Le fluorure de lithium est un conducteur essentiel- 
lement cationique. En pénétrant dans le cristal, la 
pointe repousse la matière. Les lacunes de cations 


étant les défauts élémentaires de beaucoup les plus 
mobiles, les cations sont évacués au plus vite et le pic 
de courant positif observé partant de l’électrode la 
plus proche de la piqûre, dans les expériences sans 
tension extérieure, peut ainsi être expliqué. Cet effet 


t 


Fic. 4. — Variation transitoire du courant de conduc- 
tibilité de KCI à 6000. a : La pointe est brusquement 
enfoncée et on laisse en place dans le cristal, b: La 
pointe est retirée. 


est petit (de l’ordre de 106 A). Il est inobservable dans 
le cas des expériences avec tension extérieure. Dans ce 
cas, au contraire, circule dans le cristal un courant 
permanent important. La brusque déformation du 
cristal augmente sa résistivité, qui ne reprend sa 
valeur normale que lorsque la perturbation a disparu, 
c’est-à-dire après remise à l’équilibre thermique de 
échantillon. Cette remise à l’équilibre est beaucoup 
plus rapide lorsqu'on retire la pointe de la cavité qu’elle 
a imprimée dans le cristal. On constate d’ailleurs 
qu'après retrait de la pointe, le fond de la cavité 
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reprend un profil beaucoup plus arrondi que le profil 
de la pointe. 

Une étude complète dans LiF, NaCIl et KCI, en fonc- 
tion de la température, est en cours. 


Lettre reçue le 29 juillet 1958. 
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CALCUL DE L’INTENSITÉ 
DES BANDES INDUITES PAR LA PRESSION 
DANS LES MÉLANGES GAZEUX POLAIRES 
INFLUENCE DES MOMENTS QUADRUPOLAIRES 
ET DES FORCES DE DISPERSION 


Par LE. GALATRY, 
Laboratoire des Hautes Pressions, Bellevue (Seine-et-Oise). 


Dans les spectres infrarouges des mélanges com- 
primés d’un gaz polaire { avec un gaz non polaire 2, 
tous deux diatomiques, deux régions spectrales parti- 
culières sont le siège d’une absorption induite par 
‘l'augmentation de densité : celle de la fréquence fonda- 
mentale de vibration de la molécule 2 [1] et celle de la 
somme des fréquences fondamentales de vibration des 
molécules 1 et 2 [2]. Les coefficients d'absorption 
totaux À, et À, de ces bandes, réduits à un amagat 
de chaque gaz, ont été calculés à l’aide d’un modèle 
proposé antérieurement [3], en considérant que 


N° 10 


l'absorption étudiéé est due aux couples de molé- 
cules (1-2) présents dans le mélange. On a négligé 
l'influence de la quantification des degrés de liberté 
d'orientation, mais on a tenu compte de l’orientation 


préférentielle mutuelle des deux molécules déterminée 


par l’anisotropie des forces de dispersion. 

Les calculs ont été effectués dans le cas des mélanges 
HCI-H, et HBr-H, pour lesquels on possède des 
données expérimentales [4]. On a obtenu (en 
CMS am) 


HCLH, : À, — 21.10? 


(exp 82 70) 
A2 = 1,9.107# (exp. : 3,3.1074) 
HBr-H,: 41, —0,6.10 4\V(exp.: 12-1074); 


L'accord peut être considéré comme assez satis- 
faisant, eu égard aux approximations faites et à 
l’imprécision des valeurs adoptées pour les dérivées des 


moments quadrupolaires et des polarisabilités. Les À 


calculs montrent que l’effet des moments quadru- 
polaires peut contribuer au quart de la valeur des 
résultats finaux alors que les effets d’orientation 
n’entrent que pour le vingtième environ dans les 
mêmes résultats. 

Les détails de ces calculs paraîtront prochainement 
dans le Journal des Recherches du C. N. R.S. 


Lettre reçue le 30 juillet 1958. 
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